Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,612,329 A→G 7.7% intergenic (+40/‑82) ytjA → / → yjjU uncharacterized protein/putative patatin‑like family phospholipase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,612,3290AG 7.7% 256.3 / 5.4 85intergenic (+40/‑82)ytjA/yjjUuncharacterized protein/putative patatin‑like family phospholipase
Reads supporting (aligned to +/- strand):  ref base A (35/41);  new base G (3/4);  total (39/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTC  >  NC_000913/4612188‑4612470
                                                                                                                                             |                                                                                                                                              
ggtggtCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACt                                                                                                                                        <  4:218640/149‑1 (MQ=255)
 gtggtCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTg                                                                                                                                       <  1:837780/149‑1 (MQ=255)
   ggtCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCC‑GCTGGACTGGc                                                                                                                                     <  1:101228/149‑1 (MQ=255)
     tCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCtt                                                                                                                                   <  1:404180/149‑1 (MQ=255)
     tCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCtt                                                                                                                                   >  2:560096/1‑149 (MQ=255)
      cTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGCGTTTGTTCATGGGCCCAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCGCT‑GGACTGGCtat                                                                                                                                  <  4:448794/149‑3 (MQ=255)
      cTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACGCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCGCT‑GGACTGGCttt                                                                                                                                  <  2:732964/149‑1 (MQ=255)
      cTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCGCT‑GGACTGGCttt                                                                                                                                  <  3:654986/149‑1 (MQ=255)
       tGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCtttt                                                                                                                                 <  3:109197/149‑1 (MQ=255)
         gCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGc                                                                                                                               <  4:356909/149‑1 (MQ=255)
            ggTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGt                                                                                                                            <  2:480181/149‑1 (MQ=255)
              tACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGttt                                                                                                                          >  4:339497/1‑149 (MQ=255)
                    tGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGa                                                                                                                    >  4:386509/1‑149 (MQ=255)
                    tGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGa                                                                                                                    <  1:888333/149‑1 (MQ=255)
                           gcAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAAtaataa                                                                                                             >  3:314993/1‑149 (MQ=255)
                            cAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATaaa                                                                                                            <  3:562639/149‑1 (MQ=255)
                             aGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAAt                                                                                                           >  2:674834/1‑149 (MQ=255)
                              gCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCttt                                                                                                                                  >  1:34010/1‑125 (MQ=255)
                              gCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCttt                                                                                                                                  <  2:34010/125‑1 (MQ=255)
                              gCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAAtt                                                                                                          <  1:10188/149‑1 (MQ=255)
                              gCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAAtt                                                                                                          <  3:230947/149‑1 (MQ=255)
                                tAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGc                                                                                                        <  1:721922/149‑1 (MQ=255)
                                tAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGc                                                                                                        >  4:503560/1‑149 (MQ=255)
                                tAAAATTGTCTTTGGCGTCGGGATTATTCTATTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGG‑GGACTGTATTATGCGGTTTGAGCGACTAatatatcac                                                                                                        <  2:222336/149‑4 (MQ=255)
                                 aaaaTTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCg                                                                                                       <  3:386509/149‑1 (MQ=255)
                                   aaTTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGtt                                                                                                     >  3:406307/1‑149 (MQ=255)
                                          tttGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGt                                                                                              <  4:202078/149‑1 (MQ=255)
                                           ttGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTc                                                                                             <  4:782125/149‑1 (MQ=255)
                                            tGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCa                                                                                            <  3:335777/149‑1 (MQ=255)
                                               cgtcgGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTc                                                                                         <  2:302768/149‑1 (MQ=255)
                                                gtcgGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCt                                                                                        >  2:213333/1‑149 (MQ=255)
                                                    ggATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTca                                                                                    >  1:337010/1‑149 (MQ=255)
                                                      attattCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTcaca                                                                                  >  3:385705/1‑149 (MQ=255)
                                                       ttattCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCacaa                                                                                 <  3:550400/149‑1 (MQ=255)
                                                               ttCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACa                                                                         <  1:481963/149‑1 (MQ=255)
                                                                  cTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGc                                                                                                        <  2:516174/115‑1 (MQ=255)
                                                                  cTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGc                                                                                                        >  1:516174/1‑115 (MQ=255)
                                                                   tGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACa                                                                         <  3:279764/145‑1 (MQ=255)
                                                                   tGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACa                                                                         >  4:279764/1‑145 (MQ=255)
                                                                   tGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGaa                                                                     <  1:674834/149‑1 (MQ=255)
                                                                   tGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑CGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGaa                                                                     <  3:461281/149‑1 (MQ=255)
                                                                     gTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGaa                                                                                                                   >  3:413738/1‑101 (MQ=255)
                                                                     gTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGaa                                                                                                                   <  4:413738/101‑1 (MQ=255)
                                                                     gTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAgg                                                                   >  4:205030/1‑149 (MQ=255)
                                                                     gTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAgg                                                                   >  3:525527/1‑149 (MQ=255)
                                                                     gTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAgg                                                                   >  1:400311/1‑149 (MQ=255)
                                                                          tttGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCagag                                                              <  2:400311/149‑1 (MQ=255)
                                                                          tttGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCagag                                                              <  1:213333/149‑1 (MQ=255)
                                                                             gTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTg                                                           >  1:372980/1‑149 (MQ=255)
                                                                              ttCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTgg                                                          <  3:277816/149‑1 (MQ=255)
                                                                                 aTGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGc                                                       >  1:873310/1‑149 (MQ=255)
                                                                                    taacGAAATCGACCCTAGACTTGAGAGCATTAAGAAAATTTTATCCAAAGCCAGAACAGC‑GGACTGAGTCTTGCGGGTCTAGCGAATAATAAACTACCTCACTTTGCGAATCGTTCACGACGTAATCAGTAAGGAACAGGTGGagcagc                                                    <  2:327513/146‑1 (MQ=255)
                                                                                        gAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGCGGGGGGGCAGCGAAt                                                >  2:452406/1‑149 (MQ=255)
                                                                                         aaaaCGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCATC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTTCTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATa                                               <  2:837466/149‑1 (MQ=255)
                                                                                           aaCGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATAcc                                             >  4:720519/1‑149 (MQ=255)
                                                                                                cccTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTaa                                        >  4:752289/1‑149 (MQ=255)
                                                                                                 ccTAGATTTCAGACCGTTAAGCAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGACTAATACCTTGCGTTCCTTTGTCCTTCTTTCCCTACAAAAACAGGACGGCAGAGGTGGGGCAGCGCCTACCTGTAAc                                       >  4:34771/1‑149 (MQ=255)
                                                                                                 ccTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAAc                                       <  3:809585/149‑1 (MQ=255)
                                                                                                  cTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACg                                      >  4:684084/1‑149 (MQ=255)
                                                                                                   tAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGc                                     >  2:152480/1‑149 (MQ=255)
                                                                                                    aGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAAc                                                                          <  1:616346/111‑1 (MQ=255)
                                                                                                    aGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAAc                                                                          >  2:616346/1‑111 (MQ=255)
                                                                                                    aGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCt                                    >  3:356861/1‑149 (MQ=255)
                                                                                                            gACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAata                            <  4:406307/149‑1 (MQ=255)
                                                                                                             aCCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGCAatat                           >  2:46075/1‑149 (MQ=255)
                                                                                                              ccATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATAtt                          <  1:152480/149‑1 (MQ=255)
                                                                                                              ccATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATAtt                          >  1:296275/1‑149 (MQ=255)
                                                                                                              ccATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATAtt                          <  1:46075/149‑1 (MQ=255)
                                                                                                              ccATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATAtt                          <  1:864661/149‑1 (MQ=255)
                                                                                                               cATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTg                         >  2:233513/1‑149 (MQ=255)
                                                                                                                aTTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTgc                        <  3:684084/149‑1 (MQ=255)
                                                                                                                   aaGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGcc                     <  1:233513/149‑1 (MQ=255)
                                                                                                                      aaaaTTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGtt                  <  3:5057/149‑1 (MQ=255)
                                                                                                                                ccAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAgg        <  4:525527/149‑1 (MQ=255)
                                                                                                                                ccAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAgg        >  3:101181/1‑149 (MQ=255)
                                                                                                                                 aaaaGCCAGTCCGCT‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATa                                               >  2:402024/2‑109 (MQ=255)
                                                                                                                                 aaaaGCCAGTCCGCT‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATa                                               <  1:402024/108‑1 (MQ=255)
                                                                                                                                  aaaGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGcc      >  4:12903/1‑149 (MQ=255)
                                                                                                                                  aaaGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGcc      >  2:500841/1‑149 (MQ=255)
                                                                                                                                  aaaGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAAAGCTTGGTAATATTGCGCCGTTGTCGCTAAGGcc      >  4:273540/1‑149 (MQ=255)
                                                                                                                                   aaGCCAGTCCGCT‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCg     >  2:497200/1‑149 (MQ=255)
                                                                                                                                   aaGCCAGTCCGCT‑GGACTGGCTTTTGCGGTTTTAGCGAAGAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGGTTGGTAATATTGCGCCGTTGTCGCTAAGGCCg     >  4:76286/1‑149 (MQ=255)
                                                                                                                                   aaGCCAGTCC‑GCTGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCg     >  3:863271/1‑149 (MQ=255)
                                                                                                                                     gCCAGTCTAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGtt   >  4:246082/1‑149 (MQ=255)
                                                                                                                                      ccAGTTCAGC‑GGTCTGTCTTTTGCGGTTTTAGCGAAGAATAAATTTCGTTACTTTGTCTTTCTATAACAACAAAATCAGGAAGGCAGAGGTGGGGCATAGTATACCTGTAACTCTTGGTAATATTGCTCCGTTGTCGCTAAGCCCGTTc  <  2:551582/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                              
GGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGC‑GGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTC  >  NC_000913/4612188‑4612470

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑