Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,287,673 A→G 5.5% intergenic (‑51/+109) tyrT ← / ← purU tRNA‑Tyr/formyltetrahydrofolate hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,287,6730AG 5.5% 269.1 / 4.9 94intergenic (‑51/+109)tyrT/purUtRNA‑Tyr/formyltetrahydrofolate hydrolase
Reads supporting (aligned to +/- strand):  ref base A (47/40);  new base G (3/2);  total (50/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.46e-01

TGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGAC  >  NC_000913/1287526‑1287817
                                                                                                                                                   |                                                                                                                                                
tgtgGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGcatc                                                                                                                                                <  1:1227033/150‑1 (MQ=255)
   ggtgAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGcgtcatt                                                                                                                                             >  1:498871/1‑149 (MQ=255)
    gtgAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCatat                                                                                                                                            >  1:1411625/1‑150 (MQ=255)
    gtgAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCatat                                                                                                                                            <  1:259315/150‑1 (MQ=255)
               tggtggGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGAcgcg                                                                                                                                 <  1:744060/150‑1 (MQ=255)
                   gggggAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCGTCATTTGATATGATGcgcccc                                                                                                                             <  2:498871/150‑3 (MQ=21)
                   gggggAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCGTCATTTGATATGATGcgcccc                                                                                                                             <  1:1661499/150‑3 (MQ=21)
                   gggggAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCATTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGcgccgc                                                                                                                             <  2:919017/150‑1 (MQ=255)
                    ggggAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCt                                                                                                                            >  2:1474444/1‑150 (MQ=255)
                    ggggAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCt                                                                                                                            >  1:1092113/1‑150 (MQ=255)
                        aaGGATTCGAACCTTCGAAGTCGATGACGGTAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTaa                                                                                                                        <  1:177595/150‑1 (MQ=255)
                         aGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTaaa                                                                                                                       >  1:412731/1‑150 (MQ=255)
                          ggATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCCCCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGCGGCGCATCATATCAAATGACGCGCCGCTGTAAAg                                                                                                                      >  2:382770/1‑150 (MQ=255)
                          ggATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAg                                                                                                                      >  1:608822/1‑150 (MQ=255)
                          ggATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCGTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAg                                                                                                                      >  2:848747/1‑150 (MQ=255)
                           gATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAgt                                                                                                                     >  2:1291500/1‑150 (MQ=255)
                            aTTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTTCTGGCCTGCTACCTTATCGGGAAGCGGGGCGCATAATAACAAATGACGCGCCGCAGTAAAgtg                                                                                                                    >  2:694325/1‑150 (MQ=25)
                             ttCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAgtgt                                                                                                                   <  1:1143448/150‑1 (MQ=255)
                               cGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTa                                                                                                                 >  1:365148/1‑150 (MQ=255)
                                gAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATAATATAAAAAGACGCGCCGCTGTAAAGTGTAAc                                                                                                                >  2:189511/1‑150 (MQ=255)
                                  aCCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGt                                                                                                              >  1:1118299/1‑150 (MQ=255)
                                    cTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTg                                                                                                            >  1:643051/1‑150 (MQ=255)
                                     ttCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTga                                                                                                           <  2:325291/150‑1 (MQ=255)
                                     ttCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTga                                                                                                           >  1:1315561/1‑150 (MQ=255)
                                        gAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGaa                                                                                                        <  1:1051483/150‑1 (MQ=255)
                                             cGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATTCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATg                                                                                                   <  1:1125741/150‑1 (MQ=255)
                                               aTGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCGGTTGACTGGCCTGCTCCCGTATCGGGAAGCGGGGCGCATGAGATCACAGGACGCCCCGCTGGAAATTGTTAAGTTGAGAAAAATgca                                                                                                 >  2:271842/1‑148 (MQ=25)
                                                tGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAAc                                                                                                <  1:845591/150‑1 (MQ=255)
                                                tGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAAc                                                                                                <  1:166262/150‑1 (MQ=255)
                                                 gACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTATTTACAGGCCAGCACCCATATCGGGAAGCGGGTCGCATAATTACAATGTACGCGCCTCAGTAATCTGTTCCTTTGAGAAAAATaaaat                                                                                               >  2:1617161/1‑148 (MQ=16)
                                                   cGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGaaaaa                                                                                                     >  2:87208/1‑142 (MQ=255)
                                                   cGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGaaaaa                                                                                                     <  1:87208/142‑1 (MQ=255)
                                                   cGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTgg                                                                                             >  2:1308242/1‑150 (MQ=255)
                                                     gCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGtt                                                                                           >  1:1700537/1‑150 (MQ=255)
                                                        gATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCg                                                                                        <  2:804108/150‑1 (MQ=255)
                                                            tACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAAt                                                                                    >  1:987531/1‑150 (MQ=255)
                                                             aCAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAAtt                                                                                   >  1:1608276/1‑150 (MQ=255)
                                                               aGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAAtttt                                                                                 <  1:702665/150‑1 (MQ=255)
                                                                 tCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCa                                                                               <  1:1366822/150‑1 (MQ=255)
                                                                     cTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCg                                                                           >  1:936517/1‑150 (MQ=255)
                                                                     cTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTCGCGTAATTTTCATCCg                                                                           <  2:425078/150‑1 (MQ=255)
                                                                         cTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAAc                                                                       >  2:586625/1‑150 (MQ=255)
                                                                         cTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAAc                                                                       >  2:1400984/1‑150 (MQ=255)
                                                                         cTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAAc                                                                       <  2:412731/150‑1 (MQ=255)
                                                                          tttGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACg                                                                      <  1:189511/150‑1 (MQ=255)
                                                                              gCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGAtt                                                                  >  2:392677/1‑150 (MQ=255)
                                                                               ccGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTa                                                                 >  2:1465781/1‑150 (MQ=255)
                                                                               ccGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTa                                                                 <  1:1319020/150‑1 (MQ=255)
                                                                                cGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGTTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTAGAGAAAAATGAAATGGTTGCGTAATTTACATCCGTAACGGATTaa                                                                >  1:690475/1‑150 (MQ=255)
                                                                                cGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTaa                                                                >  2:1493429/1‑150 (MQ=255)
                                                                                 gCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTaaa                                                               <  1:842325/150‑1 (MQ=255)
                                                                                    tgggAACCCCACCAGGGTGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGATGCGGGGCGCATCTTATCAATTGTCGCGCCGCTGTAAAGTGTTTCGTTGTGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGt                                                            <  1:815698/149‑1 (MQ=255)
                                                                                    cGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGt                                                            >  2:886224/1‑150 (MQ=255)
                                                                                       gAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAAc                                                         <  2:802771/150‑1 (MQ=255)
                                                                                          ccccACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCat                                                      >  2:1179587/1‑149 (MQ=255)
                                                                                          ccccACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAg                                                      <  2:199014/150‑1 (MQ=255)
                                                                                          ccccACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAg                                                      >  2:1546962/1‑150 (MQ=255)
                                                                                           cccACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGt                                                     <  2:1411625/150‑1 (MQ=255)
                                                                                           cccACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGt                                                     >  1:10659/1‑150 (MQ=255)
                                                                                             caccacGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTa                                                   <  2:936517/150‑1 (MQ=255)
                                                                                               ccacGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTAtt                                                 >  1:714446/1‑150 (MQ=255)
                                                                                                  cGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTAttttt                                              >  2:131101/1‑150 (MQ=255)
                                                                                                  cGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTAttttt                                              <  1:1670070/150‑1 (MQ=255)
                                                                                                   ggggTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTg                                             <  2:608822/150‑1 (MQ=255)
                                                                                                     ggTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTTTTTTTGCt                                           <  2:696258/150‑1 (MQ=255)
                                                                                                      gTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTg                                          <  1:1546962/150‑1 (MQ=255)
                                                                                                      gTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTg                                          <  2:564949/150‑1 (MQ=255)
                                                                                                        aaTGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGc                                        <  1:1541551/150‑1 (MQ=255)
                                                                                                         aTGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCg                                       <  2:1296589/150‑1 (MQ=255)
                                                                                                         aTGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCg                                       <  2:1118299/150‑1 (MQ=255)
                                                                                                         aTGCTATTTACTGGCCTTCTCCCTTATCGGGCAGCGGGGCGCTTCATATCAAATGTCGCGCCGCTGTAAAGTGTTAAGTTGAGAAAAATGACCAGGTTGCGTAATTTTCATCCGTACCGGATAAAAGGTAACCAGTTATTTTTGCTGGCg                                       <  2:450587/150‑1 (MQ=255)
                                                                                                             tttttACTTGCCTTCTCAGTTTTATGGAATCGGGGAGCATTATATCAAATGACGAGACGCTGTAAAGTGTTACGTTGAGAAATATGAACTGTTTGCGTAATTTTCATCCTTAACGGATTAAAGGTAACGAGTTATTTTTGCTGGCGATTa                                   <  1:24773/150‑1 (MQ=255)
                                                                                                             tttttACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTa                                   <  1:694325/150‑1 (MQ=255)
                                                                                                                 tACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTTTTTTTGCTGGCGATTAAAGa                               >  1:620208/1‑150 (MQ=255)
                                                                                                                 tACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGa                               >  1:307832/1‑150 (MQ=255)
                                                                                                                  aCTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGaa                              >  2:1034950/1‑150 (MQ=255)
                                                                                                                     ggCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGaataa                           >  1:208579/1‑150 (MQ=255)
                                                                                                                      gCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATGAAATCACGCGCAGCTGTAAAGAGATAAGTTGATAAAAAAGAACTGGTTGCGTAATTTTCATCAGTATCGGATTAAAGTATACCAGTTATATTTCCTCGAGATTAAAGaaaaat                          >  2:160795/1‑150 (MQ=255)
                                                                                                                         tGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGAGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGt                       >  2:1212055/1‑150 (MQ=255)
                                                                                                                          gCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGtt                      >  2:65127/1‑150 (MQ=255)
                                                                                                                           cTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTc                     <  2:987531/150‑1 (MQ=255)
                                                                                                                           cTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTc                     <  1:1617161/150‑1 (MQ=255)
                                                                                                                            tCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCg                    >  1:8108/1‑150 (MQ=255)
                                                                                                                             cccTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGa                   >  2:67588/1‑150 (MQ=255)
                                                                                                                               gttATCGGGTAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACTGATTAAAGGTAACCAGTTATTTTTGCTGGGGATTAAAGAATAATCGTTCGAtt                 <  1:1172813/149‑1 (MQ=255)
                                                                                                                                ttATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTa                >  1:936121/1‑150 (MQ=255)
                                                                                                                                  aTCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTAcc              <  1:1224161/150‑1 (MQ=255)
                                                                                                                                   tCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCg             <  1:1493429/150‑1 (MQ=255)
                                                                                                                                    cgtgAAGCGGGGCGATTCATATGATATGACGCGCCGCTGTAAAGTGTTACGTTGAGTATAATGAACTGGTTTCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGTTTAAAGAATAATCGTTCGATTACCGt            <  1:154326/147‑1 (MQ=255)
                                                                                                                                        aaGCGGGGCGCATCATATCAAATGACGCTCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAAAAAACCAGTTATTATTGATGGCGATTTAAGAATAATCGTTCGATTACCGTAAAc        >  2:264878/1‑150 (MQ=255)
                                                                                                                                           cGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACaaa     <  1:1047178/150‑1 (MQ=255)
                                                                                                                                           agcggCGCGTCATTTGATATGATGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACaaa     >  2:517677/4‑150 (MQ=255)
                                                                                                                                           agcggCGCGTCATTTGATATGATGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACaaa     >  1:58728/4‑150 (MQ=255)
                                                                                                                                           agcggCGCGTCATTTGATATGATGCGCCGCTGTAAAGTGTTACGGTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACaaa     >  1:138386/4‑150 (MQ=255)
                                                                                                                                             gggCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAg  >  2:704883/1‑149 (MQ=255)
                                                                                                                                              ggCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGAc  >  1:806737/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGAC  >  NC_000913/1287526‑1287817

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