Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 905,689 G→A R108R (CGC→CGT)  b0868 ← putative nucleotide di‑P‑sugar epimerase or dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913905,6890GA87.3% 630.1 / 54.6 222R108R (CGC→CGT) b0868putative nucleotide di‑P‑sugar epimerase or dehydratase
Reads supporting (aligned to +/- strand):  ref base G (18/10);  new base A (96/97);  total (114/108)
Fisher's exact test for biased strand distribution p-value = 1.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01

Only 100 of 231 total aligned reads displayed.

TATCGCGATGGTGGTGATAATCGAAGTACAGGGAGGGGGAAGAGATATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGGCGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTAAAGCTGGAGCAGTGCCACAGCGTATCAAT  >  NC_000913/905598‑905778
                                                                                           |                                                                                         
  tCGCGATGGTGGTGATAATCGAAGTACAGGGAGGGGGAAGAGATATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAgacgg                                                                                        >  2:1523219/1‑93 (MQ=255)
  tCGCGATGGTGGTGATAATCGAAGTACAGGGAGGGGGAAGAGATATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAgacgg                                                                                        >  2:1895980/1‑93 (MQ=255)
  tCGCGATGGTGGTGATAATCGAAGTACAGGGAGGGGGAAGAGATATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAgacgg                                                                                        <  2:3204456/93‑1 (MQ=255)
  tCGCGATGGTGGTGATAATCGAAGTACAGGGAGGGGGAAGAGATATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAgacgg                                                                                        >  1:522203/1‑93 (MQ=255)
   cgcgATGGTGGTGATAATCGAAGTACAGGGAGGGGGAAGAGATATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAgacggc                                                                                       <  1:1963793/93‑1 (MQ=255)
    gcgATGGTGGTGATAATCGAAGTACAGGGAGGGGGAAGAGATATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGacggcg                                                                                      <  1:129183/93‑1 (MQ=255)
    gcgATGGTGGTGATAATCGAAGCACAGGGAGGGGGAAGAGATATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGgcggcg                                                                                      >  1:670791/1‑93 (MQ=255)
     cgATGGTGGTGATAATCGAAGTACAGGGAGGGGGAAGAGATATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGa                                                                                     <  2:223979/93‑1 (MQ=255)
       aTGGTGGTGATAATCGAAGTACAGGGAGGGGGAAGAGATATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGt                                                                                   >  1:2685991/1‑93 (MQ=255)
        tggtggTGATAATCGAAGTACAGGGAGGGGGAAGAGATATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTg                                                                                  >  1:4004492/1‑93 (MQ=255)
          gtggtgATAATCGAAGTACAGGGAGGGGGAAGAGATATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCTCCCAGGCGGCGAGTGgc                                                                                <  1:2429036/93‑1 (MQ=255)
          gtggtgATAATCGAAGTACAGGGAGGGGGAAGAGATATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGgc                                                                                <  1:4103125/93‑1 (MQ=255)
           tggtgATAATCGAAGTACAGGGAGGGGGAAGAGATATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGgcg                                                                               >  1:2296098/1‑93 (MQ=255)
            ggtgATAATCGAAGTACAGGGAGGGGGAAGAGATATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGgcgc                                                                              <  1:6691308/93‑1 (MQ=255)
               gATAATCGAAGTACAGGGAGGGGGAAGAGATATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGaa                                                                           <  2:3411549/93‑1 (MQ=255)
                aTAATCGAAGTACAGGGAGGGGGAAGAGATATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAAc                                                                          <  2:6261381/93‑1 (MQ=255)
                   aTCGAAGTACAGGGAGGGGGAAGAGATATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGtt                                                                       >  2:6096566/1‑93 (MQ=255)
                   aTCGAAGTACAGGGAGGGGGAAGAGATATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGtt                                                                       >  2:6019521/1‑93 (MQ=255)
                    tCGAAGTACAGGGAGGGGGAAGAGATATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTa                                                                      <  2:2862685/93‑1 (MQ=255)
                    tCGAAGCACAGGGAGGGGGAAGAGATATGAATAAAGTTACGCACCCCCCAGGCGACAGCCCATTCACCCAGGCGGCGAGTGGCGCGAACGTTa                                                                      <  1:4245263/93‑1 (MQ=255)
                    tCGAAGCACAGGGAGGGGGAAGAGATATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGGCGGCGAGTGGCGCGAACGTTa                                                                      >  1:2758580/1‑93 (MQ=255)
                    tCGAAGCACAGGGAGGGGGAAGAGATATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGGCGGCGAGTGGCGCGAACGTTa                                                                      >  1:2388318/1‑93 (MQ=255)
                       aaGTACAGGGAGGGGGAAGAGATATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGcc                                                                   >  2:4172399/1‑93 (MQ=255)
                       aaGTACAGGGAGGGGGAAGAGATATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGcc                                                                   >  1:2012572/1‑93 (MQ=255)
                       aaGCACAGGGAGGGGGAAGAGATATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGGCGGCGAGTGGCGCGAACGTTAGcc                                                                   <  2:4617132/93‑1 (MQ=255)
                        aGTACAGGGAGGGGGAAGAGATATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCa                                                                  <  1:5022612/93‑1 (MQ=255)
                         gTACAGGGAGGGGGAAGAGATATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAg                                                                 <  1:4319772/93‑1 (MQ=255)
                          tACAGGGAGGGGGAAGAGATATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGa                                                                <  1:7427102/93‑1 (MQ=255)
                           aCAGGGAGGGGGAAGAGATATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGAt                                                               <  1:5565654/93‑1 (MQ=255)
                               ggagggGGAAGAGATATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGaa                                                           <  1:1700796/93‑1 (MQ=255)
                                 agggGGAAGAGATATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAg                                                         <  2:5360245/93‑1 (MQ=255)
                                   ggggAAGAGATATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGcc                                                       <  2:875223/93‑1 (MQ=255)
                                   ggggAAGAGATATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGcc                                                       >  1:3703615/1‑93 (MQ=255)
                                    gggAAGAGATATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCt                                                      <  1:4663467/93‑1 (MQ=255)
                                      gAAGAGATATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCtgt                                                    <  2:1844857/93‑1 (MQ=255)
                                      gAAGAGATATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCtgt                                                    <  2:3760452/93‑1 (MQ=255)
                                       aaGAGATATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCtgtt                                                   <  2:6688100/93‑1 (MQ=255)
                                       aaGAGATATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCtgtt                                                   <  2:63311/93‑1 (MQ=255)
                                         gagaTATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCtgttgt                                                 >  1:5493151/1‑93 (MQ=255)
                                         gagaTATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCtgttgt                                                 >  1:2136774/1‑93 (MQ=255)
                                          agaTATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGGCGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTtgtg                                                >  2:1077250/1‑93 (MQ=255)
                                           gaTATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTgtgt                                               >  1:6123280/1‑93 (MQ=255)
                                            atatGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGGCGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTc                                              >  1:6542925/1‑93 (MQ=255)
                                            atatGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTc                                              >  2:2043577/1‑93 (MQ=255)
                                             tatGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTcc                                             >  2:4939724/1‑93 (MQ=255)
                                             tatGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTcc                                             >  1:5880391/1‑93 (MQ=255)
                                               tGAATAAAGTTACGCACACCCCAGGCAACACCCCATTCACCCAGACGCCGAGGGCCGCGAACGTAGGCCAGACGGAAAGCCTGTTGTGCcccc                                           >  1:4780889/1‑93 (MQ=255)
                                                  aTAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAg                                        >  1:4938505/1‑93 (MQ=255)
                                                  aTAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAg                                        >  1:1415333/1‑93 (MQ=255)
                                                  aTAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAg                                        >  1:4208492/1‑93 (MQ=255)
                                                    aaaGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAggg                                      <  1:7175886/93‑1 (MQ=255)
                                                      aGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTg                                    >  1:6027897/1‑93 (MQ=255)
                                                       gTTACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGa                                   >  2:3852922/1‑93 (MQ=255)
                                                       gTTACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGa                                   <  2:1969707/93‑1 (MQ=255)
                                                        ttACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAg                                  <  1:2490100/93‑1 (MQ=255)
                                                        ttACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAg                                  >  2:3443882/1‑93 (MQ=255)
                                                        ttACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAg                                  >  2:5233049/1‑93 (MQ=255)
                                                         tACGCACACCCCAGGCGACAGCCCATTCACCCAGGCGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAgg                                 >  1:5238709/1‑93 (MQ=255)
                                                         tACGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAgg                                 <  1:1744024/93‑1 (MQ=255)
                                                          aCGCACACCCCAGGCGACAGCCCATTCACCCAGGCGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGt                                <  2:6646299/93‑1 (MQ=255)
                                                           cGCACACCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTa                               >  1:6142518/1‑93 (MQ=255)
                                                            gCACACCCCAGGCGACAGCCCATTCACCCAGGCGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTaa                              <  1:859138/93‑1 (MQ=255)
                                                             cacaCCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTaaa                             <  2:5115747/93‑1 (MQ=255)
                                                              acacCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTAAAg                            >  1:2180668/1‑93 (MQ=255)
                                                              acacCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTAAAg                            >  2:2277075/1‑93 (MQ=255)
                                                                acCCCAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTAAAGCt                          <  2:1443058/93‑1 (MQ=255)
                                                                 ccccAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTAAAGCTg                         <  1:6932101/93‑1 (MQ=255)
                                                                  cccAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTAAAGCTgg                        >  2:2683293/1‑93 (MQ=255)
                                                                  cccAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTAAAGCTgg                        >  2:6043720/1‑93 (MQ=255)
                                                                   ccAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTAAAGCTGGa                       >  1:3066153/1‑93 (MQ=255)
                                                                    cAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTAAAGCTGGAg                      <  2:4009955/93‑1 (MQ=255)
                                                                    cAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTAAAGCTGGAg                      >  1:3466983/1‑93 (MQ=255)
                                                                    cAGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTAAAGCTGGAg                      >  1:5118563/1‑93 (MQ=255)
                                                                     aGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTAAAGCTGGAGc                     >  2:6915761/1‑93 (MQ=255)
                                                                     aGGCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTAAAGCTGGAGc                     <  1:4508394/93‑1 (MQ=255)
                                                                      ggCGACAGCCCATTCACCCAGGCGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTAAAGCTGGAGCa                    >  2:4231118/1‑93 (MQ=255)
                                                                      ggCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTAAAGCTGGAGCa                    >  2:7330845/1‑93 (MQ=255)
                                                                      ggCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTAAAGCTGGAGCa                    <  2:7209194/93‑1 (MQ=255)
                                                                      ggCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTAAAGCTGGAGCa                    >  2:4471636/1‑93 (MQ=255)
                                                                      ggCGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTAAAGCTGGAGCa                    >  2:5821212/1‑93 (MQ=255)
                                                                        cGACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTAAAGCTGGAGCAGt                  <  1:1039366/93‑1 (MQ=255)
                                                                         gACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTAAAGCTGGAGCAGTg                 >  1:508040/1‑93 (MQ=255)
                                                                         gACAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTAAAGCTGGAGCAGTg                 >  1:4653408/1‑93 (MQ=255)
                                                                          aCAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTAAAGCTGGAGCAGTGc                >  1:1535849/1‑93 (MQ=255)
                                                                          aCAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTAAAGCTGGAGCAGTGc                >  2:3262985/1‑93 (MQ=255)
                                                                           cAGCCCATTCACCCAGACGGCGGGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTAAAGCTGGAGCAGTGcc               >  1:5999999/1‑93 (MQ=255)
                                                                           cAGCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTAAAGCTGGAGCAGTGcc               >  2:6720603/1‑93 (MQ=255)
                                                                             gCCCATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTAAAGCTGGAGCAGTGCcac             >  2:5938532/1‑93 (MQ=255)
                                                                              cccATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTAAAGCTGGAGCAGTGCcaca            >  1:4528000/1‑93 (MQ=255)
                                                                              cccATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTAAAGCTGGAGCAGTGCcaca            <  1:756236/93‑1 (MQ=255)
                                                                               ccATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTAAAGCTGGAGCAGTGCCACAg           <  2:3054643/93‑1 (MQ=255)
                                                                                cATTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTAAAGCTGGAGCAGTGCCACAGc          <  2:2217076/93‑1 (MQ=255)
                                                                                 aTTCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTAAAGCTGGAGCAGTGCCACAGCg         >  1:2280713/1‑93 (MQ=255)
                                                                                   tCACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTAAAGCTGGAGCAGTGCCACAGCGTa       >  1:4720313/1‑93 (MQ=255)
                                                                                    cACCCAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTAAAGCTGGAGCAGTGCCACAGCGTAt      >  2:5543549/1‑93 (MQ=255)
                                                                                       ccAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTAAAGCTGGGGCAGTGCCACAGCGTATCaa   <  2:4643332/93‑1 (MQ=255)
                                                                                        cAGGCGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTAAAGCTGGAGCAGTGCCACAGCGTATCAAt  >  2:2514210/1‑93 (MQ=255)
                                                                                        cAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTAAAGCTGGAGCAGTGCCACAGCGTATCAAt  <  1:3926933/93‑1 (MQ=255)
                                                                                        cAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTAAAGCTGGAGCAGTGCCACAGCGTATCAAt  <  1:5378168/93‑1 (MQ=255)
                                                                                        cAGACGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTAAAGCTGGAGCAGTGCCACAGCGTATCAAt  <  2:4574753/93‑1 (MQ=255)
                                                                                           |                                                                                         
TATCGCGATGGTGGTGATAATCGAAGTACAGGGAGGGGGAAGAGATATGAATAAAGTTACGCACACCCCAGGCGACAGCCCATTCACCCAGGCGGCGAGTGGCGCGAACGTTAGCCAGATCGAAAGCCTGTTGTGTCCCCCAGGGTGAGGTAAAGCTGGAGCAGTGCCACAGCGTATCAAT  >  NC_000913/905598‑905778

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