Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 935,131 T→A 35.0% D636E (GAT→GAA ftsK → DNA translocase at septal ring sorting daughter chromsomes

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913935,1310TA35.0% 127.9 / 126.2 120D636E (GAT→GAAftsKDNA translocase at septal ring sorting daughter chromsomes
Reads supporting (aligned to +/- strand):  ref base T (36/42);  new base A (17/25);  total (53/67)
Fisher's exact test for biased strand distribution p-value = 5.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.67e-01

Only 100 of 120 total aligned reads displayed.

TCGTGAACTGGCGTCTTACGGTATTAAGCTGCCCTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAGTTCGCCCAGACACAGCAGCAACGCTATGGCGAACAGTATCA  >  NC_000913/935035‑935226
                                                                                                |                                                                                               
tcgtGAACTGGCGTCTTACGGTATTAAGCTGCCCTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAg                                                                                              <  2:292857/100‑1 (MQ=255)
tcgtGAACTGGCGTCTTACGGTATTAAGCTGCCCTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAg                                                                                              <  1:273075/100‑1 (MQ=255)
 cgtGAACTGGCGTCTTACGGTATTAAGCTGCCCTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGt                                                                                             <  1:2735869/100‑1 (MQ=255)
       cTGGCGTCTTACGGTATTAAGCTGCCCTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACg                                                                                       >  1:920329/1‑100 (MQ=255)
       cTGGCGTCTTACGGTATTAAGCTGCCCTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACg                                                                                       <  2:1935403/100‑1 (MQ=255)
       cTGGCGTCTTACGGTATTAAGCTGCCCTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGAACAGTACAACg                                                                                       >  1:271350/1‑100 (MQ=255)
        tGGCGTCTTACGGTATTAAGCTGCCCTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACga                                                                                      <  1:2529796/100‑1 (MQ=255)
          gCGTCTTACGGTATTAAGCTGCCCTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGAACAGTACAACgatg                                                                                    <  2:2623668/100‑1 (MQ=255)
            gTCTTACGGTATTAAGCTGCCCTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACgatgat                                                                                  <  2:33498/100‑1 (MQ=255)
            gTCTTACGGTATTAAGCTGCCCTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACgatgat                                                                                  >  2:1222733/1‑100 (MQ=255)
             tCTTACGGTATTAAGCTGCCCTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGatgatg                                                                                 <  1:3523215/100‑1 (MQ=255)
              cTTACGGTATTAAGCTGCCCTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTATGGCGATCAGTACAACGAtgatga                                                                                >  1:2819660/1‑100 (MQ=255)
               ttACGGTATTAAGCTGCCCTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGaa                                                                               >  2:538584/1‑100 (MQ=255)
                tACGGTATTAAGCTGCCCTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGAACAGTACAACGATGATGaaa                                                                              >  1:672808/1‑100 (MQ=255)
                 aCGGTATTAAGCTGCCCTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGAAAt                                                                             <  1:860884/100‑1 (MQ=255)
                  cGGTATTAAGCTGGCCTCACAGCGTGCGGCGGAAGAAAAAGCGCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGAACAGTACAACGATGATGAAATc                                                                            >  2:2504460/1‑100 (MQ=255)
                  cGGTATTAAGCTGCCCTCACAGCGTGCGGCGGAAGAACAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGAAATc                                                                            >  1:916960/1‑100 (MQ=255)
                  cGGTATTAAGCTGCCCTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGAAATc                                                                            >  2:1762808/1‑100 (MQ=255)
                  cGGTATTAAGCTGCCCTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGAAATc                                                                            <  1:1194467/100‑1 (MQ=255)
                    gTATTAAGCTGCCCTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGAAATcga                                                                          >  2:2584998/1‑100 (MQ=255)
                    gTATTAAGCTGCCCTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGAACAGTACAACGATGATGAAATcga                                                                          >  2:1287648/1‑100 (MQ=255)
                      aTTAAGCTGCCCTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGAAATcgatg                                                                        <  2:21708/100‑1 (MQ=255)
                      aTTAAGCTGCCCTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGAACAGTACAACGATGATGAAATcgatg                                                                        <  1:1688883/100‑1 (MQ=255)
                      aTTAAGCCGCCATCACATCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAACTCAATCAGTACTATTCTTGCGAACAGTACAACGATGATGAAATcgatg                                                                        <  1:2493618/100‑1 (MQ=255)
                       ttAAGCTGCCCTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGAAATcgatgc                                                                       >  1:1994370/1‑100 (MQ=255)
                       ttAAGCTGCCCTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGAACAGTACAACGATGATGAAATcgatgc                                                                       <  2:3229397/100‑1 (MQ=255)
                        tAAGCTGCCCTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGAAATcgatgcg                                                                      <  1:408962/100‑1 (MQ=255)
                        tAAGCTGCCCTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGAAATcgatgc                                                                       >  2:2806116/1‑99 (MQ=255)
                          aGCTGCCCTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGTCCAGCGCAATCAGTACGATTCTGGCGAACAGTACAACGATGATGAAATcgatgcgat                                                                    <  1:2598534/100‑1 (MQ=255)
                          aGCTGCCCTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGAAATcgatgcgat                                                                    <  2:1714299/100‑1 (MQ=255)
                           gCTGCCCTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGAAATcgatgcgatg                                                                   >  2:2060983/1‑100 (MQ=255)
                           gCTGCCCTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGAACAGTACAACGATGATGAAATcgatgcgatg                                                                   <  2:3186797/100‑1 (MQ=255)
                            cTGCCCTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGAACAGTACAACGATGATGAAATCgatgcgatgc                                                                  <  2:329255/100‑1 (MQ=255)
                            cTGCCCTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAAGCAGTACGATTCTGGCGATCAGTACAACGATGATGAAATCtatgcgatgc                                                                  <  1:1621585/100‑1 (MQ=255)
                            cTGCCCTCACAGCGTGCGGCGGAAGAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGAACAGTACAACGATGATGAAATCGATACGATgca                                                                 >  1:2527481/1‑100 (MQ=255)
                             tGCCCTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGAAATCGATGCGATgca                                                                 <  1:3077395/100‑1 (MQ=255)
                              gCCCTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGAAATCGATGCGATgcag                                                                >  1:3601731/1‑100 (MQ=255)
                              gCCCTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGAACAGTACAGCGATGATGAAATCGCTGCGATgcag                                                                >  1:910217/1‑100 (MQ=255)
                               cccTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCATTACAACGATGATGAAATTGATGCGATgcatc                                                               <  2:2492635/100‑3 (MQ=255)
                               cccTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGAAATCGATGCGATgcagc                                                               >  2:2201921/1‑100 (MQ=255)
                                ccTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGAACAGTACAACGATGATGAAATCGATGCGATgcagca                                                              <  2:3396208/100‑1 (MQ=255)
                                 cTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGAAATCGATGCGATGcagcag                                                             >  1:544373/1‑100 (MQ=255)
                                 cTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGAAATCGATGCGATGcagcag                                                             >  1:2786563/1‑100 (MQ=255)
                                 cTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGAAATCGATGCGATGcagcag                                                             <  2:2589230/100‑1 (MQ=255)
                                 cTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGAAATCGATGCGATGcagcag                                                             <  1:1500965/100‑1 (MQ=255)
                                 cTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGAACAGTACAACGATGATGCAATCGATGCGATgcagca                                                              >  1:3279253/1‑99 (MQ=255)
                                     caGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATg                                                         >  1:2219466/1‑100 (MQ=255)
                                         gggcggcgGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGAACAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACt                                                     <  1:1324686/98‑1 (MQ=255)
                                         gTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACt                                                     <  2:2680303/100‑1 (MQ=255)
                                         gTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGAACAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACt                                                     <  2:3556122/100‑1 (MQ=255)
                                         gTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGAACAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACt                                                     <  1:3160461/100‑1 (MQ=255)
                                             ggcggAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGAACAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCa                                                 >  1:32674/1‑100 (MQ=255)
                                             ggcggAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGAACAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCa                                                 <  2:525908/100‑1 (MQ=255)
                                                ggAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGAACAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGt                                              >  2:2054460/1‑100 (MQ=255)
                                                 gaataaAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGAACAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTgg                                                   <  1:3268001/94‑1 (MQ=255)
                                                 gaagaaAAAGCCCGTGAAGCCCAGCGCAATCAGTACTATTCTGGCGAACAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTc                                             >  2:2615703/1‑100 (MQ=255)
                                                 gaagaaAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGAACAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACg                                               <  1:1543378/98‑1 (MQ=255)
                                                  aagaaAAAGCCCGTAAAGCCCAGCGCAATCAGTACGATTCTGGCGAACAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTc                                             <  2:715628/99‑1 (MQ=255)
                                                   agaaAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAg                                           <  1:1349151/100‑1 (MQ=255)
                                                   agaaAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAg                                           >  2:835541/1‑100 (MQ=255)
                                                    gaaAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGAACAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCaagcgtga                                          >  2:1260538/1‑93 (MQ=255)
                                                     aaaaaGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAGtt                                         <  1:1390819/100‑1 (MQ=255)
                                                     aaaaaGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAGtt                                         >  2:2357808/1‑100 (MQ=255)
                                                     aaaaaGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAGtt                                         >  2:1223569/1‑100 (MQ=255)
                                                     aaaaaGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGAACAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAGtt                                         >  1:2525556/1‑100 (MQ=255)
                                                      aaaaGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGAACAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAGTTc                                        >  2:2693170/1‑100 (MQ=255)
                                                        aaGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAGTTCGc                                      <  1:3649012/100‑1 (MQ=255)
                                                        aaGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGAACAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAGTTCGc                                      <  2:1017937/100‑1 (MQ=255)
                                                         tctccGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGAACAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAGTTCGcc                                     <  2:1874451/97‑1 (MQ=255)
                                                         aGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAGTTCGcc                                     <  1:2728001/100‑1 (MQ=255)
                                                         aGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAGTTCGcc                                     <  2:657364/100‑1 (MQ=255)
                                                           cccGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAGTTCGCCCa                                   >  1:2405499/1‑100 (MQ=255)
                                                           cccGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAGTTCGCCCa                                   <  1:2201908/100‑1 (MQ=255)
                                                            ccGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGAACAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAGTTCGccaga                                  >  1:1874439/1‑97 (MQ=255)
                                                              gTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAGTTCGCCCAGac                                <  2:152415/100‑1 (MQ=255)
                                                              gTGAAGCCCAGCGCAATCAGTACGATTCTGGCGAACAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAGTTCGCCCAGac                                <  2:3124203/100‑1 (MQ=255)
                                                              gTGAAGCCCAGCGCAATCAGTACGATTCTGGCGAACAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAGTTCGCCCAGac                                <  2:2735041/100‑1 (MQ=255)
                                                               tGAAACCCAGCGCAATCAGTACGCTTCTGGCGATCAGTACAGCGATGATGAAATCGATGCGATGCAGCAGGATGCACTGGCACGTCAGTTAGCCCAGaca                               >  1:381641/1‑100 (MQ=255)
                                                                gAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAGTTCGCCCAGacac                              <  1:3005964/100‑1 (MQ=255)
                                                                 aaGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAGTTCGCCCAGAcaca                             <  1:2286325/100‑1 (MQ=255)
                                                                  aGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAGTTCGCCCAGACAcag                            <  2:1610846/100‑1 (MQ=255)
                                                                    cccAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAGTTCGCCCAGACAcagca                          >  1:3115383/1‑100 (MQ=255)
                                                                     ccAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAGTTCGCCCAGACAcagcag                         >  2:1195120/1‑100 (MQ=255)
                                                                      cAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAGTTCGCCCAGACACagcagc                        >  1:2070406/1‑100 (MQ=255)
                                                                       aGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAGTTCGCCCAGACACAgcagca                       >  2:2821683/1‑100 (MQ=255)
                                                                           cAATCAGTACGATTCTGGCGAACAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAGTTCGCCCAGACACAGCAGCAACGc                   <  2:562713/100‑1 (MQ=255)
                                                                               cAGTACGATTCTGGCGATCAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAGTTCGCCCAGACACAGCAGCAACGCTATg               <  2:62454/100‑1 (MQ=255)
                                                                               cAGTACGATTCTGGCGAACAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAGTTCGCCCAGACACAGCAGCAACGCTATg               >  2:2664868/1‑100 (MQ=255)
                                                                               cAGTACGATTCTGGCGAACAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAGTTCGCCCAGACACAGCAGCAACGCTATg               <  2:1940141/100‑1 (MQ=255)
                                                                                 gTACGATTCTGGCGAACAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAGTTCGCCCAGACACAGCAGCAACGCTATGGc             <  1:1323903/100‑1 (MQ=255)
                                                                                  tACGATTCTGGCGAACAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAGTTCGCCCAGACACAGCAGCAACGCTATGGCg            <  2:672813/100‑1 (MQ=255)
                                                                                     gATTCTGGCGATCAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAGTTCGCCCAGACACAGCAGCAACGCTATGGCGAAc         <  2:544376/100‑1 (MQ=255)
                                                                                     gATTCTGGCGATCAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAGTTCGCCCAGACACAGCAGCAACGCTATGGCGAAc         <  2:2356350/100‑1 (MQ=255)
                                                                                         cTGGCGATCAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAGTTCGCCCAGACACAGCAGCAACGCTATGGCGAACAGTa     <  1:3498689/100‑1 (MQ=255)
                                                                                         cTGGCGATCAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAGTTCGCCCAGACACAGCAGCAACGCTATGGCGAACAGTa     >  1:539524/1‑100 (MQ=255)
                                                                                         cTGGCGATCAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAGTTCGCCCAGACACAGCAGCAACGCTATGGCGAACAGTa     <  2:695492/100‑1 (MQ=255)
                                                                                         cTGGCGATCAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAGTTCGCCCAGACACAGCAGCAACGCTATGGCGAACAGTa     >  2:3498716/1‑100 (MQ=255)
                                                                                           ggCGATCAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAGTTCGCCCAGACACAGCAGCAACGCTATGGCGAACAGTATc   >  2:396431/1‑100 (MQ=255)
                                                                                            gCGATCAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAGTTCGCCCAGACACAGCAGCAACGCTATGGCGAACAGTATCa  >  1:3601279/1‑100 (MQ=255)
                                                                                            gCGATCAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAGTTCGCCCAGACACAGCAGCAACGCTATGGCGAACAGTATCa  <  1:327249/100‑1 (MQ=255)
                                                                                                |                                                                                               
TCGTGAACTGGCGTCTTACGGTATTAAGCTGCCCTCACAGCGTGCGGCGGAAGAAAAAGCCCGTGAAGCCCAGCGCAATCAGTACGATTCTGGCGATCAGTACAACGATGATGAAATCGATGCGATGCAGCAGGATGAACTGGCACGTCAGTTCGCCCAGACACAGCAGCAACGCTATGGCGAACAGTATCA  >  NC_000913/935035‑935226

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