Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 7,690 T→C G90G (GGA→GGG)  yaaJ ← inner membrane transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009137,6900TC90.6% 619.0 / 45.3 213G90G (GGA→GGG) yaaJinner membrane transport protein
Reads supporting (aligned to +/- strand):  ref base T (12/8);  new base C (103/90);  total (115/98)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.25e-01

Only 100 of 227 total aligned reads displayed.

GTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATGA  >  NC_000913/7600‑7782
                                                                                          |                                                                                            
gTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGcccc                                                                                            >  1:6559253/1‑93 (MQ=255)
gTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGcccc                                                                                            >  1:89426/1‑93 (MQ=255)
gTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGcccc                                                                                            <  1:2596624/93‑1 (MQ=255)
 tCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCa                                                                                           >  2:5464320/1‑93 (MQ=255)
  cACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAg                                                                                          <  1:5676578/93‑1 (MQ=255)
   aCGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAgg                                                                                         <  1:630190/93‑1 (MQ=255)
     gTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTc                                                                                       >  1:1372471/1‑93 (MQ=255)
        cTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCAc                                                                                    <  1:6446842/93‑1 (MQ=255)
         tttATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCAcc                                                                                   >  2:4122878/1‑93 (MQ=255)
          ttATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACcg                                                                                  <  2:6460609/93‑1 (MQ=255)
          ttATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACcg                                                                                  <  2:5523295/93‑1 (MQ=255)
            ataAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACcggc                                                                                <  1:6283632/93‑1 (MQ=255)
             taAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACcggcg                                                                               >  1:473055/1‑93 (MQ=255)
             taAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACcggcg                                                                               >  1:2856185/1‑93 (MQ=255)
               aaGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGt                                                                             >  1:1120601/1‑93 (MQ=255)
                  ttGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAAt                                                                          <  2:3921776/93‑1 (MQ=255)
                  ttGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAAt                                                                          >  1:1222426/1‑93 (MQ=255)
                  ttGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAAt                                                                          >  2:1644060/1‑93 (MQ=255)
                   tgtgCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATa                                                                         <  2:3629539/93‑1 (MQ=255)
                   tgtgCCAGAGAAAATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATa                                                                         <  1:2095632/93‑1 (MQ=255)
                       ccAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCa                                                                     <  2:5245217/93‑1 (MQ=255)
                       ccAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCa                                                                     >  1:6067940/1‑93 (MQ=255)
                        cAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAg                                                                    <  2:3840804/93‑1 (MQ=255)
                         agagAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAgc                                                                   >  1:5951454/1‑93 (MQ=255)
                         agagAACATTAGGCAAACGAGGTCACCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAgc                                                                   <  1:1062494/93‑1 (MQ=255)
                          gagaACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAgcg                                                                  >  2:6518418/1‑93 (MQ=255)
                          gagaACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAgcg                                                                  >  1:2677051/1‑93 (MQ=255)
                           agaACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCCCCGGCGGTAATAGCCAgcgc                                                                 <  1:2956772/93‑1 (MQ=255)
                           agaACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAgcgc                                                                 >  1:4558785/1‑93 (MQ=255)
                            gaACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCa                                                                >  2:5120695/1‑93 (MQ=255)
                              aCATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAAc                                                              <  2:2735794/93‑1 (MQ=255)
                              aCATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAAc                                                              >  2:5968927/1‑93 (MQ=255)
                              aCATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGAGCAAc                                                              >  1:6190124/1‑93 (MQ=255)
                               cATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACg                                                             <  2:1574590/93‑1 (MQ=255)
                                aTTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGc                                                            <  1:320003/93‑1 (MQ=255)
                                 ttCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGcc                                                           >  2:2211094/1‑93 (MQ=255)
                                 ttCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGcc                                                           >  1:1438887/1‑93 (MQ=255)
                                  tCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCg                                                          <  1:4881771/93‑1 (MQ=255)
                                  tCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCg                                                          >  1:955059/1‑93 (MQ=255)
                                  tCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCg                                                          >  1:2556229/1‑93 (MQ=255)
                                   cGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCgg                                                         <  2:2940699/93‑1 (MQ=255)
                                    ggCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCTGCGGGAATAGCCAGCGCAACGCCGGc                                                        >  2:1205409/1‑93 (MQ=255)
                                     gCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGcc                                                       <  1:6637793/93‑1 (MQ=255)
                                     gCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGcc                                                       <  2:5780101/93‑1 (MQ=255)
                                     gCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGca                                                       >  2:902513/1‑92 (MQ=255)
                                       aaaCGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAg                                                     <  2:3224754/93‑1 (MQ=255)
                                       aaaCGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAg                                                     <  2:4938905/93‑1 (MQ=255)
                                        aaCGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCA‑GTCCACAGGCGGTAATAGCCAGCGCAACGCCGGCCAGGt                                                   >  1:3469600/1‑93 (MQ=255)
                                          cGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGtt                                                  >  2:6085048/1‑93 (MQ=255)
                                           gAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTg                                                 <  2:434189/93‑1 (MQ=255)
                                            aGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGc                                                >  1:3933168/1‑93 (MQ=255)
                                            aGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGc                                                >  2:957760/1‑93 (MQ=255)
                                            aGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGc                                                <  2:6593515/93‑1 (MQ=255)
                                            aGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCACAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGc                                                >  1:537768/1‑93 (MQ=255)
                                             ggTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGcc                                               >  1:2138799/1‑93 (MQ=255)
                                                  ccATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTAc                                          >  1:4143301/1‑93 (MQ=255)
                                                  ccATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTAc                                          >  2:194858/1‑93 (MQ=255)
                                                   cATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTAcc                                         <  2:2615226/93‑1 (MQ=255)
                                                    aTGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTAccc                                        >  2:3133063/1‑93 (MQ=255)
                                                    aTGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTAccc                                        >  2:4044737/1‑93 (MQ=255)
                                                     tGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCa                                       <  2:1884231/93‑1 (MQ=255)
                                                         gATAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCAcgcg                                   >  2:2602/1‑93 (MQ=255)
                                                          aTAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCAcgcgc                                  <  2:5341429/93‑1 (MQ=255)
                                                           tAAACGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCAcgcgcg                                 >  1:3598230/1‑93 (MQ=255)
                                                            aaaCGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACgcgcgc                                >  1:491691/1‑93 (MQ=255)
                                                            aaaCGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACgcgcgc                                >  2:3345866/1‑93 (MQ=255)
                                                             aaCGCGGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGcgcc                               <  2:5098678/93‑1 (MQ=255)
                                                                 cgGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCaa                           <  1:5847742/93‑1 (MQ=255)
                                                                 cgGCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCaa                           <  2:3543440/93‑1 (MQ=255)
                                                                  ggCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAg                          >  1:6334006/1‑93 (MQ=255)
                                                                   gCAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGa                         >  1:4201541/1‑93 (MQ=255)
                                                                    cAACCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGAc                        >  1:2265339/1‑93 (MQ=255)
                                                                     aaCCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACt                       <  1:3693782/93‑1 (MQ=255)
                                                                     aaCCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACt                       >  2:4566049/1‑93 (MQ=255)
                                                                      aCCCACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTg                      >  1:736685/1‑93 (MQ=255)
                                                                       cccATATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTgg                     <  1:5788619/93‑1 (MQ=255)
                                                                       cccACATCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTgg                     <  1:4403349/93‑1 (MQ=255)
                                                                         cacaTCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTa                   <  2:4655837/93‑1 (MQ=255)
                                                                         cacaTCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTa                   >  2:1508286/1‑93 (MQ=255)
                                                                         cacaTCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTa                   <  1:1860612/93‑1 (MQ=255)
                                                                         cacaTCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTa                   <  2:3532430/93‑1 (MQ=255)
                                                                           caTCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTaca                 >  1:3055640/1‑93 (MQ=255)
                                                                           caTCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTaca                 >  2:4419236/1‑93 (MQ=255)
                                                                            aTCCAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTacac                <  2:764981/93‑1 (MQ=255)
                                                                               cAATAGACGCCCCCAGCGCCACCGGCGGTAATAGCCGTGGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAAt             <  2:1767445/93‑1 (MQ=255)
                                                                               cAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAAt             >  2:297888/1‑93 (MQ=255)
                                                                               cAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCCAGACTGGTACACAAt             >  1:3670185/1‑93 (MQ=255)
                                                                               cAAAAGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAAt             >  2:483789/1‑93 (MQ=255)
                                                                                aaaaGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATg            >  2:1244719/1‑93 (MQ=255)
                                                                                aaaaGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATg            >  2:5688176/1‑93 (MQ=255)
                                                                                aaaaGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATg            <  2:11724/93‑1 (MQ=255)
                                                                                aaaaGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATg            >  1:3906029/1‑93 (MQ=255)
                                                                                aaaaGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATg            >  2:1298112/1‑93 (MQ=255)
                                                                                 aaaGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGa           >  1:1643980/1‑93 (MQ=255)
                                                                                  aaGACGGCTCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGAc          <  1:1180490/93‑1 (MQ=255)
                                                                                  aaGACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGAc          <  2:6056321/93‑1 (MQ=255)
                                                                                    gACGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTAGACAATGACTg        >  1:3476146/1‑93 (MQ=255)
                                                                                     aCGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGa       >  1:2234886/1‑93 (MQ=255)
                                                                                      cGGCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTTCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGaa      <  2:6411742/93‑1 (MQ=255)
                                                                                       ggCCCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGaaa     <  2:4991994/93‑1 (MQ=255)
                                                                                          |                                                                                            
GTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATGA  >  NC_000913/7600‑7782

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