Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,287,687 C→T 5.1% intergenic (‑65/+95) 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,6870CT 5.1% 319.2 / 3.0 98intergenic (‑65/+95)tyrT/purUtRNA‑Tyr/formyltetrahydrofolate hydrolase
Reads supporting (aligned to +/- strand):  ref base C (53/40);  new base T (3/2);  total (56/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

TGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTG  >  NC_000913/1287541‑1287836
                                                                                                                                                  |                                                                                                                                                     
tggtggGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGAcgcg                                                                                                                                                    <  1:744060/150‑1 (MQ=255)
    gggggAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCGTCATTTGATATGATGcgcccc                                                                                                                                                <  1:1661499/150‑3 (MQ=21)
    gggggAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCGTCATTTGATATGATGcgcccc                                                                                                                                                <  2:498871/150‑3 (MQ=21)
    gggggAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCATTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGcgccgc                                                                                                                                                <  2:919017/150‑1 (MQ=255)
     ggggAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCt                                                                                                                                               >  1:1092113/1‑150 (MQ=255)
     ggggAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCt                                                                                                                                               >  2:1474444/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:166262/150‑1 (MQ=255)
                                 tGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAAc                                                                                                                   <  1:845591/150‑1 (MQ=255)
                                  gACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTATTTACAGGCCAGCACCCATATCGGGAAGCGGGTCGCATAATTACAATGTACGCGCCTCAGTAATCTGTTCCTTTGAGAAAAATaaaat                                                                                                                  >  2:1617161/1‑148 (MQ=16)
                                    cGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGaaaaa                                                                                                                        <  1:87208/142‑1 (MQ=255)
                                    cGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGaaaaa                                                                                                                        >  2:87208/1‑142 (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                                                                                    <  1:1319020/150‑1 (MQ=255)
                                                                ccGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTa                                                                                    >  2:1465781/1‑150 (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                                                                        >  1:10659/1‑150 (MQ=255)
                                                                            cccACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGt                                                                        <  2:1411625/150‑1 (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                        >  1:58728/4‑150 (MQ=255)
                                                                                                                            agcggCGCGTCATTTGATATGATGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACaaa                        >  2:517677/4‑150 (MQ=255)
                                                                                                                            agcggCGCGTCATTTGATATGATGCGCCGCTGTAAAGTGTTACGGTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACaaa                        >  1:138386/4‑150 (MQ=255)
                                                                                                                              gggCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAgc                      >  2:704883/1‑149 (MQ=255)
                                                                                                                               ggCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGAc                     >  1:806737/1‑150 (MQ=255)
                                                                                                                                        tatCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGCCGCGCTGTGc            >  1:1638671/1‑150 (MQ=255)
                                                                                                                                          tCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCa          >  2:648560/1‑150 (MQ=255)
                                                                                                                                          tCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCa          >  2:807448/1‑150 (MQ=255)
                                                                                                                                           cAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAg         >  2:773516/1‑150 (MQ=255)
                                                                                                                                            aaaTGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGt        >  2:1517593/1‑150 (MQ=255)
                                                                                                                                            aaaTGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGt        >  1:1303655/1‑150 (MQ=255)
                                                                                                                                               tGACGAGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTAATTTTTGCTGGAGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCAGTGACAGTACt     >  2:221198/1‑150 (MQ=255)
                                                                                                                                                  cgcgCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTg  <  1:1465781/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTG  >  NC_000913/1287541‑1287836

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: