Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 985,856 A→G 5.6% intergenic (‑147/+38) aspC ← / ← ompF aspartate aminotransferase/outer membrane porin F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913985,8560AG 5.6% 334.3 / 6.8 108intergenic (‑147/+38)aspC/ompFaspartate aminotransferase/outer membrane porin F
Reads supporting (aligned to +/- strand):  ref base A (47/54);  new base G (2/4);  total (49/59)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 109 total aligned reads displayed.

ATGACGAGGTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTT  >  NC_000913/985708‑986003
                                                                                                                                                    |                                                                                                                                                    
 tGACGAGGTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGt                                                                                                                                                   >  1:1655010/1‑150 (MQ=255)
    atagtttCCATTATTGTTACAGAAGGGAAGTCCGCTATAATGGTAATTTTAGATTTACATATTTCCAACTATTTCTGATTAAATTGAATTCTTTTTTTTCAAGATGTGTGATTTCTCTGGAATGTATAAAAAAAATAGGACAAAAGTCCt                                                                                                                                                <  2:1125062/145‑1 (MQ=255)
     gAGGTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTg                                                                                                                                               <  1:898410/150‑1 (MQ=255)
     gAGGTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTg                                                                                                                                               <  2:128945/150‑1 (MQ=255)
      aGGTTCCATTATGGTTACAGTAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGt                                                                                                                                              <  1:788898/150‑1 (MQ=255)
        gTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTcc                                                                                                                                                 <  2:518302/145‑1 (MQ=255)
        gTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTcc                                                                                                                                                 >  1:518302/1‑145 (MQ=255)
         ttCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGtttt                                                                                                                                           >  2:1619413/1‑150 (MQ=255)
          tCCATTATGGTTACAGAAGCGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGttttt                                                                                                                                          >  1:375728/1‑150 (MQ=255)
           ccATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGtttttt                                                                                                                                         >  1:1223276/1‑150 (MQ=255)
            cATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGCCGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTc                                                                                                                                        >  1:1838323/1‑150 (MQ=255)
            cATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACTTTGGTCCTGttttttt                                                                                                                                        <  2:1822899/150‑2 (MQ=255)
             aTTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACTCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAATGTATAAAAAAAACAGGACTTTGGTCCTGtttttttt                                                                                                                                       <  2:805866/150‑3 (MQ=255)
             aTTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACTTTGGTCCTGtttttttt                                                                                                                                       <  1:844619/150‑3 (MQ=255)
             aTTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCg                                                                                                                                       <  2:1095359/150‑1 (MQ=255)
             aTTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCg                                                                                                                                       <  2:1457427/150‑1 (MQ=255)
               tATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGc                                                                                                                                     <  2:799859/150‑1 (MQ=255)
               tATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGc                                                                                                                                     <  2:1506452/150‑1 (MQ=255)
                aTGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGGTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTATTTCGGCa                                                                                                                                    >  2:1404603/1‑150 (MQ=255)
                aTGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCa                                                                                                                                    <  2:1946623/150‑1 (MQ=255)
                 tGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACATAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAAAGGGCCAAAGGGCGTGTTTTTTGGGctt                                                                                                                                   >  1:1881770/1‑148 (MQ=255)
                   gTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCAttt                                                                                                                                 >  1:497730/1‑150 (MQ=255)
                      aCAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAAc                                                                                                                              <  2:1193225/150‑1 (MQ=255)
                         gAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACaaa                                                                                                                           >  1:432607/1‑150 (MQ=255)
                            gggAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAagag                                                                                                                        <  2:1351736/150‑1 (MQ=255)
                             ggAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAgg                                                                                                                       <  1:1461885/150‑1 (MQ=255)
                              gAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGgt                                                                                                                      >  2:1338/1‑150 (MQ=255)
                               aaGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAATACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGgtg                                                                                                                     <  2:375728/150‑1 (MQ=255)
                                 gTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGgt                                                                                                                      >  1:1650298/1‑147 (MQ=255)
                                 gTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGtgtg                                                                                                                   >  1:1821850/1‑150 (MQ=255)
                                 gTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACAAAAGTCATGTTTTTTCGGAATTTAAAAAAGAGGtgtg                                                                                                                   >  2:1225310/1‑150 (MQ=255)
                                   ccGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCt                                                                                                                 <  2:1675823/150‑1 (MQ=255)
                                   ccGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCt                                                                                                                 >  1:194570/1‑150 (MQ=255)
                                        aTCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAg                                                                                                            <  2:504860/150‑1 (MQ=255)
                                        aTCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAg                                                                                                            <  2:1390017/150‑1 (MQ=255)
                                        aTCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAg                                                                                                            <  2:497730/150‑1 (MQ=255)
                                          cAGTGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGaa                                                                                                          >  2:169688/1‑150 (MQ=255)
                                            gggTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACt                                                                                                        <  1:1619413/150‑1 (MQ=255)
                                              gTAACTGGAGTTTTACAAAATTCCTACTATTACTTATGAAATCGCAGGGTGTTTTTGCAAGACTTGAGTTTGCTCTTGAAGGTATAAAAAATACAGGACCAATGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTTCTATTAGAACTgg                                                                                                      <  2:650787/150‑1 (MQ=255)
                                               tAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGt                                                                                                     >  1:114242/1‑150 (MQ=255)
                                                aaCGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTa                                                                                                    <  2:62866/150‑1 (MQ=255)
                                                aaCGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTa                                                                                                    <  2:1049722/150‑1 (MQ=255)
                                                 aCGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTaa                                                                                                   >  2:1573029/1‑150 (MQ=255)
                                                   gggAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAAc                                                                                                 <  1:880786/150‑1 (MQ=255)
                                                     gagaTTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGa                                                                                               >  1:1279417/1‑150 (MQ=255)
                                                      agaTTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGAt                                                                                              >  1:1632645/1‑150 (MQ=255)
                                                       gaTTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAGGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATa                                                                                             >  2:1486417/1‑150 (MQ=255)
                                                       gaTTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATa                                                                                             <  2:1655010/150‑1 (MQ=255)
                                                           tACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTa                                                                                                                                >  2:1649600/1‑111 (MQ=255)
                                                           tACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTa                                                                                                                                <  1:1649600/111‑1 (MQ=255)
                                                             cAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCcaca                                                                                       >  1:68664/1‑150 (MQ=255)
                                                                      aaCTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGgtg                                                                              <  2:970442/150‑1 (MQ=255)
                                                                           ttACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTgtcgtc                                                                         <  2:737609/150‑1 (MQ=255)
                                                                                  tgaAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCt                                                                  >  1:316448/1‑150 (MQ=255)
                                                                                   gaAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTa                                                                 >  2:1422588/1‑150 (MQ=255)
                                                                                     aaaCGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACg                                                               <  1:1573029/150‑1 (MQ=255)
                                                                                     aaaCGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACg                                                               <  1:1877143/150‑1 (MQ=255)
                                                                                      aaCGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGc                                                              <  2:316448/150‑1 (MQ=255)
                                                                                           aGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGt                                                         >  1:577451/1‑150 (MQ=255)
                                                                                              cTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTttg                                                      >  1:554815/1‑150 (MQ=255)
                                                                                                 tttttGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTttgttg                                                   >  1:224284/1‑150 (MQ=255)
                                                                                                    ttGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCa                                                <  2:1998618/150‑1 (MQ=255)
                                                                                                     tGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAg                                               >  2:207832/1‑150 (MQ=255)
                                                                                                     tGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGCGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAg                                               <  2:320648/150‑1 (MQ=255)
                                                                                                       cAAGACGTGAGATTGCTCTGGAGGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGaa                                             <  1:1486417/150‑1 (MQ=255)
                                                                                                         aGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATc                                           <  1:1422588/150‑1 (MQ=255)
                                                                                                         aGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGATATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATc                                           >  1:1896823/1‑150 (MQ=255)
                                                                                                         aGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGATATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATc                                           >  1:1896827/1‑150 (MQ=255)
                                                                                                             gTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATc                                       <  1:69704/150‑1 (MQ=255)
                                                                                                             gTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATc                                       <  2:761602/150‑1 (MQ=255)
                                                                                                              tGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTttgttg                                                   >  2:529700/1‑137 (MQ=255)
                                                                                                              tGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTttgttg                                                   <  1:529700/137‑1 (MQ=255)
                                                                                                              tGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCt                                      <  1:1565581/150‑1 (MQ=255)
                                                                                                               gagaTTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTg                                     <  1:1198033/150‑1 (MQ=255)
                                                                                                               gagaTTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTg                                     >  2:1784659/1‑150 (MQ=255)
                                                                                                                    tGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTTGTCTGAACCTACGCCGAGTTTGTTGTCAGAATCGATCTGGTtga                                <  1:1406103/150‑1 (MQ=255)
                                                                                                                      ctctGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTtgatg                              >  2:679453/1‑150 (MQ=255)
                                                                                                                       tctGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTtgatga                             <  1:1308349/150‑1 (MQ=255)
                                                                                                                       tctGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTtgatga                             >  2:1308349/1‑150 (MQ=255)
                                                                                                                        ctGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTgatgat                            >  2:768319/1‑150 (MQ=255)
                                                                                                                         tGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGatgatg                           <  1:1155735/150‑1 (MQ=255)
                                                                                                                         gggAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGatgatg                           <  2:1914350/149‑1 (MQ=255)
                                                                                                                          ggAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCGCAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGt                          <  2:194570/150‑1 (MQ=255)
                                                                                                                          ggAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGt                          <  2:577451/150‑1 (MQ=255)
                                                                                                                          ggAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGt                          <  2:1632645/150‑1 (MQ=255)
                                                                                                                             aGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGt                       <  2:1465949/150‑1 (MQ=255)
                                                                                                                              ggTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTc                      <  2:68664/150‑1 (MQ=255)
                                                                                                                               gTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCa                     <  1:1404603/150‑1 (MQ=255)
                                                                                                                                tataAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCaa                    >  1:1091378/1‑150 (MQ=255)
                                                                                                                                  taAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACa                  >  1:1365887/1‑150 (MQ=255)
                                                                                                                                   aaaaaaaaCAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACAt                 >  2:826335/1‑150 (MQ=255)
                                                                                                                                      aaaaaCAGGACTTTGGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAgg              >  1:1999501/1‑150 (MQ=255)
                                                                                                                                      aaaaaCAGGACTTTGGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAgg              >  2:1083033/1‑150 (MQ=255)
                                                                                                                                      aaaaaCAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAgg              <  1:679453/150‑1 (MQ=255)
                                                                                                                                        aaaCAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTg            <  2:1279417/150‑1 (MQ=255)
                                                                                                                                         aaCAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTttgttg                                                   <  1:1650873/110‑1 (MQ=255)
                                                                                                                                         aaCAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTttgttg                                                   >  2:1650873/1‑110 (MQ=255)
                                                                                                                                             ggAACAAAGTCCTGTTTTTTAGGAATTTAAAAAAGAGTTGTTCTATTAGAACTGGTAAACGATAAAAACAGCAACGGTGTCGTCTGCACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACAt       >  2:1766865/1‑150 (MQ=255)
                                                                                                                                                 cAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGttt   >  2:1654322/1‑150 (MQ=255)
                                                                                                                                                  aaaGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGtttt  >  2:389561/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ATGACGAGGTTCCATTATGGTTACAGAAGGGAAGTCCGCTATCAGGGTAACGGGAGATTTACAAAATTCCAACTATTACTGATGAAAACGCAGGCTGTTTTTGCAAGACGTGAGATTGCTCTGGAAGGTATAAAAAAAACAGGACCAAAGTCCTGTTTTTTCGGCATTTAACAAAGAGGTGTGCTATTAGAACTGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTT  >  NC_000913/985708‑986003

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