Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,608,721 A→C 8.5% intergenic (+72/+30) zntA → / ← tusA zinc, cobalt and lead efflux system/mnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,608,7210AC 8.5% 1468.1 / 30.3 520intergenic (+72/+30)zntA/tusAzinc, cobalt and lead efflux system/mnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base A (303/160);  new base C (2/54);  total (305/215)
Fisher's exact test for biased strand distribution p-value = 1.29e-20
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 545 total aligned reads displayed.

TGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGC  >  NC_000913/3608572‑3608870
                                                                                                                                                      |                                                                                                                                                     
  gctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaaa                                                                                                                                                     >  2:5167600/1‑150 (MQ=255)
  gctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaaa                                                                                                                                                     >  1:3381198/1‑150 (MQ=255)
    tggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAgg                                                                                                                                                   >  1:3928248/1‑150 (MQ=255)
    tggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAgg                                                                                                                                                   <  2:4122594/150‑1 (MQ=255)
           gCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCTAATGCGTTAAGTTTGTTGCGCCGGAGATACGGCAAACCTATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGTTGAGGGTTAGGGTGAGGGGAAAAGGCGGCATc                                                                                                                                            <  2:4099293/150‑1 (MQ=38)
             tGGCTGATAGGGGGGCGAGGGTGGTGGTGACAGCTATTGCTTTAATTTTGTTGCGCAGTATATAAGGCAAACCGATCGAAACTTTGAGCGCGATGGGTCCCCTCGCCCCTGTGGGGAGAGGTTTCGGGTGTGGGGAAAAGGCGGCATCGa                                                                                                                                          <  2:1935495/150‑1 (MQ=255)
             tGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGa                                                                                                                                          >  2:7016721/1‑150 (MQ=255)
                 aGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGcc                                                                                                                                      >  2:6692021/1‑150 (MQ=255)
                   aTAGGGGGGCGCCGGTGTTGGTTACTGGGCTTGGGTTATTTTTTTTGCGCTGGTGATAAGGAAAACCTTACGCAAATTTGGGCGGGTCGGGCCCCCCCCCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCaa                                                                                                                                    <  2:4602245/150‑1 (MQ=12)
                   aTACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCaa                                                                                                                                    >  2:6459427/1‑150 (MQ=255)
                       gggggCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAg                                                                                                                                >  1:6841297/1‑150 (MQ=255)
                          ggggagggTGCTGGTGACAGCGAATGCGTTAAGATTGTGGCGCAGGATATAAGGCAACCCGATCGAAACTTTGTGCGCGATCGGCCCCCTCGCCCCTCTGGGGTGTGGGTTTGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGccc                                                                                                                             <  2:2459270/144‑1 (MQ=255)
                                  tGCTGGTGACTGCGAATGCGTTATGATTGTTGCGCAGGAGATAAGGCAACCCTATCGCAACATTGAGGGCGATGGGTCCCCTCGCCCCTGTGGGGTGTGGGTTAGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAAc                                                                                                                     <  2:7599559/150‑1 (MQ=255)
                                    cTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCg                                                                                                                   >  2:2900417/1‑150 (MQ=255)
                                     tGGTGACTGCGAATGCGTTAATTTTGTGGCGCAGGATATAAGCCAAACCGAGCGCAACCTTGCGCGCGATCGGCCCCCTCGCCCCTCTGGGGAGTGGGTTCGGGTGTGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGc                                                                                                                  <  1:8226731/150‑1 (MQ=255)
                                     gggTGCCCGCGATTGCGTTATGTTTGTTGCGCTGGAGATATGGAAAACCGATCGAAACTTTGAGCGCGTTCGGTCCCCGCGCCCCTGTGGGTAGTGGTTTTGGTTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGc                                                                                                                  <  2:347884/149‑1 (MQ=255)
                                      ggTGACTGCGATTGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGCCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGcc                                                                                                                 <  2:6882978/150‑1 (MQ=255)
                                           cAGCTAATGGGTTAAGATTGTTGCGCAGGATATACGGCAAACCGATCGCAACTTTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGTAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTAc                                                                                                            <  2:2726303/150‑1 (MQ=38)
                                            aGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACg                                                                                                           >  2:8375564/1‑150 (MQ=255)
                                             gCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGa                                                                                                          >  1:6370000/1‑150 (MQ=255)
                                              cGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGaa                                                                                                         >  1:5571942/1‑150 (MQ=255)
                                              cGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGaa                                                                                                         >  1:1867385/1‑150 (MQ=255)
                                               gAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGact                                                                                                        >  2:2301965/1‑148 (MQ=255)
                                                     gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                                                  >  2:6072963/1‑150 (MQ=255)
                                                      ttAAGTTTGTTGCGCTGGATTTACGGCAAACCGATCGCAACTTTGAGCGCGATCGGTCCCCTCGCCCCTGTGGGGAGCGGGTTTGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATaa                                                                                                 <  2:2764645/150‑1 (MQ=255)
                                                        taatttgttgcCGCAGGAGAAAAGCCAACCCTATCCAAAATTTGTGCGCGATCGGCCCCCCCCCCCCTCTGGGGTGTGGGTTTGGGTTCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAAc                                                                                                <  1:5843734/146‑1 (MQ=255)
                                                         aGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGa                                                                                              >  1:899166/1‑150 (MQ=255)
                                                           aTTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATa                                                                                            >  1:8115590/1‑150 (MQ=255)
                                                                ggcgcCGGTGATTAGGCAACCCGAGCTCAACTTTGCGGGCGAGGGGTCCCCGCGCCCCTGGGGGGTGTGGGTTTGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCa                                                                                       <  1:3302125/149‑1 (MQ=255)
                                                                   cgcAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGt                                                                                     >  1:3114301/1‑150 (MQ=255)
                                                                     cAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTcc                                                                                   >  2:4976790/1‑150 (MQ=255)
                                                                        gcgaTATGGCAACCCGATCGCAACCTTGCGCGCGTTGGGCCCCCTCGCCCCTCGGGGGAGTGGGTTCGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATc                                                                                <  2:5840006/148‑1 (MQ=255)
                                                                        gagaTAAGGCAAACCGATCGCAACATTGAGCGCGACCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATc                                                                                >  2:5744125/1‑150 (MQ=255)
                                                                          gaTAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCg                                                                              >  1:5173666/1‑150 (MQ=255)
                                                                          gaTAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCg                                                                              >  1:843978/1‑150 (MQ=255)
                                                                             aaGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGtct                                                                           >  1:6223009/1‑150 (MQ=255)
                                                                             aaGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGtct                                                                           >  2:2126181/1‑150 (MQ=255)
                                                                               ggaaaaCCGATCGAAACTTTGAGCGCGATCGGTCCCCTCGCCCCTCGGGGGTGAGGTTTAGGGTGAGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTctct                                                                         <  2:7730506/147‑1 (MQ=255)
                                                                               ggCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTctct                                                                         >  2:923971/1‑150 (MQ=255)
                                                                                gCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                                                                        >  2:6095331/1‑150 (MQ=255)
                                                                                gCAAACCGATCGAAACATTGCGCGCGATGGGCCCCCTCGCCCCTGGGGGGAGTGGGTTTGGGTGTGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                                                                        <  1:4307393/150‑1 (MQ=255)
                                                                                  aaccctctCGAAACATTGCGCGCGATCGGTCCCCTCGCCCCTCTGGGGTGTGGTTTTGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                                                      <  1:4948907/143‑1 (MQ=255)
                                                                                  aaaCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                                                      >  1:286316/1‑150 (MQ=255)
                                                                                  aaaCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                                                      >  2:6506494/1‑150 (MQ=255)
                                                                                  aaaCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCAACCGTCTCTTTa                                                                      >  2:4723402/1‑150 (MQ=255)
                                                                                      cGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTATACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCaa                                                                  >  1:6626473/1‑150 (MQ=255)
                                                                                       gATCGCAACTTTGAGCGCGCTCGGTCCCCTCGCCCCTGTGGGGTGAGGGTTCGGGTGTGGGTACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                                                 <  1:3480447/150‑1 (MQ=255)
                                                                                       gATCGCAACATTGCGCGCGATCGGTCCCCTCGCCCCTGTGGGGTGTGGGTTAGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                                                 <  1:1085835/150‑1 (MQ=255)
                                                                                       gATCGAAACTTTGCGCGCGATGGGTCCCCTCGCCCCTCTGGGTTGTGGTTTTGGTTTGGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                                                 <  2:3158406/150‑1 (MQ=255)
                                                                                         tCGAAACATTGCGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGTGGGTTCGGGTGAGGGCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCa                                                               <  1:1011928/150‑1 (MQ=39)
                                                                                          cGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCACCCAg                                                              >  1:5685907/1‑150 (MQ=255)
                                                                                            cAACATTGCGCGCGATCGGTCCCCTCGCCCCTGTGGGGAGAGGTTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGtt                                                            <  2:3329671/150‑1 (MQ=255)
                                                                                                aTTGCGCGCGAGCGGTCCCCCCGCCCCTGTGGGGTGTGGTTTCGGGTGAGGGCACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtg                                                        <  1:7893880/150‑1 (MQ=255)
                                                                                                 ttGCGCGCTATGGGCCCCCTCGCCCCTCGGGGGCGTGGTTTTGGGTGTGGGTACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtgt                                                       <  1:1528807/150‑1 (MQ=255)
                                                                                                 ttGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtgt                                                       >  2:2681585/1‑150 (MQ=255)
                                                                                                 tggAGCGCGATCGGTCCCCTCGCCCCTGTGGGGAGTGGGTTTGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtgt                                                       <  1:7282820/148‑1 (MQ=255)
                                                                                                  tGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGtt                                                      >  2:8426152/1‑150 (MQ=255)
                                                                                                    cgcgcgcTCGGCCCCCTCCCCCCGGGGGGGTGTGGGTTTGGGTGAGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTcc                                                    <  1:3713995/149‑1 (MQ=255)
                                                                                                    aGCGCGATCGGTCCCCTCGCCCCTCTGGGGTGTGGGTTCGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTcc                                                    <  1:4597491/150‑1 (MQ=255)
                                                                                                     gcgcGATGGGCCCCCTCGCCCCTGTGGGGTTTGGTTTTGGGTGCGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCa                                                   <  1:7208281/150‑1 (MQ=255)
                                                                                                     gcgcGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCa                                                   >  2:420095/1‑150 (MQ=255)
                                                                                                     gcgcGAGCGGCCCCCTCGCCCCTGGGGGGTGTGGGTTCGGGTGTGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCa                                                   <  1:7395087/150‑1 (MQ=255)
                                                                                                      cgcgATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCAt                                                  >  1:7761533/1‑150 (MQ=255)
                                                                                                       gcgATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATa                                                 >  1:6766887/1‑150 (MQ=255)
                                                                                                       gcgATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATa                                                 >  1:1213647/1‑150 (MQ=255)
                                                                                                        cgATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaa                                                >  1:2733789/1‑150 (MQ=255)
                                                                                                        cgATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaa                                                >  1:4199858/1‑150 (MQ=255)
                                                                                                        cgAGCGGTCCCCGCGCCCCTCGGGGGTGTGTGTTTGGGTGCGGGGACCAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCGTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaa                                                <  1:262238/150‑1 (MQ=255)
                                                                                                         gATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaaa                                               >  2:1765589/1‑150 (MQ=255)
                                                                                                         gATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaaa                                               >  2:1018773/1‑150 (MQ=255)
                                                                                                         gATCGGCCCCCTCGCCCCTCGGGGGAGCGGGTTTGGGTGTGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaaa                                               <  1:3263771/150‑1 (MQ=255)
                                                                                                          aTCGGTCCCCTCGCCCCTCTGGGGAGTGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAg                                              <  1:4073229/150‑1 (MQ=255)
                                                                                                              gTCCCCTCCCCCCTCGGGGTTTTGGTTTTGGGTGAGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTAc                                          <  2:8000488/150‑1 (MQ=255)
                                                                                                                ccccGCGCCCCTGGGGGTTGTGGGTTTGGGTGTGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaa                                        <  2:3698750/150‑1 (MQ=255)
                                                                                                                 cccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaaa                                       >  2:1371640/1‑150 (MQ=255)
                                                                                                                 cccTCGCCCCTGTGGGTTGTGGTTTTGGGTGTGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaaa                                       <  2:1634603/150‑1 (MQ=255)
                                                                                                                    tCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAcc                                    >  1:7058427/1‑150 (MQ=255)
                                                                                                                      gCCCTTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa                                  <  2:5595379/150‑1 (MQ=255)
                                                                                                                      gCCCCTGTGGGGTGTGGGTTCGGGTGAGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa                                  <  1:7441648/150‑1 (MQ=255)
                                                                                                                      gCCCCTCTGGGGAGTGGGTTAGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa                                  <  1:7005632/150‑1 (MQ=255)
                                                                                                                         ccTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGa                               >  2:8106027/1‑150 (MQ=255)
                                                                                                                          ctctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGaa                              >  1:3645827/1‑150 (MQ=255)
                                                                                                                              ggggAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATc                          >  1:3765489/1‑150 (MQ=255)
                                                                                                                                ggAGTGGTTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATcgc                        <  1:6078542/150‑1 (MQ=255)
                                                                                                                                ggAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATcgc                        >  1:1428061/1‑150 (MQ=255)
                                                                                                                                    agGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGt                    >  1:4871103/1‑150 (MQ=255)
                                                                                                                                     gggTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGACTATCGCGGGTa                   >  1:6408121/1‑150 (MQ=255)
                                                                                                                                        ttAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTg                >  1:4487190/1‑150 (MQ=255)
                                                                                                                                          tgggTGCGGGCACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGc              <  1:2966621/149‑1 (MQ=255)
                                                                                                                                          aGGGTGAGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGc              <  2:8300814/150‑1 (MQ=255)
                                                                                                                                             gTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCgg           >  2:3964404/1‑150 (MQ=255)
                                                                                                                                              tgtggggAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATATCGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGa          <  1:5110092/147‑1 (MQ=255)
                                                                                                                                               gAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGAt         <  2:6324243/150‑1 (MQ=255)
                                                                                                                                                tggggACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGAtc        <  2:3015680/149‑1 (MQ=255)
                                                                                                                                                cggggAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGAtc        <  1:8388191/149‑1 (MQ=255)
                                                                                                                                                   gtaaaaGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGAtcgtc     <  2:7312993/148‑1 (MQ=255)
                                                                                                                                                   ggAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGAtcgtc     <  2:4322865/150‑1 (MQ=255)
                                                                                                                                                    gAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTATAGCAACCAGTTCGTGTTCCATAAAGGAACAAAACCCAGGAATATCGCGGGTAGTGGCCGGAtcgtcg    >  1:2338558/1‑150 (MQ=255)
                                                                                                                                                     aaaaGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCgg   <  2:5340268/150‑1 (MQ=255)
                                                                                                                                                     aaaaGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCgg   >  2:4444775/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGC  >  NC_000913/3608572‑3608870

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: