Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,748,779 (A)6→7 coding (310/720 nt) b1671 ← putative oxidoreductase Fe‑S subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,748,7731.A99.3% 522.6 / ‑3.5 147coding (316/720 nt)b1671putative oxidoreductase Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base A (65/81);  total (66/81)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

Only 100 of 156 total aligned reads displayed.

TACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCG‑AAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAATGTGTGCTATCGATAAACGGCTTCCCTGAGCGGGGACGTGG  >  NC_000913/1748681‑1748858
                                                                                              |                                                                                      
tACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCg                                                                                         >  2:79392/1‑93 (MQ=255)
 aCCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCG‑a                                                                                       >  2:618210/1‑93 (MQ=255)
      cACGATCCCCTGCTCATCACGCCACGACGCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCG‑aaaaa‑a                                                                                 >  1:1936107/1‑93 (MQ=255)
         gATCCCCTGCTCATCACGCCACGACGCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCG‑aaaaa‑aagt                                                                              <  2:2637692/93‑4 (MQ=255)
          aTCCCCTGCTCATCACGCCACGACGCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCG‑aaaaa‑aagtg                                                                             >  1:3465827/1‑89 (MQ=255)
          aTCCCCTGCTCATCACGCCACGACGCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCG‑aaaaa‑aagtg                                                                             >  2:1127013/1‑89 (MQ=255)
             cccTGCTCCTCACGCCCCGACGCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGGCTGGGAG‑AAAAA‑AGCGaca                                                                           >  2:1600839/1‑91 (MQ=255)
             cccTGCTCATCACGCCACGACGCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATa                                                                           >  1:4672931/1‑93 (MQ=255)
              ccTGCTCATCACGCCACGACGCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATAc                                                                          >  2:582087/1‑93 (MQ=255)
               cTGCTCATCACGCCACGACGCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACt                                                                         >  1:2041091/1‑93 (MQ=255)
                cgCTCATCACGCCACGACGCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTg                                                                        <  1:1286401/92‑1 (MQ=255)
                 gCTCATCACGCCACGACGCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGc                                                                       <  1:1333634/93‑1 (MQ=255)
                   tcatcaCGCCACGACGCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGt                                                                     >  2:2841103/1‑93 (MQ=255)
                    catcaCGCCACGACGCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGtt                                                                    <  2:611244/93‑1 (MQ=255)
                    catcaCGCCACGACGCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGtt                                                                    >  1:2724833/1‑93 (MQ=255)
                     atcaCGCCACGACGCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGttt                                                                   >  1:46667/1‑93 (MQ=255)
                       caCGCCACGACGCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCg                                                                 <  1:65209/93‑1 (MQ=255)
                        aCGCCACGACGCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGt                                                                <  1:1746837/93‑1 (MQ=255)
                         cGCCACGACGCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGtt                                                               <  2:4107432/93‑1 (MQ=255)
                         cGCCACGACGCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGtt                                                               <  2:357402/93‑1 (MQ=255)
                         cGCCACGACGCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGtt                                                               <  2:914571/93‑1 (MQ=255)
                          gCCACGACGCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTg                                                              >  2:2087967/1‑93 (MQ=255)
                          gCCACGACGCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTg                                                              >  2:4026923/1‑93 (MQ=255)
                          gCCACGACGCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTg                                                              <  2:4467245/93‑1 (MQ=255)
                           ccACGACGCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGt                                                             <  2:2650030/93‑1 (MQ=255)
                            cACGACGCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTc                                                            <  2:1924933/93‑1 (MQ=255)
                            cACGACGCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTc                                                            <  2:2427792/93‑1 (MQ=255)
                             acgacgCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCg                                                           >  2:1745973/1‑93 (MQ=255)
                                acgCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGttt                                                        >  1:938155/1‑92 (MQ=255)
                                acgCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTa                                                        >  1:3885091/1‑93 (MQ=255)
                                acgCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTa                                                        >  2:3149377/1‑93 (MQ=255)
                                 cgCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTAt                                                       <  1:2507552/93‑1 (MQ=255)
                                 cgCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTAt                                                       <  1:4333796/93‑1 (MQ=255)
                                  gCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATc                                                      <  2:4727091/93‑1 (MQ=255)
                                  gCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATc                                                      <  1:2347608/93‑1 (MQ=255)
                                  gCCCCCGTCGGGCAAACGTCGATGAACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATc                                                      >  1:563899/1‑93 (MQ=255)
                                   cccccGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCg                                                     >  1:2527769/1‑93 (MQ=255)
                                     cccGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGt                                                   <  1:3616267/93‑1 (MQ=255)
                                      ccGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTg                                                  >  2:2741162/1‑93 (MQ=255)
                                      ccGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTg                                                  >  1:1120047/1‑93 (MQ=255)
                                      ccGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTg                                                  >  2:733962/1‑93 (MQ=255)
                                      ccGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTg                                                  >  1:384314/1‑93 (MQ=255)
                                      ccGTCGGGCAAACGTCGAGGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTg                                                  <  2:1468822/93‑1 (MQ=255)
                                       cGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTgc                                                 >  1:1237077/1‑92 (MQ=255)
                                        gTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGAc                                                <  2:2703245/93‑1 (MQ=255)
                                        gTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGAc                                                >  2:1698917/1‑93 (MQ=255)
                                        ctCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTTATACTGCGTTTCGTTGTCGTTATCGGTGAc                                                <  2:3461183/92‑1 (MQ=255)
                                          cGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTg                                              >  2:4455375/1‑93 (MQ=255)
                                           gggCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTgg                                             >  1:1261041/1‑93 (MQ=255)
                                           gggCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTgg                                             >  1:2776668/1‑93 (MQ=255)
                                            ggCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGCCATGACT‑GGCG‑AAAAATAGAGATACTGCGTTTCGTTGTCGTTATAGGTGTCTGGa                                            >  2:3769511/1‑93 (MQ=255)
                                             gCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGaa                                           <  1:3055413/93‑1 (MQ=255)
                                             gCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGaa                                           <  2:1453095/93‑1 (MQ=255)
                                              cAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGGCTGGAAt                                          >  1:1241834/1‑93 (MQ=255)
                                              cAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTAGTTATCGGTGACTGGAAt                                          <  1:2274664/93‑1 (MQ=255)
                                                 aCGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAAtgtg                                       >  1:2124765/1‑93 (MQ=255)
                                                 aCGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAAtgtg                                       >  1:1320708/1‑93 (MQ=255)
                                                 aCGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAAtgtg                                       <  2:2061649/93‑1 (MQ=255)
                                                  cGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAAtgtgt                                      <  2:4269561/93‑1 (MQ=255)
                                                  cGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAAtgtgt                                      <  2:2390354/93‑1 (MQ=255)
                                                  cGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAAtgtgt                                      >  1:1842912/1‑93 (MQ=255)
                                                     cGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAATGTGTGCt                                   <  2:1087294/93‑1 (MQ=255)
                                                      gATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAATGTGTGCTa                                  >  2:3303594/1‑93 (MQ=255)
                                                       aTGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAATGTGTGCTAt                                 <  2:1161366/93‑1 (MQ=255)
                                                        tGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAATGTGTGCTATc                                <  2:446759/93‑1 (MQ=255)
                                                             ggAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAATGTGTGCTATCGATaa                           >  1:4083155/1‑93 (MQ=255)
                                                             ggAGCATCTTCGCAGTACTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAATGTGTGCTATCGATaa                           >  1:423190/1‑93 (MQ=255)
                                                              gAGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAATGTGTGCTATCGATaaa                          <  1:544460/93‑1 (MQ=255)
                                                               aGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAATGTGTGCTATCGATAAAc                         <  2:4459606/93‑1 (MQ=255)
                                                               aGCATCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAATGTGTGCTATCGATAAAc                         <  2:3557843/93‑1 (MQ=255)
                                                                  aTCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAATGTGTGCTATCGATAAACGGc                      >  1:1763954/1‑93 (MQ=255)
                                                                  aTCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAATGGGTGCTATCGATAAACGGc                      >  1:3521922/1‑93 (MQ=255)
                                                                   tCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAATGTGTGCTATCGATAAACGGCt                     >  1:3554001/1‑93 (MQ=255)
                                                                   tCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAATGTGTGCTATCGATAAACGGCt                     >  1:2921346/1‑93 (MQ=255)
                                                                   tCTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAATGTGTGCTATCGATAAACGGCt                     <  1:1799917/93‑1 (MQ=255)
                                                                    cTTCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAATGTGTGCTATCGATAAACGGCtt                    <  2:964966/93‑1 (MQ=255)
                                                                      tCGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAATGTGTGCTATCGATAAACGGCTTcc                  >  1:1660768/1‑93 (MQ=255)
                                                                       cGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAATGTGTGCTATCGATAAACGGCTTccc                 <  1:53376/93‑1 (MQ=255)
                                                                       cGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAATGTGTGCTATCGATAAACGGCTTccc                 <  1:3284050/93‑1 (MQ=255)
                                                                       cGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAATGTGTGCTATCGATAAACGGCTTccc                 <  2:3757029/93‑1 (MQ=255)
                                                                       cGCAGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAATGTGTGCTATCGATAAACGGCTTccc                 <  1:1816205/93‑1 (MQ=255)
                                                                        gCAGTGCTGACATGACA‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAATGTGTGCTATCGATAAACGGCTTCCCt                >  2:3131664/1‑93 (MQ=255)
                                                                          aGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAATGTGTGCTATCGATAAACGGCTTCCCTGa              <  1:2211318/93‑1 (MQ=255)
                                                                          aGTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAATGTGTGCTATCGATAAACGGCTTCCCTGa              >  2:1885728/1‑93 (MQ=255)
                                                                           gTGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAATGTGTGCTATCGATAAACGTCTTCCCTGAg             >  1:707314/1‑93 (MQ=255)
                                                                            tGCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAATGTGTGCTATCGATAAACGGCTTCCCTGAGc            <  1:925464/93‑1 (MQ=255)
                                                                             gCTGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAATGTGTGCTATCGATAAACGGCTTCCCTGAGCg           <  1:2647326/93‑1 (MQ=255)
                                                                               tGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAATGTGTGCTATCGATAAACGGCTTCCCTGAGCggg         >  2:143365/1‑93 (MQ=255)
                                                                               tGACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTAGTTATCGGTGACTGGAATGTGTGCTATCGATAAACGGCTTCCCTGAGCggg         <  2:2755449/93‑1 (MQ=255)
                                                                                gACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAATGTGTGCTATCGATAAACGGCTTCCCTGAGCgggg        <  2:965215/93‑1 (MQ=255)
                                                                                gACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAATGTGTGCTATCGATAAACGGCTTCCCTGAGCgggg        <  2:335804/93‑1 (MQ=255)
                                                                                gACATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAATGTGTGCTATCGATAAACGGCTTCCCTGAGCgggg        <  1:2539914/93‑1 (MQ=255)
                                                                                 aCATGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAATGTGTGCTATCGATAAACGGCTTCCCTGAGCGGGGa       <  1:119425/93‑1 (MQ=255)
                                                                                   aTGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAATGTGTGCTATCGATAAACGGCTTCCCTGAGCGGGGACg     >  2:2333656/1‑93 (MQ=255)
                                                                                   aTGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAATGTGTGCTATCGATAAACGGCTTCCCTGAGCGGGGACg     <  1:2241749/93‑1 (MQ=255)
                                                                                    tGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAATGTGTGCTATCGATAAACGGCTTCCCTGAGCGGGGACgg    >  1:4369197/1‑92 (MQ=255)
                                                                                    tGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAATGTGTGCTATCGATAAACGGCTTCCCTGAGCGGGGACGt    >  2:4697232/1‑93 (MQ=255)
                                                                                    tGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAATGTGTGCTATCGATAAACGGCTTCCCTGAGCGGGGACGt    <  1:4093099/93‑1 (MQ=255)
                                                                                    tGACT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAATGTGTGCTATCGATAAACGGCTTCCCTGAGCGGGGACGt    >  1:2908844/1‑93 (MQ=255)
                                                                                     gtcT‑GGCGAAAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAATGGGGGCTATCGATAAACGGCTTCCCTGAGCGGGGACGTg   <  1:3579592/91‑1 (MQ=255)
                                                                                              |                                                                                      
TACCCGCACGATCCCCTGCTCATCACGCCACGACGCCCCCGTCGGGCAAACGTCGATGCACGGAGCATCTTCGCAGTGCTGACATGACT‑GGCG‑AAAAA‑AGTGATACTGCGTTTCGTTGTCGTTATCGGTGACTGGAATGTGTGCTATCGATAAACGGCTTCCCTGAGCGGGGACGTGG  >  NC_000913/1748681‑1748858

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: ‑