Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 7,646 G→A T105I (ACC→ATC)  yaaJ ← inner membrane transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009137,6460GA99.3% 957.0 / ‑4.6 271T105I (ACC→ATC) yaaJinner membrane transport protein
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (150/119);  total (152/119)
Fisher's exact test for biased strand distribution p-value = 5.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.84e-01

Only 100 of 281 total aligned reads displayed.

CGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTG  >  NC_000913/7556‑7735
                                                                                          |                                                                                         
   gcCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGcc                                                                                      <  1:2284589/93‑1 (MQ=255)
    ccATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCa                                                                                     <  2:6850871/93‑1 (MQ=255)
    ccATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCa                                                                                     <  2:128214/93‑1 (MQ=255)
    ccATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCa                                                                                     <  1:3489259/93‑1 (MQ=255)
    ccATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCa                                                                                     <  1:7643707/93‑1 (MQ=255)
     cATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCAt                                                                                    >  2:7244837/1‑93 (MQ=255)
     cATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCAt                                                                                    <  2:1668359/93‑1 (MQ=255)
       tataCCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGc                                                                                  >  1:6950640/1‑93 (MQ=255)
       tataCCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGc                                                                                  <  2:2134841/93‑1 (MQ=255)
       tataCCATGCCGGTCCGCCACGAAACTGCCCATTGACGACACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATAGCCATGc                                                                                  >  2:6115383/1‑93 (MQ=255)
        ataCCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGcc                                                                                 >  1:2992644/1‑93 (MQ=255)
         taCCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCg                                                                                >  2:1341680/1‑93 (MQ=255)
         taCCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCg                                                                                >  2:587601/1‑93 (MQ=255)
          tccACGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGa                                                                               <  2:5696561/92‑1 (MQ=255)
             aTGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATaa                                                                            >  1:4852650/1‑93 (MQ=255)
               gCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAAc                                                                          <  2:4261986/93‑1 (MQ=255)
               gCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAAc                                                                          <  1:1974949/93‑1 (MQ=255)
                 cGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAAcgc                                                                        <  1:6798386/93‑1 (MQ=255)
                  ggTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGACCGCCATGCCGATAAAcgcg                                                                       >  1:2796801/1‑93 (MQ=255)
                   gcccGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCgg                                                                      <  1:3795372/91‑1 (MQ=255)
                   gTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCgg                                                                      <  1:6240233/93‑1 (MQ=255)
                   gTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCgg                                                                      <  1:7421011/93‑1 (MQ=255)
                    tCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGc                                                                     >  1:2635238/1‑93 (MQ=255)
                    tCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGc                                                                     >  1:6729077/1‑93 (MQ=255)
                        ccACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAAcc                                                                 <  1:3698387/93‑1 (MQ=255)
                          aCGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCca                                                               <  2:1185650/93‑1 (MQ=255)
                          aCGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCca                                                               >  1:7543606/1‑93 (MQ=255)
                          aCGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCca                                                               >  1:2425920/1‑93 (MQ=255)
                            gAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCcaca                                                             <  2:4459017/93‑1 (MQ=255)
                            gAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCcaca                                                             <  2:7775913/93‑1 (MQ=255)
                             aaaCTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACAt                                                            >  2:7376949/1‑93 (MQ=255)
                              aaCTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATc                                                           <  1:1273083/93‑1 (MQ=255)
                               aCTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATcc                                                          >  2:7487268/1‑93 (MQ=255)
                               aCTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATcc                                                          <  1:4125318/93‑1 (MQ=255)
                               aCTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATcc                                                          >  1:6830742/1‑93 (MQ=255)
                               aCTGCCAATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATCCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATcc                                                          <  2:257550/93‑1 (MQ=255)
                                 tGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCaa                                                        >  2:5856690/1‑93 (MQ=255)
                                  gCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCaca                                                       >  1:1917723/1‑91 (MQ=255)
                                   cccATTGACGGCACGTTCTTGATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCaaaa                                                      <  1:4243813/93‑1 (MQ=255)
                                     cATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGa                                                    <  1:3081962/93‑1 (MQ=255)
                                      aTTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGAc                                                   <  2:2250209/93‑1 (MQ=255)
                                       ttGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACg                                                  >  2:5184674/1‑93 (MQ=255)
                                       ttGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACg                                                  >  1:533203/1‑93 (MQ=255)
                                       ttGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACg                                                  <  2:1538545/93‑1 (MQ=255)
                                       ttGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACg                                                  >  1:4049706/1‑93 (MQ=255)
                                        tGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACgg                                                 <  1:2084891/93‑1 (MQ=255)
                                        tGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACgg                                                 <  2:4159756/93‑1 (MQ=255)
                                        tGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACgg                                                 >  2:7331129/1‑93 (MQ=255)
                                        tGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACgg                                                 >  1:1924421/1‑93 (MQ=255)
                                         gACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGc                                                <  2:3831105/93‑1 (MQ=255)
                                         gACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGc                                                >  1:6437459/1‑93 (MQ=255)
                                          aCGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCt                                               >  1:7822404/1‑93 (MQ=255)
                                            gTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTcc                                             >  2:6290548/1‑93 (MQ=255)
                                             tCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCa                                            <  1:913250/93‑1 (MQ=255)
                                             tCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCa                                            >  1:4953304/1‑93 (MQ=255)
                                             tCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCa                                            <  2:1082681/93‑1 (MQ=255)
                                             tCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCa                                            <  1:5826583/93‑1 (MQ=255)
                                              cACGTTCTTTATAAAGTTGTGCCAGAGAATATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAg                                           <  1:2921399/93‑1 (MQ=255)
                                               aCGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAgg                                          <  2:4524143/93‑1 (MQ=255)
                                                cGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGt                                         >  2:7624243/1‑93 (MQ=255)
                                                 gTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTc                                        >  1:4437496/1‑93 (MQ=255)
                                                  ttCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAAGCCACATCCAAAAGACGGCTCCAGGTca                                       >  2:5372857/1‑92 (MQ=255)
                                                  ttCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTcc                                       >  1:1098490/1‑93 (MQ=255)
                                                   tCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCa                                      >  1:3224433/1‑93 (MQ=255)
                                                   tCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCa                                      >  1:5721379/1‑93 (MQ=255)
                                                    cTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCAc                                     >  2:3602302/1‑93 (MQ=255)
                                                     tttATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCAcc                                    >  2:5761644/1‑93 (MQ=255)
                                                       tataAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACcgg                                  >  2:234762/1‑93 (MQ=255)
                                                       tataAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACcgg                                  >  1:1577236/1‑93 (MQ=255)
                                                       tataAAGTTGTGCCAGAGAACATACGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACcgg                                  >  1:1852499/1‑93 (MQ=255)
                                                        ataAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCAAGGCCAACcgcg                                 >  1:6461114/1‑91 (MQ=255)
                                                         taAAGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACcggcg                                >  1:7903372/1‑93 (MQ=255)
                                                          aaaGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCCTGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACcggcgg                               <  2:2729839/93‑1 (MQ=255)
                                                           aaGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGt                              >  2:7589963/1‑93 (MQ=255)
                                                            aGTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTa                             >  1:4694314/1‑93 (MQ=255)
                                                             gTTGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTaa                            <  1:730895/93‑1 (MQ=255)
                                                              ttGTGCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAAt                           <  2:6808699/93‑1 (MQ=255)
                                                                  gCCAGAGAACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGcc                       >  2:3982899/1‑93 (MQ=255)
                                                                      gagaACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAgcg                   >  2:7132798/1‑93 (MQ=255)
                                                                      gagaACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAgcg                   >  1:2708093/1‑93 (MQ=255)
                                                                       agaACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAgcgc                  >  1:7321831/1‑93 (MQ=255)
                                                                       agaACATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAgcgc                  <  1:6079736/93‑1 (MQ=255)
                                                                          aCATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAGCGCCAc               >  2:3731819/1‑93 (MQ=255)
                                                                          aCATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAGCGCAAc               <  2:3414896/93‑1 (MQ=255)
                                                                          aCATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAGCGCAAc               <  2:1600826/93‑1 (MQ=255)
                                                                          aCATTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGAAGGCTCCAGGTCGACCGGTGGTAATAGCCAGGGCAAc               >  2:1736179/1‑93 (MQ=255)
                                                                            aTTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGc             <  2:6333021/93‑1 (MQ=255)
                                                                            aTTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGc             >  2:727848/1‑93 (MQ=255)
                                                                            aTTCGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGc             >  2:5541146/1‑93 (MQ=255)
                                                                               cGGCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCgg          >  2:3573945/1‑93 (MQ=255)
                                                                                 gCAAACGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGcc        >  2:2460070/1‑93 (MQ=255)
                                                                                   aaaCGAGATCGCCATGCCTATAAACGCGTCAACCCACATCCAAAAGACGGCTCCAGGTCCACCTGCGGTAATAGCCAGCGCAACGCCGGCCAg      >  2:6011432/1‑93 (MQ=255)
                                                                                   aaaCGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAg      <  2:1869732/93‑1 (MQ=255)
                                                                                    aaCGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCCgg     >  2:296103/1‑93 (MQ=255)
                                                                                    aaCGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAgg     >  1:1638619/1‑93 (MQ=255)
                                                                                    aaCGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAgg     <  1:5248654/93‑1 (MQ=255)
                                                                                    aaCGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAgg     <  1:5655992/93‑1 (MQ=255)
                                                                                    aaCGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCACCCGGGGGTAATAGCCAGCGCAACGCCGGCCAgg     >  1:5972733/1‑93 (MQ=255)
                                                                                     aCGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGt    <  2:487102/93‑1 (MQ=255)
                                                                                      cGAGATCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGtt   <  2:121604/93‑1 (MQ=255)
                                                                                          |                                                                                         
CGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTG  >  NC_000913/7556‑7735

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