Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 422,429 T→G 21.3% intergenic (+70/‑89) proY → / → malZ proline‑specific permease/maltodextrin glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913422,4290TG21.3% 148.9 / 5.5 87intergenic (+70/‑89)proY/malZproline‑specific permease/maltodextrin glucosidase
Reads supporting (aligned to +/- strand):  ref base T (13/53);  new base G (19/2);  total (32/55)
Fisher's exact test for biased strand distribution p-value = 6.97e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGC  >  NC_000913/422282‑422578
                                                                                                                                                   |                                                                                                                                                     
ttCGCGTGCATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCATGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGCTCCGGCTAATCCGAAACTTATCCCCCCACGgcc                                                                                                                                                     >  2:1524766/1‑150 (MQ=255)
 tCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGtccg                                                                                                                                                    >  2:376100/1‑149 (MQ=255)
 tCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGtccc                                                                                                                                                    >  2:436224/1‑149 (MQ=255)
   gcgTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGtcctcc                                                                                                                                                  <  1:1770118/150‑1 (MQ=255)
   gcgTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGtcctcc                                                                                                                                                  <  1:164984/150‑1 (MQ=255)
       ggATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTccc                                                                                                                                              <  1:874891/150‑1 (MQ=255)
          ttgttgTGCTGTTGATTGGCTGGATGTTTAACCGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTGCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCa                                                                                                                                           <  2:33379/150‑1 (MQ=255)
            gttgtGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCGCCCCCcaa                                                                                                                                         >  2:1062813/1‑150 (MQ=255)
             ttgtGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGGCGCCGGCGGCCGCCCCCCAat                                                                                                                                        >  1:1654210/1‑150 (MQ=255)
              tgtgCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAata                                                                                                                                       >  2:1574076/1‑150 (MQ=255)
               gtgCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAatat                                                                                                                                      <  2:401634/150‑1 (MQ=255)
               gtgCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAatat                                                                                                                                      <  2:1845276/150‑1 (MQ=255)
                  cTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTGCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCt                                                                                                                                   <  1:570086/150‑1 (MQ=255)
                  cTGTTGATTGGCTGGAAGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCTCCACCCCCAACTATcag                                                                                                                                   >  1:418831/1‑148 (MQ=255)
                    gTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCGTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTc                                                                                                                                 <  1:1697463/150‑1 (MQ=255)
                     ttGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCa                                                                                                                                <  2:919823/150‑1 (MQ=255)
                     ttGATTGGCTGGATGTTTAAACGCCGCCACGATCGGCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCa                                                                                                                                <  2:1627868/150‑1 (MQ=255)
                        aTTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCGCGGGCGGCTCCCGCAAATATACTGGAgat                                                                                                                             >  1:696490/1‑150 (MQ=255)
                        aTTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAgat                                                                                                                             <  2:1024026/150‑1 (MQ=255)
                              tGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAg                                                                                                                       <  1:1679341/150‑1 (MQ=255)
                              tGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAg                                                                                                                       <  1:1104145/150‑1 (MQ=255)
                                gATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTg                                                                                                                     <  2:588659/150‑1 (MQ=255)
                                 aTGTTTAAACGCCTCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGa                                                                                                                    <  2:1821415/150‑1 (MQ=255)
                                    tttAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGATGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCGCCGTCCTCCGCCCACAAAAATAATTCACAGGATGAGTGagac                                                                                                                 >  2:1195161/1‑147 (MQ=255)
                                    tttAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCGCGGCCTCCGCCCGCAAAGAGAATACAGAGGATGAGTGATcg                                                                                                                 >  1:492966/1‑149 (MQ=255)
                                      tAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTg                                                                                                               <  2:854357/150‑1 (MQ=255)
                                       aaaCGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCGTGCTCCCCCAAATATCCTTCAGCTGATGAGTgctacggc                                                                                                              >  2:1341624/1‑143 (MQ=255)
                                          cgccgcCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGGCGCCCCCAAAGATCCGTCCGATGATGAGGGTACCTGCagt                                                                                                           >  2:1807185/1‑148 (MQ=255)
                                              gccACGATCGTCAGCTGGCTGAAACCCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGGTATCCGGGCCCGGCCTCAGCCCCCCAATATAGTTCAGCTGATGAGTGATACTGCAGTAggg                                                                                                       >  2:112262/1‑149 (MQ=255)
                                              gccACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGc                                                                                                       <  1:689152/150‑1 (MQ=255)
                                              gccACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGc                                                                                                       <  2:62457/150‑1 (MQ=255)
                                              gaaaCGATCGTCAGCTGGCTGAAAAACAGTAATCCTTGCGTTCTGTAAGCAGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGc                                                                                                       <  2:1165542/147‑1 (MQ=255)
                                                 aCGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTAt                                                                                                    <  1:1472318/150‑1 (MQ=255)
                                                  cGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCGTCCCCCCAATAACATTCAGATTATGAGTGATCCTGCAGTTGGGTATg                                                                                                   >  1:1589058/1‑150 (MQ=255)
                                                    aTCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGc                                                                                                 <  2:1658175/150‑1 (MQ=255)
                                                        tCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATTGCAAgg                                                                                             <  2:1116089/150‑1 (MQ=255)
                                                        tCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAgg                                                                                             <  2:412451/150‑1 (MQ=255)
                                                            ctggctgAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGAt                                                                                         <  2:16894/150‑1 (MQ=255)
                                                             tggctgAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATc                                                                                        <  2:1550317/150‑1 (MQ=255)
                                                                ctgAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCGCCCCCCAAAATACTTCAGATGAGGAGGGATCCGGCAGTAGGGGTTGGCAAGGGGaggcgg                                                                                     >  2:894164/1‑145 (MQ=255)
                                                                    aaaCCAGTAATCCTTGCGTTCTGTAAGCCGGTTCATGAGGAAATCCTGATCAGGCTATTCAGATAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTTATCAGAtttt                                                                                 <  2:635066/150‑1 (MQ=255)
                                                                     aaCCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAATGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTtc                                                                                <  1:427363/150‑1 (MQ=255)
                                                                      aCCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCGCGGCCTGCTCCACCAAATATACATAATATCATGTTTGATCCTGCATTTGGCTATGGCTTGTTGATCAGAGTTtca                                                                               >  2:532085/1‑150 (MQ=255)
                                                                      aCCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTtca                                                                               <  1:1718949/150‑1 (MQ=255)
                                                                      aCCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTtca                                                                               <  2:1498035/150‑1 (MQ=255)
                                                                      aCCAGTAATCCTTGCGTTCTGAAAGCCTGATCAGGAGGTAATCCTGATCCGGCTATTCCGAAAGTTATCCGCGCGCGTCCTCCTCCCACAATTATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTtca                                                                               <  1:116005/150‑1 (MQ=255)
                                                                        cAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCGCCCGGCCTCCGCCCCCCACTATCATTGAGATGTGGTGTGATGCTGGGTTAGGGTTTGGGATGGTGATGAGCttttgttg                                                                             >  2:752325/1‑146 (MQ=255)
                                                                            aaTCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAg                                                                         <  1:1471043/150‑1 (MQ=255)
                                                                              tCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAggg                                                                       <  2:27185/150‑1 (MQ=255)
                                                                               ccTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCGCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAgggg                                                                      >  2:1026019/1‑150 (MQ=255)
                                                                                cTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCGCCGGCCGCCGCCCCCACAAATCCATCAGAAGAGGAGTGATGCTGCATTAGGCTATGGCAAGGTGATCATATTTTAATCAGAggggt                                                                     >  2:1768682/1‑149 (MQ=255)
                                                                                     gTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGGCTCCCCCAAATATCCTTCAGATGATGTGTGATGCTGCATTAGGCTATGGGTAGGTGGTGAGATTTTCATCACAGGGGTGTTat                                                                >  2:1595815/1‑150 (MQ=255)
                                                                                     gTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCGCCCCCACATATCCTTCAGAGGATGAGGGATGGTGCGTTAGGCTTTGGGAAGGGGATGAGATTTTCATCACAGGGGAATTat                                                                >  2:1545449/1‑150 (MQ=255)
                                                                                     gTTCTGTAAGCCGGTTCAGGAGGAAATACTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGAGCAGATTTTCATCACAGGGGAAtttt                                                                <  1:1209532/150‑3 (MQ=255)
                                                                                      ttCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTatg                                                               <  2:660152/150‑1 (MQ=255)
                                                                                      ttCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTatg                                                               <  1:200363/150‑1 (MQ=255)
                                                                                       tCTGTAAGCCGGGTGAGGTGGAAATCCTGATCCGGCTATTCCGAAAGATATCCGGCCCCGGCCACCTCCCCCAAATATCCTTCACATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGATTTatga                                                              <  2:1445114/150‑1 (MQ=255)
                                                                                        cTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCGCCTCCCCCAAAAATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATTGCAAGATGATCAGATTTTCATAACAGGGGAATAattat                                                             >  2:1684729/1‑150 (MQ=255)
                                                                                         tGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCGCCACCCCCAAATATACTTCAGATGATGAGTGATCCTCCAGTATGATATGGCAAGGTGATCAGATTTTCATAACACGGTAATTatgttt                                                            >  2:1001660/1‑147 (MQ=255)
                                                                                             aGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCGGCCCCCACACAGCCTGAAGATGATGATGGCTCCAGCAGTAGGCTAGGGCTAGGTGATCAGAGTTTCATCACCGGGGAATTATGATGTTaa                                                        >  2:1413489/1‑150 (MQ=255)
                                                                                              gCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCGGCGGCCGCCGCCCCCAAATATGCTACAGATGATGAGTGAGCGTGCAGTAGGCTATGGGAAGGTGATCAGATTTTGATCAGAGGGGAATTATGATGTTtaa                                                       >  2:797969/1‑150 (MQ=255)
                                                                                                cGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCTGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATg                                                     >  2:669476/1‑150 (MQ=255)
                                                                                                   ttCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCAt                                                  <  2:1406312/150‑1 (MQ=255)
                                                                                                   ttCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCAt                                                  <  1:139960/150‑1 (MQ=255)
                                                                                                    tCAGGAGGAAATCCTGATCAGGCTATTCCGAAAGTTATCCGCCCCCTTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTATGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATg                                                 <  1:1344973/150‑1 (MQ=255)
                                                                                                            aaaTCCTGATCCGGCTATTCCGAAAGTTATCCGGCCCCGTCGTCCGCCACCAAATATCATGGAGATGATGAGGTATCCTGGGGTAGGGTTGGTCAAGGTGAGGGGTTTTTCTTAGCAGGGGTATTTTGATGTTTATTGCTTTGCagatac                                         >  2:351136/1‑145 (MQ=255)
                                                                                                            aaaTCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCGCCGACACCAAAGATCCGTCAGATTATGAGTGTTACTGCAGTAGGGGATGGCATGGTGATAGTAGTTAAATCGAAGTGGATTTTTGAGTTTAAATGACTGGCCCCTGc                                         >  2:356527/1‑150 (MQ=255)
                                                                                                             aaTCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGcc                                        <  2:345891/150‑1 (MQ=255)
                                                                                                                 cTGATCCGGCTATTCCGAAAGTTATCCGCCCGCGGCCGCCGCCCCCAAAAATCATTCAGATGATGAGTGGGCGTGCGTTAGGCTATGGGAAGGTGATGAGATTTTGATAACAGGGGAATTATGATGTTAAATTCAATGAACCTGCCGGTg                                    >  2:1791174/1‑150 (MQ=255)
                                                                                                                  tGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCGCCCCCAAAGATCCTTCAGATGAGGGGTGATCCGGCAGTAGGCTATGGCTAGGTGATCAGATTTTCATCACATGGGGATGATGATGTTAAATGCAAGGCACCTGCCGGGGc                                   >  2:430983/1‑150 (MQ=255)
                                                                                                                   gATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGcc                                  <  1:752325/150‑1 (MQ=255)
                                                                                                                    aTCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGCGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGccc                                 <  1:759583/150‑1 (MQ=255)
                                                                                                                     tCCGGCTATTCCGAAAGTTATCCGCCCCTGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGcccc                                <  2:1698509/150‑1 (MQ=255)
                                                                                                                     tCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGcccc                                <  2:1017353/150‑1 (MQ=255)
                                                                                                                      ccGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCGGCAGGTGccatg                               <  2:1838205/150‑4 (MQ=255)
                                                                                                                      ccGGCTATTCCGAAAGTTATCCGCCCCCGGCCTCCGCCCCCAACAAACCTACAGATGAAGAGTGATCCTGCATTAGGCTAAGGCAAGCGGATCAGAGATACATCACAGGGGAATTATGATGTTAAATGAATGGCAACAGCCGGTGccccc                               >  2:1100573/1‑150 (MQ=255)
                                                                                                                       cGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCGGCAGGTGccatgc                              <  1:1638244/150‑5 (MQ=255)
                                                                                                                       cGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGAGGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCGGCAGGTGccatgc                              <  2:170716/150‑5 (MQ=255)
                                                                                                                           tATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATCATGAGTGATCCTGCATTAGGCTACGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCACCATTTg                          <  1:913943/150‑1 (MQ=255)
                                                                                                                           tATTCCGAAAGTTATCCGCCCCCGGCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTg                          <  2:114831/150‑1 (MQ=255)
                                                                                                                            aTTCCGAAAGTTATCCGCCCCCGTCCTGCGCCCCCAAATAACCTTCAGATGATGAGTGATCCTGCATTAGGCGATGGGAAGGGGATCAGTTTTTGATCACAGGGGAATTATGATGTTTAATGGATGGACCCGGCCGGGGcccccctgtgg                         >  2:1261328/1‑144 (MQ=255)
                                                                                                                            aTTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGt                         <  2:85978/150‑1 (MQ=255)
                                                                                                                             ttCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGtt                        <  1:1265810/150‑1 (MQ=255)
                                                                                                                              tCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTa                       <  2:78064/150‑1 (MQ=255)
                                                                                                                               ccGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTaa                      <  2:1165474/150‑1 (MQ=255)
                                                                                                                                 gAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTGAGAAGATGAGGTATCCTTCATTTGGATATGTCAAGGTGATAAGATTTTAATCACAAGGTAATAATGATTTTAAATTAAAGGAACCTGCATGTGCCCCAAAATTTTAAAc                    >  2:1682290/1‑150 (MQ=255)
                                                                                                                                 gAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAACATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTTATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAAc                    <  1:1413489/150‑1 (MQ=255)
                                                                                                                                     gTTATCCGCCCCCGTCCTCCGCCCCCAAATAACCTAAAGATGATTAGTGATACTGGATTAGGCTATTGGAAGGTGATCAGATTTTCTTCACAGGGGAATTATGATGTTAATTTCATGGCAACTGACGGTGCGCCCAGTTGTTAAAcaaac                >  2:1573167/1‑149 (MQ=255)
                                                                                                                                           cGCCCCCGTCCTCCTCCCCCAAATTTCCTTAAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTGATCACAGGGGAATTTTGATTTTAAATTCATGGCACCTGACGTTGCCCCCATTTGTTAAAcaaagcaaata          >  2:90499/1‑148 (MQ=255)
                                                                                                                                           cGCCCACGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAATTGCATGGCACCTGCCGGTGCCCCATTTGTTAAACAAAGCAAAGAt         <  1:387097/150‑1 (MQ=255)
                                                                                                                                              ccccGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCa       <  1:603231/150‑1 (MQ=255)
                                                                                                                                                 cGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACt    <  1:376100/150‑1 (MQ=255)
                                                                                                                                                   tcctcctccCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGTGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCTCCCATTTGTTAAACAAAGCAAAGATCAACTGc  <  1:1524766/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGC  >  NC_000913/422282‑422578

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