Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,719,429 G→A intergenic (‑385/‑49) yiaF ← / → yiaG DUF3053 domain‑containing protein YiaF/putative DNA‑binding transcriptional regulator YiaG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,719,4290GA93.8% 132.1 / 0.9 49intergenic (‑385/‑49)yiaF/yiaGDUF3053 domain‑containing protein YiaF/putative DNA‑binding transcriptional regulator YiaG
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (23/22);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.91e-01

AATTTGTTGTTTCATTGTTAAAAAATGATAGTGCCGTTGCCTTTTTTGTAGCCAGAGATTTACGCCTTCGCAGAAGGTTAGATTTTCGAAATCTTTACAAAACAGGAGTGGTAAGAGATATATCGGCATTTATCAGATCGTTAGGTGGCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATGAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTGGAACAGATTGTTTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGTATCAGTCGCCATGGTAAAGGAAT  >  NC_000913/3719185‑3719673
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
aaTTTGTTGTTTCATTGTTAAAAAATGATAGTGCCGTTGCCTTTTTTGTAGCCAGAGATTTACGCCTTCGCAGAAGGTTAGATTTTCGAAATCTTTACAAAACAGGAGTGGTAAGAGATATATCGGCATTTATCAGATCGTTAGGTGGCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCCCCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTc                                                                                                                                                                                                                                                >  2:109333/1‑251 (MQ=255)
 aTTTGTTGTTTCATTGTTAAAAAATGATAGTGCCGTTGCCTTTTTTGTAGCCAGAGATTTACGCCTTCGCAGAAGGTTAGATTTTCGAAATCTTTATAAAACAGGAGTGGTAAGAGATATATCGGCATTTATCAGATCGTTAGGTGGCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCg                                                                                                                                                                                                                                               >  2:306844/1‑251 (MQ=255)
  tttGTTGTTTCATTGTTAAAAAATGATAGTGCCGTTGCCTTTTTTGTAGCCAGAGATTTACGCCTTCGCAGAAGGTTAGATTTTCGAAATCTTTACAAAACAGGAGTGGTAAGAGATATATCGGCATTTATCAGATCGTTAGGTGGCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCg                                                                                                                                                                                                                                               >  1:157757/1‑250 (MQ=255)
           tCATTGTTAAAAAATGATAGTGCCGTTGCCTTTTTTGTAGCCAGAGATTTACGCCTTCGCAGAATGTTAGATTTTCGAAATCTTTACAATACAGGAGTGGTAAGAGATATATCGGCATTTATCAGATCGTTAGGTAGCTATAAGTCACGTTAAATGACAGTCTGAATCGAATCATATCCAGAGCATGCCCTGACTTCCCCCCTCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTc                                                                                                                                                                                                                                      <  1:294415/250‑1 (MQ=255)
            cATTGTTAAAAAATGAGAGTGCCGTTGCCTTTTTTGTAGCCAGAGATTTACGCCTTCGCAGAAGGTTAGAGTTTCGAAATCTTTACAAAACAGGAGTGGTAAGAGATATATCGGCATTTATCAGATCGTTAGGTGGCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCGTCGAGGCAATTCa                                                                                                                                                                                                                                     <  2:248392/250‑1 (MQ=255)
             aTTGTTAAAAAATGATAGTGCCGTTGCCTTTTTTGTAGCCAGAGATTTACGACTTCGCAGAAGGTTAGATTTTCGAAATCTTTACAAAACAGGGGTGGTAAGAGATATATAGGCATTTATCCGATCGATAGGTGGCTATCAGTCACGCTAAATGACAGGCTGACTCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATCAGCTTCGATGCAATTCACg                                                                                                                                                                                                                                   >  1:211668/1‑251 (MQ=255)
                      aaaTGATAGTGCCGTTGCCTTTTTTGTAGCCAGAGATTTACGCCTTCGCAGAAGGTTAGATTTTCGAAATCTTTACAAAACAGGAGTGGTAAGAGATATATCGGCATTTATCAGATCGTTAGGTGGCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGc                                                                                                                                                                                                                            >  1:339560/1‑249 (MQ=255)
                                gCCGTTGCCTTTTTTGTAGCCAGAGATTTACGCCTTCGCAGAAGGTTAGATTTTCGAAATCTTTACAAAACAGGAGTGGTAAGAGATATATCGGCATTTATCAGATCGTTAGGTGGCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTga                                                                                                                                                                                                                <  1:397883/251‑1 (MQ=255)
                                  cGTTGCCTTTTTTGTAGCCAGAGATTTACGCCTTCGCAGAAGGTTAGATTTTCGAAATCTTTACAAAACAGGAGTGGTAAGAGATATATCGGCATTTATCAGATCGTTAGGTGGCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTgagg                                                                                                                                                                                                              >  1:223876/1‑251 (MQ=255)
                                     tGCCTTTTTTGTAGCCAGAGATTTACGCCTTCGCAGAAGGTTAGATTTTCGAAATCTTTACAAAACAGGAGTGGTAAGAGATATATCGGCATTTATCAGATCGTTAGGTGGCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGt                                                                                                                                                                                                           >  1:344057/1‑251 (MQ=255)
                                           ttttGTAGCCACATATTTACGCTTTCGCATAATGTTAGAGTTTCGAAATCTTTACAAAACAGGAGTGGTAGGAGATATATCGGCATTTATCAGAGCGTTAGGTGGCTATAAGGCATGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGCCGTCAGCCCGCTGTGGCTGCTTTTCCCGACTATGCTTAAGAAGCTGCGATGCAATGCGCTAGCGCGCAGTGTGATTTGAGGCGTTTTCa                                                                                                                                                                                                      <  2:211668/250‑1 (MQ=255)
                                            tttGTAGCCAGAGATTTACGCCTTCGCAGAAGGTTAGATTTTCGAAATCTTTACAAAACAGGAGTGGTAAGAGATATATCGGCATTTATCAGATCGTTAGGTGGCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCaa                                                                                                                                                                                                     <  1:90608/250‑1 (MQ=255)
                                               gTAGCCAGAGATTTACGCCTTCGCAGAAGGTTAGATTTTCGAAATCTTTACAAAACAGGAGTGGTAAGAGATATATCGGCATTTATCAGATCGTTAGGTGGCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATgg                                                                                                                                                                                                  <  2:223876/250‑1 (MQ=255)
                                                tAGCCAGAGATTTACGCCTTCGCAGAAGGTTAGATTTTCGAAATCTTTACAAAACAGGAGTGGTAAGAGATATATCGGCATTTATCAGATCGTTAGGTGGCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGa                                                                                                                                                                                                 >  1:450650/1‑250 (MQ=255)
                                                    cagcgaTTTACGCCTTCGCAGAAGGTTAGATTTTCGAAATCTTTACAAAACAGGAGTGGTAAGAGATATATCGGCATTTATCAGATCGTTAGGTTGCTATAAGTCACGCTAAAGGACAGGCTGAATCGAATCATAGCCAGAGCTTGCCCTGACTTCTCCCCGCTGTGTCGGCTTTTCCCTACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGTAAtata                                                                                                                                                                                            <  2:184220/247‑1 (MQ=255)
                                                        gaTTTACGCCTTCGCAGAAGGTTAGATTTTCGAAAACTTTACAAAACAGGAGTGGTAAGAGATATATCGGCATTTATCAGATCGTTAGGTGGCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGa                                                                                                                                                                                        <  2:52423/251‑1 (MQ=255)
                                                            tACGCCTTCGCAGACGGTTAGATTTTCGAAATCTTTACAAAACAGGCGTGGTAAGAGATATATCGGCATTTATCAGATCGTTAGGTGGCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCCGTGTGATTTGAGGAGTTTTCAATGGATTATAAAGATcc                                                                                                                                                                                     >  2:185665/1‑250 (MQ=255)
                                                             aCGCCTTCGCAGAAGGTTAGATTTTCGAAATCTTTACAAAACAGGAGTGGTAAGAGATATATCGGCATTTATCAGATCGTTAGGTGGCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCa                                                                                                                                                                                    >  1:397842/1‑250 (MQ=255)
                                                                    cgcggAAGGTTAGATTTTCGAAATCTTTACAAAACAGGAGTGGTAAGAGATATATCGGCATTTATCAGATCGTTAGGTGGCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATGAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATg                                                                                                                                                                             <  2:238086/246‑1 (MQ=255)
                                                                          aGGTTAGATTTTCGAAATCTTTACAAAACAGGAGTGGTAAGAGATATATCGGCATGTATCAGATCGTTAGGTGGCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCGGTGTCTGCTTTGCCCGACTATTCTTAATGAGCTTCGATGCAATTCACGATCCCGCAGTGTGATGTGAGGAGTGTTCAATGGAATATAAAGATCCAATGCATGAGCTGtt                                                                                                                                                                      <  2:244046/251‑1 (MQ=255)
                                                                           ggTTAGAGTTTCGAAATCGTTACAAAACAGGAGTGGTAAGAGATATATCGGCATTTATCAGATCGTTAGGTGGCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCTCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTg                                                                                                                                                                     <  2:35591/251‑1 (MQ=255)
                                                                                                aCAAAACAGGAGTGGTAAGAGATATATCGGCATTTATCAGATCGTTAGGTGGCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTGGAACAGATTGtt                                                                                                                                                >  1:11765/1‑251 (MQ=255)
                                                                                                     aCAGGAGTGGTAAGAGATATATCGGCATTTATCAGATCGTTAGGTGGCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTGGAACAGATTGTTTTTaa                                                                                                                                           <  1:358940/251‑1 (MQ=255)
                                                                                                      cAGGAGTGGTAAGAGATATATCGGCATTTATCAGATCGTTAGGTGGCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTGGAACAGATTGTTTTTaaa                                                                                                                                          <  1:496053/251‑1 (MQ=255)
                                                                                                                aaGAGATATATCGGCATTTATCAGATCGTTAGGTGGCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCATTTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGACTATAAGGATCCCATGCATGAGCTGTTGAGCAGCCTGGAACAGATTGTTTTTAAAGATGaaaagc                                                                                                                                >  1:79966/1‑248 (MQ=255)
                                                                                                                aaGAGATATATCGGCATTTATCAGATCGTTAGGTGGCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTGGCACAGATTGATTTTAAAGATGAAACGc                                                                                                                                >  2:33881/1‑251 (MQ=255)
                                                                                                                    gaTATATCGGCATTTATCAGATCGTTAGGTGGCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTTACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTGGAACAGATTGTTTTTAAAGATGAAACGCAGaa                                                                                                                            <  2:11765/251‑1 (MQ=255)
                                                                                                                           cGGCATTTATCAGATCGTTAGGTGGCTATAAGTCACGCTAAATGACAGGCTTAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTGGAACAGATTGTTTTTAAAGATGAAACGCAGAAAATTAc                                                                                                                      <  1:177120/250‑1 (MQ=255)
                                                                                                                               aTTTATCAGATCGTTAGGTGGCTATAAGTCACGCTAAATGACAGTCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACACCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTGGAACAGATTGTTTTTAAAGATGAAACGCAGAAAATTCCCCTg                                                                                                                  >  2:112459/1‑250 (MQ=255)
                                                                                                                                        aTCGTTAGGTGGCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTGGAACAGATTGTTTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGa                                                                                                         >  2:389233/1‑250 (MQ=255)
                                                                                                                                                gTGGCTATAAGTCACGCTAAATTACAGGCTGAATCGAATCATAGCCATAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGTCTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTGTAACAGATTGTTTCTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTcc                                                                                                <  2:378450/251‑1 (MQ=255)
                                                                                                                                                     tataAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAat                                                                                                                                                                                               >  2:455993/1‑151 (MQ=255)
                                                                                                                                                       taAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTGGAACAGATTGTTTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCg                                                                                         >  2:490003/1‑251 (MQ=255)
                                                                                                                                                                     tGACAGGCCGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTGGAACAGATTGTTTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGtt                                                                            >  1:349740/1‑250 (MQ=255)
                                                                                                                                                                      gACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTGGAACAGATTGTTTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTa                                                                           <  1:161047/250‑1 (MQ=255)
                                                                                                                                                                          ggCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTGGAACAGATTGTTTTTAAAGATGAAAc                                                                                                                                  <  1:81214/191‑1 (MQ=255)
                                                                                                                                                                          ggCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTGGAACAGATTGTTTTTAAAGATGAAAc                                                                                                                                  >  2:81214/1‑191 (MQ=255)
                                                                                                                                                                                  cgtagcATAGCCAAAGCATTCCCTGTCGTCACACCGCTGTTTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTGGAACAGATTGTTTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAg                                                              <  2:346533/246‑1 (MQ=255)
                                                                                                                                                                                  cgaatcATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATGAGCTTCGCTGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCATTGCATGAGCTGTTGAGCAGCCTGGAACAGATTGTTTTTAACGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCCGTTACGAAAAGGGAc                                                                >  1:204837/1‑249 (MQ=255)
                                                                                                                                                                                        aTAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTGGAACAGATTGTTTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTaa                                                        <  1:267718/251‑1 (MQ=255)
                                                                                                                                                                                                   aTGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTgg                                                                                                                                                           <  2:432317/141‑1 (MQ=255)
                                                                                                                                                                                                   aTGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTgg                                                                                                                                                           >  1:432317/1‑141 (MQ=255)
                                                                                                                                                                                                      cccTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTGGAACAGATTGTTTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTc                                          <  1:33881/251‑1 (MQ=255)
                                                                                                                                                                                                         tGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCCATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCCATGAATGAGCTGTTGAGCAGCATGGCACAGATTGTTTTTAAAGATGAAACGCAGAAACTTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAACATCGATGATTTCGcc                                       >  1:27052/1‑251 (MQ=255)
                                                                                                                                                                                                           aCTTCCCCCCGCTGGGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTGGAACAGATTGTTTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCg                                     <  2:450650/251‑1 (MQ=255)
                                                                                                                                                                                                                aCCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTGGAACAGATTGTTTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGttt                                >  1:105078/1‑251 (MQ=255)
                                                                                                                                                                                                                          gtCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTGGAACAGATTGTTTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGTATc                      <  1:306844/251‑1 (MQ=255)
                                                                                                                                                                                                                          gtCTGCTTTTCCCGACTATTCTTAATAAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTGGAACAGATTGTTTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGTATc                      >  1:96493/1‑251 (MQ=255)
                                                                                                                                                                                                                                              cTTAATAAGCTTCGATGCAATTCACGATACCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTGGAACAGATTGTTTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAAAAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCTCCCGGGTTTGGGGCGTATCAGTAGCCATGGTAAAGGAAt  >  2:299129/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
AATTTGTTGTTTCATTGTTAAAAAATGATAGTGCCGTTGCCTTTTTTGTAGCCAGAGATTTACGCCTTCGCAGAAGGTTAGATTTTCGAAATCTTTACAAAACAGGAGTGGTAAGAGATATATCGGCATTTATCAGATCGTTAGGTGGCTATAAGTCACGCTAAATGACAGGCTGAATCGAATCATAGCCAGAGCATGCCCTGACTTCACCCCGCTGTGTCTGCTTTTCCCGACTATTCTTAATGAGCTTCGATGCAATTCACGATCCCGCAGTGTGATTTGAGGAGTTTTCAATGGAATATAAAGATCCAATGCATGAGCTGTTGAGCAGCCTGGAACAGATTGTTTTTAAAGATGAAACGCAGAAAATTACCCTGACGCACAGAACAACGTCCTGTACCGAAATTGAGCAGTTACGAAAAGGGACAGGATTAAAAATCGATGATTTCGCCCGGGTTTTGGGCGTATCAGTCGCCATGGTAAAGGAAT  >  NC_000913/3719185‑3719673

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: