Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,532,993 A→G D303G (GAC→GGC)  ptsI → PEP‑protein phosphotransferase system enzyme I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,532,9930AG99.7% 1023.6 / ‑5.4 321D303G (GAC→GGC) ptsIPEP‑protein phosphotransferase system enzyme I
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (149/171);  total (149/172)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

Only 100 of 329 total aligned reads displayed.

CTAACATTGGTACGGTTCGTGACGTTGAAGGTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGACCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCAGTGGCTGAAGCGTGTGGCTCGCAAGCGGTTATCGTTCGTA  >  NC_000913/2532903‑2533082
                                                                                          |                                                                                         
cTAACATTGGTACGGTTCGTGACGTTGAAGGTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGcc                                                                                         <  1:5056752/93‑1 (MQ=255)
cTAACATTGGTACGGTTCGTGACGTTGAAGGTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGcc                                                                                         >  1:3478492/1‑93 (MQ=255)
 tAACATTGGTACGGTTCGTGACGTTGAAGGTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCcg                                                                                        >  1:4706393/1‑93 (MQ=255)
 tAACATTGGTACGGTTCGTGACGTTGAAGGTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCcg                                                                                        >  1:7219739/1‑93 (MQ=255)
 tAACATTGGTACGGTTCGTGACGTTGAAGGTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCcg                                                                                        <  1:2208430/93‑1 (MQ=255)
 tAACATTGGTACAGTTCGTGACGTTGAAGGTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCcg                                                                                        <  1:6028613/93‑1 (MQ=255)
  aaCATTGGTACGGTTCGTGACGTTGAAGGTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCcgc                                                                                       <  2:1564305/93‑1 (MQ=255)
  aaCATTGGTACGGTTCGTGACGTTGAAGGTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCcgc                                                                                       >  1:7348869/1‑93 (MQ=255)
    cATTGGTACGGTTCGTGACGTTGAAGGTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGa                                                                                     <  2:4671256/93‑1 (MQ=255)
     aTTGGTACGGTTCGTGACGTTGAAGGTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGAc                                                                                    <  1:3843558/93‑1 (MQ=255)
      ttGGTACGGTTCGTGACGTTGAAGGTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACg                                                                                   >  2:4192547/1‑93 (MQ=255)
        ggTACGGTTCGTGACGTTGAAGGTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCa                                                                                 <  2:1847806/93‑1 (MQ=255)
        ggTACGGTTCGTGACGTTGAAGGTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCa                                                                                 >  2:8062833/1‑93 (MQ=255)
        ggTACGGTTCGTGACGTTGAAGGTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCa                                                                                 <  2:6955198/93‑1 (MQ=255)
         gTACGGTTCGTGACGTTGAAGGTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCAc                                                                                <  2:8007884/93‑1 (MQ=255)
          tACGGTTCGTGACGTTGAAGGTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACt                                                                               <  1:8925543/93‑1 (MQ=255)
          tACGGTTCGTGACGTTGAAGGTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACt                                                                               >  1:271677/1‑93 (MQ=255)
           aCGGTTCGTGACGTTGAAGGTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTg                                                                              >  1:6149172/1‑93 (MQ=255)
            cGGTTCGTGACGTTGGAGGTGCAGAGCGTAACGACGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGc                                                                             <  1:8982993/93‑1 (MQ=255)
            cGGTTCGTGACGTTGAAGGTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGc                                                                             >  1:6290614/1‑93 (MQ=255)
            cGGTTCGTGACGTTGAAGGTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGc                                                                             >  1:4173936/1‑93 (MQ=255)
             ggTTCGTGACGTTGAAGGTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGcc                                                                            <  2:1405552/93‑1 (MQ=255)
              gTTCGTGACGTTGAAGGTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGccc                                                                           >  1:4319893/1‑93 (MQ=255)
              gTTCGTGACGTTGAAGGTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGccc                                                                           <  2:3380964/93‑1 (MQ=255)
               ttCGTGACGTTGAAGGTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCa                                                                          <  1:8757817/93‑1 (MQ=255)
               ttCGTGACGTTGAAGGTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCa                                                                          >  2:815546/1‑93 (MQ=255)
               ttCGTGACGTTGAAGGTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCa                                                                          <  2:3017294/93‑1 (MQ=255)
               ttCGTGACGTTGAAGGTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCa                                                                          <  2:6150687/93‑1 (MQ=255)
                  gTGACGTTGAAGGTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTg                                                                       >  2:8849948/1‑93 (MQ=255)
                  ctGACGTTGAAGGTGCAGAGCGTAACGGCTCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTg                                                                       <  1:5330413/92‑1 (MQ=255)
                   tGACGTTGAAGGTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTga                                                                      >  2:2410629/1‑93 (MQ=255)
                   tGACGTTGAAGGTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTga                                                                      >  2:9372741/1‑93 (MQ=255)
                   tGACGTTGAAGGTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCCGTTCATGGGCCGCGACGCACTGCCCACTga                                                                      >  1:8215121/1‑93 (MQ=255)
                     aCGTTGAAGGTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTgcag                                                                    >  1:1024420/1‑93 (MQ=255)
                     aCGTTGAAGGTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTgaag                                                                    <  2:1395382/93‑1 (MQ=255)
                        ttGAAGGTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTgaagaag                                                                 >  2:1867059/1‑93 (MQ=255)
                         tGAAGGTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTgaagaaga                                                                >  2:8638354/1‑93 (MQ=255)
                          gAAGGTGCAGAGCGTAACGGCTCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTGCATGGGCCGCGACGCACTGCCCACTgaagaagaa                                                               <  1:3527282/93‑1 (MQ=255)
                           aaGGTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAAc                                                              <  2:2838208/93‑1 (MQ=255)
                             ggTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAg                                                            >  1:7818826/1‑93 (MQ=255)
                              gTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGt                                                           <  2:7522506/93‑1 (MQ=255)
                              gTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGt                                                           >  1:8152384/1‑93 (MQ=255)
                               tGCAGAGCGTAACGGCGCTTAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGtt                                                          <  2:2603986/93‑1 (MQ=255)
                               tGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGtt                                                          <  1:9971445/93‑1 (MQ=255)
                                gCAGAGCGTAACGGCTCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGttt                                                         <  2:3916649/93‑1 (MQ=255)
                                gCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGttt                                                         >  2:4524499/1‑93 (MQ=255)
                                gCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGttt                                                         <  2:6881967/93‑1 (MQ=255)
                                gCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGttt                                                         <  2:6456563/93‑1 (MQ=255)
                                   gagCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTtgct                                                      >  1:5516187/1‑93 (MQ=255)
                                   gagCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTtgct                                                      >  1:7665837/1‑93 (MQ=255)
                                   gagCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTtgct                                                      <  1:6118155/93‑1 (MQ=255)
                                   gagCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGGACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTGtgct                                                      >  2:480164/1‑93 (MQ=255)
                                      cGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTgctgct                                                   >  2:9653740/1‑93 (MQ=255)
                                           cGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACaa                                              >  2:1307521/1‑93 (MQ=255)
                                           cGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACaa                                              <  1:4195837/93‑1 (MQ=255)
                                             gcgcTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAg                                            >  2:7295983/1‑93 (MQ=255)
                                             gcgcTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAg                                            <  1:9056309/93‑1 (MQ=255)
                                               gcTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCa                                          >  1:2919716/1‑93 (MQ=255)
                                               gcTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCa                                          <  1:7158882/93‑1 (MQ=255)
                                                 tGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCAGt                                        <  1:8561581/93‑1 (MQ=255)
                                                 tGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCAGt                                        <  2:7573752/93‑1 (MQ=255)
                                                  gAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCAGTg                                       <  2:5608551/93‑1 (MQ=255)
                                                  gAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCAGTg                                       <  1:5868939/93‑1 (MQ=255)
                                                  gAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCAGTg                                       <  2:9546438/93‑1 (MQ=255)
                                                   aaGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCAGTgg                                      <  1:1290360/93‑1 (MQ=255)
                                                   aaGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCAGTgg                                      >  1:1053284/1‑93 (MQ=255)
                                                    aGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCAGTGGc                                     >  1:9523668/1‑93 (MQ=255)
                                                    aGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCAGTGGc                                     >  1:7549960/1‑93 (MQ=255)
                                                    aGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCAGTGGc                                     <  1:6421410/93‑1 (MQ=255)
                                                    aGGCGTTGGTCTGTATCGTACTGAGTTCCTGGTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGATTGCTGCTTACAAAGCAGTGGc                                     >  2:7036157/1‑93 (MQ=255)
                                                      gCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCAGTGGCTg                                   >  2:5321804/1‑93 (MQ=255)
                                                       cGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCAGTGGCTGa                                  <  2:1087445/93‑1 (MQ=255)
                                                         ttGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGGCGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCAGTGGCTGAAg                                <  2:1008767/93‑1 (MQ=255)
                                                          tGGTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCAGTGGCTGAAGc                               <  1:9378783/93‑1 (MQ=255)
                                                            gTCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCAGTGGCTGAAGCgt                             >  2:2067721/1‑93 (MQ=255)
                                                             tCTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCAGTGGCTGAAGCgtg                            >  2:3978163/1‑93 (MQ=255)
                                                              cTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCAGTGGCTGAAGCgtgt                           <  1:9763671/93‑1 (MQ=255)
                                                              cTGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCAGTGGCTGAAGCgtgt                           <  1:9175362/93‑1 (MQ=255)
                                                               tGTATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCAGTGGCTGAAGCGtgtg                          <  1:9496016/93‑1 (MQ=255)
                                                                 tATCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCAGTGGCTGAAGCGTGTGGc                        <  2:526891/93‑1 (MQ=255)
                                                                  aTCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCAGTGGCTGAAGCGTGTGGCt                       <  2:9749283/93‑1 (MQ=255)
                                                                  aTCGTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCAGTGGCTGAAGCGTGTGGCt                       >  1:1735799/1‑93 (MQ=255)
                                                                     gTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCAGTGGCTGAAGCGTGTGGCTCGc                    >  2:5866490/1‑93 (MQ=255)
                                                                     gTACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCAGTGGCTGAAGCGTGTGGCTCGc                    >  1:5712744/1‑93 (MQ=255)
                                                                      tACTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACACAGCAGTGGCTGAAGCGTGTGGCTCTCa                   >  1:4755657/1‑93 (MQ=255)
                                                                       aCTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCAGTGGCTGAAGCGTGTGGCTCGCaa                  <  1:4550285/93‑1 (MQ=255)
                                                                       aCTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCAACTGAAGAAGAACAGATTGCTGCTTGCAAAGAAGTGGCTGAAGAGTGTGGCTCGcta                  >  1:1608926/1‑91 (MQ=255)
                                                                        cTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCAGTGGCTGAAGCGTGTGGCTCGCAAg                 <  2:6588433/93‑1 (MQ=255)
                                                                        cTGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCAGTGGCTGAAGCGTGTGGCTCGCAAg                 >  1:1792366/1‑93 (MQ=255)
                                                                         tGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCAGTGGCTGAAGCGTGTGGCTCGCAAGc                <  2:5043873/93‑1 (MQ=255)
                                                                         tGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCAGTGGCTGAAGCGTGTGGCTCGCAAGc                >  1:6632911/1‑93 (MQ=255)
                                                                         tGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCAGTGGCTGAAGCGTGTGGCTCGCAAGc                <  2:6235231/93‑1 (MQ=255)
                                                                         tGAGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCAGTGGCTGAAGCGTGTGGCTCGCAAGc                <  1:4085267/93‑1 (MQ=255)
                                                                           aGTTCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCAGTGGCTGAAGCGTGTGGCTCGCAAGCgg              >  1:5967176/1‑93 (MQ=255)
                                                                           aGTGCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCAGTGGCTGAAGCGTGTGGCTCGCAAGCgg              <  1:3207157/93‑1 (MQ=255)
                                                                              tCCTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCAGTGGCTGAAGCGTGTGGCTCGCAAGCGGTTa           <  1:2526300/93‑1 (MQ=255)
                                                                               ccTGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCAGTGGCTGAAGCGTGTGGCTCGCAAGCGGTTAt          <  2:5692000/93‑1 (MQ=255)
                                                                                 tGTTCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCAGTGGCTGAAGCGTGTGGCTCGCAAGCGGTTAtcg        >  2:8251162/1‑93 (MQ=255)
                                                                                    tCATGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCAGTGGCTGAAGCGTGTGGCTCGCAAGCGGTTAtcgttc     <  2:1221382/93‑1 (MQ=255)
                                                                                       tGGGCCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCAGTGGCTGAAGCGTGTGGCTCGCAAGCGGTTATCGTTCGTa  <  2:5956122/93‑1 (MQ=255)
                                                                                          |                                                                                         
CTAACATTGGTACGGTTCGTGACGTTGAAGGTGCAGAGCGTAACGGCGCTGAAGGCGTTGGTCTGTATCGTACTGAGTTCCTGTTCATGGACCGCGACGCACTGCCCACTGAAGAAGAACAGTTTGCTGCTTACAAAGCAGTGGCTGAAGCGTGTGGCTCGCAAGCGGTTATCGTTCGTA  >  NC_000913/2532903‑2533082

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