Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,501,888 +C coding (618/2496 nt) b2383 ← putative PTS system enzyme IIA component, enzyme I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,501,8841.C98.3% 380.3 / ‑0.9 119coding (622/2496 nt)b2383putative PTS system enzyme IIA component, enzyme I
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base C (66/51);  total (68/51)
Fisher's exact test for biased strand distribution p-value = 5.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

Only 100 of 140 total aligned reads displayed.

AATGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGT  >  NC_000913/2501792‑2501971
                                                                                             |                                                                                       
aaTGTCGCTAACACGTTCGCGCAGATAATCACTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAg                                                                                          >  1:3132950/1‑93 (MQ=255)
 aTGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGc                                                                                        >  2:1129984/1‑93 (MQ=255)
   gTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGccc                                                                                      >  2:1225221/1‑93 (MQ=255)
     cGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGccccc                                                                                    <  1:1019045/93‑2 (MQ=255)
     cGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGccccc                                                                                    <  2:1130542/93‑2 (MQ=255)
     cGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGccccc                                                                                    <  1:3856879/93‑2 (MQ=255)
      gCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGccccct                                                                                   <  1:1656344/93‑3 (MQ=255)
      gCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGccccct                                                                                   <  2:2927912/93‑3 (MQ=255)
      gCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGccccct                                                                                   >  1:1412941/1‑91 (MQ=255)
       cTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGccccctg                                                                                  <  2:1640314/93‑4 (MQ=255)
       cTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGATTGGGGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGcccccta                                                                                  >  2:3604562/1‑90 (MQ=255)
        tAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGccccctga                                                                                 >  1:1260642/1‑89 (MQ=255)
           cacGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGc                                                                               >  1:4117843/1‑93 (MQ=255)
           cacGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGc                                                                               <  1:2752570/93‑1 (MQ=255)
            acGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCt                                                                              <  2:3853588/93‑1 (MQ=255)
             cGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTg                                                                             <  1:413897/93‑1 (MQ=255)
               ttCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGGTCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGtt                                                                           <  2:3345700/93‑1 (MQ=255)
               ttCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGtt                                                                           <  2:3345699/93‑1 (MQ=255)
                tCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTc                                                                          <  2:2369846/93‑1 (MQ=255)
                 cgcgCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCt                                                                         <  2:28000/93‑1 (MQ=255)
                 cgcgCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCt                                                                         <  2:28001/93‑1 (MQ=255)
                      aGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGtcat                                                                    >  1:1617245/1‑93 (MQ=255)
                      aGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGtcat                                                                    >  2:3700297/1‑93 (MQ=255)
                       gATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGtcatc                                                                   >  2:2665480/1‑93 (MQ=255)
                        aTAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGtcatca                                                                  <  2:2887087/93‑1 (MQ=255)
                         tAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAg                                                                 <  2:85910/93‑1 (MQ=255)
                         tAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAg                                                                 <  1:3939343/93‑1 (MQ=255)
                          aaTCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAgg                                                                >  2:1552116/1‑93 (MQ=255)
                           aTCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAggg                                                               >  2:1625088/1‑92 (MQ=255)
                           aTCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGc                                                               >  2:1552144/1‑93 (MQ=255)
                           aTCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGc                                                               >  2:761111/1‑93 (MQ=255)
                           aTCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGc                                                               >  2:589361/1‑93 (MQ=255)
                           aTCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGc                                                               <  1:263184/93‑1 (MQ=255)
                              gCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGAc                                                            <  2:3640180/93‑1 (MQ=255)
                                tGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACgg                                                          >  1:342877/1‑93 (MQ=255)
                                tGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACgg                                                          >  2:165867/1‑93 (MQ=255)
                                   cagaggcagaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGCCATCAGGCGACGGata                                                       <  1:2075570/93‑1 (MQ=255)
                                    agaggcagaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGatat                                                      >  1:1353547/1‑93 (MQ=255)
                                    agaggcagaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGatat                                                      >  1:2132060/1‑93 (MQ=255)
                                    agaggcagaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGatat                                                      <  1:1897084/93‑1 (MQ=255)
                                     gaggcagaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATAtt                                                     >  1:2302279/1‑93 (MQ=255)
                                     gaggcagaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATAtt                                                     >  2:602054/1‑93 (MQ=255)
                                      aggcagaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTg                                                    >  1:235184/1‑93 (MQ=255)
                                      aggcagaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTg                                                    <  1:1138900/93‑1 (MQ=255)
                                        gcagaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGcc                                                  >  2:2601992/1‑93 (MQ=255)
                                          agaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTg                                                >  2:2042497/1‑93 (MQ=255)
                                           gaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGc                                               <  1:3619801/93‑1 (MQ=255)
                                              aGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCaaa                                            >  2:4046460/1‑93 (MQ=255)
                                               gTTTGGCGCAAACTTGCTCCATATTGCTGATGATCGCCTGCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAAt                                           <  2:1219270/93‑1 (MQ=255)
                                                tttGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAAtt                                          <  2:813098/93‑1 (MQ=255)
                                                tttGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAAtt                                          >  1:2970423/1‑93 (MQ=255)
                                                tttGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAAtt                                          >  1:2970424/1‑93 (MQ=255)
                                                tttGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAAtt                                          >  2:2462950/1‑93 (MQ=255)
                                                   ggCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATtcat                                       >  1:4216227/1‑93 (MQ=255)
                                                   gccgcAAACATGCTACATATTGCTGATGATCGCTGCCTGCAGCCCCCTGATGCTGTTCTGTCATCAGGTGACGGATATTGCCTGCAAATtcat                                       <  2:1844300/91‑1 (MQ=255)
                                                    gcgcAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATtcatc                                      >  1:4456159/1‑93 (MQ=255)
                                                     cgcAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATtcatca                                     <  2:4330643/93‑1 (MQ=255)
                                                      gcAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTcatcat                                    <  1:3224805/93‑1 (MQ=255)
                                                      gaaaaCCTGCTCCATTTTGCTGATGATCTCCTGCACCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTcatcat                                    <  1:3841579/91‑1 (MQ=255)
                                                       cAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCatcatc                                   <  2:4396440/93‑1 (MQ=255)
                                                       cAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCatcatc                                   >  1:4251941/1‑93 (MQ=255)
                                                       cAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCatcatc                                   >  1:4251942/1‑93 (MQ=255)
                                                        caaCCTGCTCCATATTGCTGATGATCGCCGTCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATcc                                  <  2:992816/92‑1 (MQ=255)
                                                        aaaCCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATcc                                  >  2:1619263/1‑93 (MQ=255)
                                                         aaCCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCt                                 >  2:778387/1‑93 (MQ=255)
                                                         aaCCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCt                                 <  2:4486933/93‑1 (MQ=255)
                                                           ccTGCTCCATATTGCTGATGATCGCCTCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGa                               <  1:2902252/93‑1 (MQ=255)
                                                           ccTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGa                               >  2:3657437/1‑93 (MQ=255)
                                                               cTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCa                           <  2:2009370/93‑1 (MQ=255)
                                                               cTCCATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCa                           <  2:3055419/93‑1 (MQ=255)
                                                                 ccATATTGCTGCTGATCGCCTCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGc                         <  2:3766572/93‑1 (MQ=255)
                                                                 ccATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGc                         >  1:2032152/1‑93 (MQ=255)
                                                                 ccATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGc                         <  2:310488/93‑1 (MQ=255)
                                                                 ccATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGc                         >  1:2523739/1‑93 (MQ=255)
                                                                  cATATTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCg                        >  1:4487325/1‑93 (MQ=255)
                                                                   atatTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGCATCAGCGa                       >  1:1011014/1‑93 (MQ=255)
                                                                    tatTGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGAc                      >  1:74908/1‑93 (MQ=255)
                                                                       tGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACaaa                   >  1:4709827/1‑93 (MQ=255)
                                                                       tGCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACaaa                   >  1:4709826/1‑93 (MQ=255)
                                                                        gCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAAt                  <  1:3559154/93‑1 (MQ=255)
                                                                        gCTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAAt                  <  1:519898/93‑1 (MQ=255)
                                                                         cTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATg                 <  2:3352315/93‑1 (MQ=255)
                                                                         cTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATg                 <  1:4724051/93‑1 (MQ=255)
                                                                         cTGATGATCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATg                 <  2:124313/93‑1 (MQ=255)
                                                                            atgatCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGgc              <  1:731431/93‑1 (MQ=255)
                                                                             tgatCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGgcg             >  2:2010445/1‑93 (MQ=255)
                                                                              gatCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGgcgc            >  2:1802744/1‑93 (MQ=255)
                                                                              gatCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGgcgc            >  2:1575628/1‑93 (MQ=255)
                                                                              gatCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGgcgc            <  1:2729020/93‑1 (MQ=255)
                                                                              gatCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGgcgc            >  2:1391540/1‑93 (MQ=255)
                                                                                tCGCCGCCCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTg          <  1:2449123/93‑1 (MQ=255)
                                                                                  gccgccCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAg        >  1:3810561/1‑93 (MQ=255)
                                                                                  gccgccCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAg        >  2:2359914/1‑93 (MQ=255)
                                                                                   ccgccCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAgg       >  1:4731813/1‑93 (MQ=255)
                                                                                   ccgccCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAgg       <  1:692996/93‑1 (MQ=255)
                                                                                     gccCCCAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGAt     >  1:3573794/1‑93 (MQ=255)
                                                                                      cccccAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATa    <  2:2986367/93‑1 (MQ=255)
                                                                                      cccccAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATa    <  1:1954808/93‑1 (MQ=255)
                                                                                      cccccAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATa    >  2:1883817/1‑93 (MQ=255)
                                                                                        cccAGCCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGt  >  2:1411089/1‑93 (MQ=255)
                                                                                             |                                                                                       
AATGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGT  >  NC_000913/2501792‑2501971

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