Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,270,143 C→T 6.5% intergenic (+541/+73) yhaC → / ← rnpB pentapetide repeats‑related protein/RNase P, M1 RNA component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,270,1430CT 6.5% 783.4 / 25.4 276intergenic (+541/+73)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base C (145/108);  new base T (5/13);  total (152/124)
Fisher's exact test for biased strand distribution p-value = 2.47e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 287 total aligned reads displayed.

ACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAG  >  NC_000913/3269997‑3270292
                                                                                                                                                  |                                                                                                                                                     
aCTGACCGATAAGCCTTTTTCTTTTGGTTATAGTGTCGTGGACAGTCATTATTCTTTCTTCACCTCAAAAATCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgtc                                                                                                                                                    <  1:5280545/150‑1 (MQ=11)
aCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgtc                                                                                                                                                    <  1:6281842/150‑1 (MQ=37)
  tGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgtcgt                                                                                                                                                  <  2:2440172/150‑1 (MQ=37)
   gACCGATCAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCAGTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTg                                                                                                                                                 <  1:5785321/150‑1 (MQ=255)
            gCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtcat                                                                                                                                        <  2:2956360/150‑1 (MQ=37)
             ccGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtcatt                                                                                                                                       <  1:484943/150‑1 (MQ=37)
                   tctttTGGGTAGAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGAGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCtt                                                                                                                                 <  1:89267/150‑1 (MQ=255)
                        tGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCACAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGc                                                                                                                            <  1:5551041/150‑1 (MQ=37)
                        tGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGc                                                                                                                            <  1:992224/150‑1 (MQ=37)
                          ggTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGccc                                                                                                                          <  1:3827933/150‑1 (MQ=37)
                            tataGTGTCGTGGACAGTCATTCAACTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACACGTGAAACTGACCGATAAGCCGCTTTATTTTGGGGATAGCGTCGTGGAAAGGAAATCAACATTCAGCCCCt                                                                                                                        >  1:4655014/1‑150 (MQ=11)
                                  gtCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                                                                                                  <  2:6275226/150‑1 (MQ=37)
                                    cGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTGGGGGATAGCGTCGTGGCCAGGCCAACATCATAAAGCCCCACCAAAAGc                                                                                                                >  2:4554087/1‑150 (MQ=11)
                                      tGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATATTTCTGCCCCTCCAAAAGCaa                                                                                                              >  1:2376548/1‑150 (MQ=37)
                                           aGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAcc                                                                                                         <  2:3402953/150‑1 (MQ=37)
                                             tcattcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCCTTCATCTTTCTGCCCCTCCAAAAGCAAAAACccg                                                                                                       >  2:1073237/1‑150 (MQ=37)
                                               attcatCTTTCTGCCCCTCCAAAAGCAAAATCCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACTGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccgcc                                                                                                     >  2:5487383/1‑150 (MQ=255)
                                                 tcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGa                                                                                                   >  1:4574245/1‑150 (MQ=37)
                                                 tcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATATTTCTGCCCCGCCAAAAGCAAAAACCCGCCGa                                                                                                   >  2:5149257/1‑150 (MQ=255)
                                                   atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCGCCAAAAGCAAAAACCCGCCGAAg                                                                                                 >  1:1141241/1‑150 (MQ=37)
                                                   atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTGTTGGGTATAGCGTCGTGGACAGTCATTCAACTTTCTGCCCCGCCAAAAGCAAAAACCCGCCGAAg                                                                                                 >  2:3734233/1‑150 (MQ=255)
                                                    tCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGc                                                                                                >  1:3418063/1‑150 (MQ=37)
                                                     cTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCg                                                                                               <  2:1205086/150‑1 (MQ=17)
                                                      tttCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCgg                                                                                              <  1:239249/150‑1 (MQ=17)
                                                        tCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCGTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGt                                                                                            >  2:6528769/1‑150 (MQ=37)
                                                         cggCCCCTCCACCAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGtt                                                                                           <  2:4799370/148‑1 (MQ=255)
                                                           gCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGGCAGTCATTCATCTTTCTGCCCCGCCCAAAGCAAAAACCCGCCGCAGCggggtgt                                                                                         >  2:5589541/1‑146 (MQ=11)
                                                            ccccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTGTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttttt                                                                                        >  2:3386395/1‑150 (MQ=37)
                                                              ccTCCACAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCCGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAc                                                                                      <  2:3466622/150‑1 (MQ=255)
                                                              ccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTGTTAc                                                                                      >  1:728178/1‑150 (MQ=37)
                                                                 ccAAAAGCAAAAACCCGCCGAAGCGGTTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTa                                                                                   <  2:5797771/150‑1 (MQ=17)
                                                                  cAAAAGCAAAAACCCGCCGAAGTGGGTTTTTACGTAAACCAGGTGAAACTGCCCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaa                                                                                  <  1:1364881/150‑1 (MQ=255)
                                                                  cAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaa                                                                                  <  2:438432/150‑1 (MQ=37)
                                                                  cAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaa                                                                                  >  2:4176096/1‑150 (MQ=37)
                                                                  cAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaa                                                                                  >  1:4042539/1‑150 (MQ=37)
                                                                   aaaaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGACGCGGGTTTTTACGTaaa                                                                                 >  1:319627/1‑150 (MQ=37)
                                                                     aaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATc                                                                               >  1:4459944/1‑150 (MQ=38)
                                                                     aaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATc                                                                               >  1:5615299/1‑150 (MQ=38)
                                                                     aaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATc                                                                               >  1:5994875/1‑150 (MQ=38)
                                                                      aGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCa                                                                              <  1:93387/150‑1 (MQ=39)
                                                                       gCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGGATAGCGTCGTGGACAGGCATTCATCGTTCTGCCCCGCCCACAGCAAAAACCCGCCGATGCGGGTTTTTTCGGAAATCAg                                                                             >  2:4040225/1‑150 (MQ=11)
                                                                        taaaaaCCCGCTTCGGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGGCAGGCATGCAACATTAAGCCCCTCCAAAAACCAAAAACCGCCGAAGCGGGGTTTTTCGTAAAtgatg                                                                            >  1:3044971/2‑146 (MQ=11)
                                                                        cacacaCCCGCCGAAGCGGGTTATTACGTACACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGATTTTACGTAAATCAgg                                                                            <  1:4227799/145‑1 (MQ=255)
                                                                        cAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAgg                                                                            >  1:145303/1‑150 (MQ=38)
                                                                        cAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTGTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAgg                                                                            >  1:1022540/1‑150 (MQ=38)
                                                                         aaaaaCCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGt                                                                           <  2:4667034/150‑1 (MQ=37)
                                                                         aaaaaCCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGGAAATCAGGt                                                                           >  1:3177876/1‑150 (MQ=37)
                                                                         aaaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGt                                                                           >  2:1997787/1‑150 (MQ=17)
                                                                           aaaCCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGa                                                                         >  2:6548407/1‑150 (MQ=37)
                                                                           aaaCCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCAGCCCCTCCAAAATCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGa                                                                         >  1:4380827/1‑150 (MQ=255)
                                                                            aaCCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGaa                                                                        >  2:5716435/1‑150 (MQ=37)
                                                                             aCCCGCCGTAGCGGGTTATTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGaaa                                                                       <  2:1842720/150‑1 (MQ=12)
                                                                               ccgccgAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAATGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACt                                                                     >  1:535421/1‑150 (MQ=37)
                                                                               ccgccgAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACt                                                                     >  2:756772/1‑150 (MQ=37)
                                                                                 gccgAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGa                                                                   >  1:3268048/1‑150 (MQ=37)
                                                                                   cgAAGCGGGTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGAcc                                                                 <  1:859019/150‑1 (MQ=37)
                                                                                      atcGGTTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGAt                                                              <  1:4747941/148‑1 (MQ=12)
                                                                                        cGGGTTTTTACGTACAGCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATaa                                                            <  2:4518257/150‑1 (MQ=21)
                                                                                        cGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTAGGGTATAGCGTCGTGGACAGTCATACATCTTTATTCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATAAGGTGAAACTGACCGATaa                                                            >  2:4253200/1‑150 (MQ=11)
                                                                                         gggTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGCAGCGGGTGTTTACGTAAATCAGGGGAAACTGACCGATAAg                                                           >  1:3001832/1‑150 (MQ=255)
                                                                                         gggTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAg                                                           >  2:490523/1‑150 (MQ=37)
                                                                                          ggTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCTTTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGc                                                          <  1:4882228/150‑1 (MQ=17)
                                                                                           gTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGGCATTCATCTTTCTGCCCCGCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATCAGcc                                                         >  2:4106215/1‑150 (MQ=255)
                                                                                             ttttACGTAATCCAGGTGAAACTGACCGATCAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCgg                                                       <  2:2344799/150‑1 (MQ=12)
                                                                                             tattACGTACACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGGGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCgg                                                       <  2:1460044/148‑1 (MQ=255)
                                                                                                tACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGtt                                                    <  1:5195698/150‑1 (MQ=40)
                                                                                                tACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCAAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGtt                                                    <  1:658824/150‑1 (MQ=255)
                                                                                                 aCGTAAACCCGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTc                                                   <  2:5061760/150‑1 (MQ=255)
                                                                                                       aCCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATGTTTCTGCCCCTGCAAAAGCAAAATACCGCCGAAGCGGGTTTTTTCGTAAATGAGGTGAAACTGTCCGATAAGCCGGGTTATGTCGt                                             >  1:5833681/1‑150 (MQ=14)
                                                                                                           ggTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGGGTTCTTTGGTGGaa                                         >  2:4862211/1‑149 (MQ=255)
                                                                                                            gggAAACTGACCGAGAAGCCGCTTTCTTTTGGGTATCGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGGGGACa                                        <  2:3234177/148‑1 (MQ=255)
                                                                                                            gTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTTAAACAGACCGATAAGCCGGGTTCTGTCGTGGACa                                        >  2:5671049/1‑150 (MQ=255)
                                                                                                            gTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATAAGGTGAAACTGACCGATAAGCCGGGGTCTGTCGTGGACa                                        >  2:2112645/1‑150 (MQ=255)
                                                                                                             tGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTATGCCCCACCAAAAGCAAAAAACCGCCGCAGCGGGTTTTTACGTAAATCAGGTGTAACAGACCGATAAGCCGGGGTCTGTCGTGGACAg                                       >  1:2629095/1‑150 (MQ=14)
                                                                                                               aaaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGtc                                     <  2:5205762/150‑1 (MQ=255)
                                                                                                                  cTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGtcatt                                  <  1:783389/150‑1 (MQ=255)
                                                                                                                   tGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGtcattc                                 <  1:2207498/150‑1 (MQ=255)
                                                                                                                     aCCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGtcattcat                               <  1:1997787/150‑1 (MQ=255)
                                                                                                                     aCCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTAATCTTTATGCCCCACCCAAAGCAAAAAACCCCCCAACCGGGTTTTTTAGTAAAAAAGGTGTAAAAGACCGAAAAGCCCGGGTATGTCGTGGACAGtcattcaa                               >  1:1668119/1‑149 (MQ=9)
                                                                                                                     aCCGATAAGCCGCTTTCGTTTGGGTATAGCGTCGTGGACAGTCAATAATCTTTCTGCCCCACCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGtcattcat                               >  1:504467/1‑150 (MQ=255)
                                                                                                                      ccGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATc                              <  1:4299524/150‑1 (MQ=255)
                                                                                                                        gATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTa                            <  2:1168004/150‑1 (MQ=255)
                                                                                                                        gATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTa                            <  1:3952288/150‑1 (MQ=255)
                                                                                                                        gATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTa                            <  1:354859/150‑1 (MQ=255)
                                                                                                                         aTAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAg                           <  1:4975311/150‑1 (MQ=255)
                                                                                                                         aTAAGCCGCTTTCTTTTGGGTATAGAGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAg                           <  1:3638814/150‑1 (MQ=255)
                                                                                                                          tAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAgg                          <  1:2301465/150‑1 (MQ=255)
                                                                                                                           aaGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGc                         <  2:816826/150‑1 (MQ=255)
                                                                                                                                 ctttctttTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCaa                   >  2:2705192/1‑150 (MQ=255)
                                                                                                                                 ctctctctTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCaa                   <  2:3586946/147‑1 (MQ=255)
                                                                                                                                  tttctttTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAAt                  >  2:6197440/1‑150 (MQ=255)
                                                                                                                                  tttctttTGGGTATAGCGTCGTGGACACTCAATCATCTTTCTGCCCCGCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTAGTGGACAGTCATTCATCTAGGCCAGCAAt                  >  2:3923514/1‑150 (MQ=255)
                                                                                                                                   ttctttTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATc                 <  1:5716435/150‑1 (MQ=255)
                                                                                                                                    tctttTGGGTATAGTGTCGTCGACAGTCAATCATCATTCTGCACCTACAAAACCAAAAACCAGCCTAAGCGTGTTTTTACTTAAATCAGGTTAAACTGACCGATAATCCGGTTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCg                <  2:221147/150‑1 (MQ=18)
                                                                                                                                      ttttGGGTATAGTGTCGTGGACAGTCATTAATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTTACCGATAAGCCGGGTTCTGTAGTGTACAGTGATTCATCTAGGCCATCAATAGCt              >  2:3063717/1‑150 (MQ=17)
                                                                                                                                      ttttGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccgcc                                                                                                     <  1:233268/63‑1 (MQ=35)
                                                                                                                                        ttGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCa            <  1:410080/150‑1 (MQ=255)
                                                                                                                                             tataGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAACAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGc       <  1:2794892/150‑1 (MQ=255)
                                                                                                                                              ataGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCt      <  2:2053888/150‑1 (MQ=255)
                                                                                                                                               taGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTc     >  2:1739924/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
ACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAG  >  NC_000913/3269997‑3270292

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