Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,894,852 T→G 13.5% intergenic (+81/+67) ygcP → / ← ygcQ putative anti‑terminator regulatory protein/putative flavoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,894,8520TG13.5% 340.7 / 46.5 133intergenic (+81/+67)ygcP/ygcQputative anti‑terminator regulatory protein/putative flavoprotein
Reads supporting (aligned to +/- strand):  ref base T (52/63);  new base G (8/10);  total (60/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.82e-01

Only 100 of 133 total aligned reads displayed.

CGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGT  >  NC_000913/2894705‑2895000
                                                                                                                                                   |                                                                                                                                                    
cgcgGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTg                                                                                                                                                    <  2:741738/150‑1 (MQ=255)
cgcgGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTg                                                                                                                                                    <  2:158588/150‑1 (MQ=255)
cgcgGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTg                                                                                                                                                    <  1:1309500/150‑1 (MQ=255)
  cgGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGAc                                                                                                                                                  <  1:1226431/150‑1 (MQ=255)
     gtgtCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGAcgcg                                                                                                                                               <  1:1270612/150‑1 (MQ=255)
       gtCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTc                                                                                                                                             >  1:366030/1‑150 (MQ=255)
        tCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCt                                                                                                                                            <  1:309582/150‑1 (MQ=255)
            ggCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATc                                                                                                                                        <  2:216900/150‑1 (MQ=255)
            ggCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCGTCAAGCGTcgcatc                                                                                                                                        <  1:442762/150‑6 (MQ=255)
               gcTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAgg                                                                                                                                     >  1:786112/1‑150 (MQ=255)
                 tttCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCCCGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGcc                                                                                                                                   >  2:1012381/1‑150 (MQ=255)
                  ttCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATAAGACGCGTCAAGCGTCGCATCAGGcgc                                                                                                                                  <  2:267594/150‑3 (MQ=255)
                   tgcaccacGAATACCGGGGTCTGGACTTTAGAGATTAAATTACTTTGACGGGATTATCGTTAGCTATTGAATTTGGTTTTATTTGTAGGCCGCACGCCAGATCGGACATTCAGCGCCTGATAAGACGCGTCAAGCGTCGCATCAGGcgct                                                                                                                                 <  2:541352/148‑4 (MQ=19)
                   tCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTa                                                                                                                                 <  1:550709/150‑1 (MQ=255)
                   tCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATAAGACGCGTCAAGCGTCGCATCAGGcgct                                                                                                                                 <  1:1299078/150‑4 (MQ=38)
                    ccaccaCGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTAc                                                                                                                                >  1:163805/1‑150 (MQ=255)
                    ccaccaCGAATACCGGGGTCTGGACCTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGcgctg                                                                                                                                <  2:1308388/150‑5 (MQ=38)
                     caccacGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACa                                                                                                                               >  1:6095/1‑150 (MQ=255)
                     caccacGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACa                                                                                                                               >  1:960211/1‑150 (MQ=255)
                     caccacGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACa                                                                                                                               <  2:1161088/150‑1 (MQ=255)
                      accacGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATAAGACGCGTCAAGCGTCGCATCAGGcgctgaa                                                                                                                              <  1:653205/150‑7 (MQ=38)
                       ccacGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAg                                                                                                                             <  2:326226/150‑1 (MQ=255)
                       ccacGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAg                                                                                                                             <  1:354705/150‑1 (MQ=255)
                         acGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTc                                                                                                                           >  2:937968/1‑150 (MQ=255)
                          cGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTcc                                                                                                                          >  2:1406070/1‑150 (MQ=255)
                               aCCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAgg                                                                                                                     <  2:710310/150‑1 (MQ=255)
                                ccGGGTGCTGGCCGCTTGCGACAAACTTACGTTGACTGGATACTCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGc                                                                                                                    <  1:668064/150‑1 (MQ=37)
                                ccGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGc                                                                                                                    >  2:700255/1‑150 (MQ=255)
                                  ggggTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCa                                                                                                                  <  2:964928/150‑1 (MQ=255)
                                          gACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTa                                                                                                          <  1:253306/150‑1 (MQ=255)
                                           aCGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAg                                                                                                         <  1:237073/150‑1 (MQ=255)
                                            cGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCGTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAgg                                                                                                        <  2:786112/150‑1 (MQ=255)
                                            cGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAgg                                                                                                        >  2:1473755/1‑150 (MQ=255)
                                               tAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAg                                                                                                     <  2:927145/150‑1 (MQ=255)
                                                aGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAgg                                                                                                    >  1:125332/1‑150 (MQ=255)
                                                aGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATAAGACGCGTCAAGCGTCGCATCAGGcgctgaa                                                                                                                              <  1:978606/124‑7 (MQ=37)
                                                 gCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGa                                                                                                   >  1:807966/1‑150 (MQ=255)
                                                  cGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGat                                                                                                  >  1:926961/1‑150 (MQ=255)
                                                   gAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGatc                                                                                                 >  1:106398/1‑150 (MQ=255)
                                                   gAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGatc                                                                                                 >  1:1273630/1‑150 (MQ=255)
                                                    aaaaaaTTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGatca                                                                                                >  2:376407/1‑150 (MQ=255)
                                                        aaTTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCg                                                                                            <  2:569431/150‑1 (MQ=255)
                                                            aCTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCaa                                                                                        >  2:234868/1‑150 (MQ=255)
                                                            aCTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATAAGACGCGTCAAGCGTCGCATCAGGcgctg                                                                                                                                <  1:720611/110‑5 (MQ=25)
                                                            aCTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATAAGACGCGTCAAGCGTCGCATCAGGcgctg                                                                                                                                >  2:720611/1‑106 (MQ=25)
                                                             cTTTGACGGGATACTCGTAAACAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGAGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAac                                                                                       >  2:1253850/1‑150 (MQ=255)
                                                             cTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAac                                                                                       >  2:1234095/1‑150 (MQ=255)
                                                              tttGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAaca                                                                                      <  1:126161/150‑1 (MQ=255)
                                                                   cGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCACGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTa                                                                                 >  1:1001727/1‑150 (MQ=255)
                                                                     ggATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTTCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTAtt                                                                               >  2:1189154/1‑150 (MQ=255)
                                                                      gATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTg                                                                              <  1:563513/150‑1 (MQ=255)
                                                                                  cAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTCCAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCACCCCGTTATTGGCAATCTGCGTg                                                                  >  2:77786/1‑150 (MQ=255)
                                                                                  cAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTg                                                                  <  2:729041/150‑1 (MQ=255)
                                                                                   aaTTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGa                                                                 >  1:45245/1‑150 (MQ=255)
                                                                                      tGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATg                                                              <  2:960211/150‑1 (MQ=255)
                                                                                        aaTTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGtt                                                            <  2:406646/150‑1 (MQ=255)
                                                                                        aaTTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGtt                                                            <  1:1234095/150‑1 (MQ=255)
                                                                                              gTTTGATTTGTAGGCCGAACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCAc                                                      >  1:1358719/1‑150 (MQ=255)
                                                                                               tttGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCAcc                                                     >  1:848178/1‑150 (MQ=255)
                                                                                               tttGATTTGTAGGCCGCACGCCACATCCGACATTAAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCAcc                                                     >  2:799429/1‑150 (MQ=255)
                                                                                                 tGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAg                                                   <  1:21138/150‑1 (MQ=255)
                                                                                                     ttGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGcc                                               >  2:1422119/1‑150 (MQ=255)
                                                                                                      tGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCt                                              <  2:1236807/150‑1 (MQ=255)
                                                                                                       gTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTc                                             <  2:1242311/150‑1 (MQ=255)
                                                                                                        tAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTcc                                            <  1:1063422/150‑1 (MQ=255)
                                                                                                        tAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTcc                                            <  1:234868/150‑1 (MQ=255)
                                                                                                          ggCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAg                                          >  1:1405723/1‑150 (MQ=255)
                                                                                                           gCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGc                                         <  2:1273630/150‑1 (MQ=255)
                                                                                                            ccGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGtt                                                            >  2:489200/1‑130 (MQ=255)
                                                                                                            ccGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGtt                                                            <  1:489200/130‑1 (MQ=255)
                                                                                                             cGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCGCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGcac                                       <  2:681008/150‑1 (MQ=255)
                                                                                                             cGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGcac                                       <  1:681013/150‑1 (MQ=255)
                                                                                                                   ccACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACttt                                 >  2:1260848/1‑150 (MQ=255)
                                                                                                                     acaTCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTcc                               <  1:490972/150‑1 (MQ=255)
                                                                                                                        tCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGt                            <  1:376407/150‑1 (MQ=255)
                                                                                                                        tCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGt                            >  2:1337818/1‑150 (MQ=255)
                                                                                                                         ccGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGtc                           <  2:658538/150‑1 (MQ=255)
                                                                                                                         ccGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTAATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGtc                           >  2:973563/1‑150 (MQ=255)
                                                                                                                          cGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGtca                          <  1:1222866/150‑1 (MQ=255)
                                                                                                                          cGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGtca                          >  2:1169659/1‑150 (MQ=255)
                                                                                                                             cATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGtcatca                       <  1:1473755/150‑1 (MQ=255)
                                                                                                                              cttgtaggcCTGATAAGACGCGTCAAGCGTCGCATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCaa                      >  2:462785/8‑150 (MQ=38)
                                                                                                                              aTTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCaa                      >  2:1087025/1‑150 (MQ=255)
                                                                                                                               ttCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAAc                     <  2:125332/150‑1 (MQ=255)
                                                                                                                                 cAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCa                   >  2:518947/1‑150 (MQ=255)
                                                                                                                                    cgcCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGc                <  1:1047492/150‑1 (MQ=255)
                                                                                                                                     gcCTGATAAGACGCGTCAAGCGTCGCATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGcc               >  2:1024642/1‑150 (MQ=255)
                                                                                                                                     gcCTGATAAGACGCGTCAAGCGTCGCATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGcc               >  2:600748/1‑150 (MQ=255)
                                                                                                                                      ccTGATAAGACGCGTCAAGCGTCGCATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCCCGCCa              >  2:1195374/1‑150 (MQ=255)
                                                                                                                                        tGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGACCCGTAGGCAGGATCATCCGGCAACACGGTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAAc            >  2:217488/1‑150 (MQ=255)
                                                                                                                                          aTGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACAt          <  1:499796/150‑1 (MQ=255)
                                                                                                                                           taagACGCGTCAAGCGTCGCATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATc         <  1:600748/147‑1 (MQ=255)
                                                                                                                                           tGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAg                             >  2:483266/1‑130 (MQ=255)
                                                                                                                                           tGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATc         >  2:758025/1‑150 (MQ=255)
                                                                                                                                           tGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATc         <  2:274949/150‑1 (MQ=255)
                                                                                                                                             cGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTg       >  2:1010230/1‑150 (MQ=255)
                                                                                                                                             agACGCGTCAAGCGTCGCATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTg       >  2:898910/2‑150 (MQ=255)
                                                                                                                                                 gCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCtg   <  2:1196898/150‑1 (MQ=255)
                                                                                                                                                 gCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCtg   <  1:1260848/150‑1 (MQ=255)
                                                                                                                                                  cTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCtgt  <  1:758025/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGT  >  NC_000913/2894705‑2895000

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