Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,521,565 (A)6→7 coding (456/1350 nt) yjhF ← KpLE2 phage‑like element; putative transporter YjhF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,521,5591.A99.4% 1462.8 / ‑0.0 343coding (462/1350 nt)yjhFKpLE2 phage‑like element; putative transporter YjhF
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base A (157/184);  total (157/186)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 390 total aligned reads displayed.

GAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGC‑AAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGATATTTTTGCTTCACGGGCAATAGCGAATAAAATAGG  >  NC_000913/4521459‑4521654
                                                                                                     |                                                                                               
gAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCggcggc                                                                                                  >  2:8784905/1‑101 (MQ=255)
  aTAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGC‑aa                                                                                               <  2:8784254/101‑1 (MQ=255)
  aTAAAGGTTCGAATGCCGACAATGATACCGTATATCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGC‑aa                                                                                               >  1:8784908/1‑101 (MQ=255)
   tAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGC‑aaa                                                                                              >  2:8784913/1‑101 (MQ=255)
    aaaGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGC‑aaaa                                                                                             <  2:8784312/101‑1 (MQ=255)
     aaGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGC‑aaaaa                                                                                            >  1:8784915/1‑101 (MQ=255)
      aGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGC‑aaaaaa                                                                                           <  1:8784624/101‑1 (MQ=255)
       ggTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGC‑aaaaaaa                                                                                          <  2:8784064/101‑2 (MQ=255)
        gTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGC‑aaaaaaac                                                                                         >  1:8784920/1‑99 (MQ=255)
          tGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGC‑a                                                                                                >  1:8784926/1‑92 (MQ=255)
           ggAATGCCGACAATGATACCGTATATCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGC‑aaaaaaaccgt                                                                                      <  2:8784437/101‑6 (MQ=255)
            gAATGCCGACAATGATACCGTATATCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCgtg                                                                                      <  1:8783947/101‑1 (MQ=255)
             aaTGCCGACAATGATACCGTATATCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCgtgt                                                                                     >  1:8784929/1‑101 (MQ=255)
             aaTGCCGACAATGATACCGTATATCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCgtgt                                                                                     >  2:8784930/1‑101 (MQ=255)
               tGCCGACAATGATACCGTATATCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGc                                                                                   <  1:8784370/101‑1 (MQ=255)
               tGCCGACAATGATACCGTATATCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGc                                                                                   <  1:8784249/101‑1 (MQ=255)
                 ccGACAATGATACCGTATATCAGTACCAATCCTACATCAGCGCTATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGa                                                                                 >  1:8784939/1‑101 (MQ=255)
                   gACAATGATACCGTATATCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGAcc                                                                               <  2:8784335/101‑1 (MQ=255)
                    aCAATGATACCGTATATCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCa                                                                              <  2:8784402/101‑1 (MQ=255)
                    aCAATGATACCGTATATCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCa                                                                              >  1:8784945/1‑101 (MQ=255)
                    aCAATGATACCGTATATCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCa                                                                              >  2:8784947/1‑101 (MQ=255)
                    aCAATGATACCGTATATCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCa                                                                              <  1:8784108/101‑1 (MQ=255)
                      aaTGATACCGTATATCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGc                                                                            >  1:8784953/1‑101 (MQ=255)
                      aaTGATACCGTATATCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGc                                                                            >  2:8784954/1‑101 (MQ=255)
                        tGATACCGTATATCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCaa                                                                          <  2:8783881/101‑1 (MQ=255)
                           tACCGTATATCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAAcc                                                                       >  2:8784962/1‑101 (MQ=255)
                            aCCGTATATCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCa                                                                      <  2:8784537/101‑1 (MQ=255)
                             ccGTATATCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCggcggc                                                                                                  <  2:8784532/72‑1 (MQ=255)
                             ccGTATATCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAg                                                                     <  2:8784224/101‑1 (MQ=255)
                             ccGTATATCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAg                                                                     <  2:8784514/101‑1 (MQ=255)
                                tataTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGaaa                                                                  <  1:8783700/101‑1 (MQ=255)
                                tataTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGaaa                                                                  >  2:8784966/1‑101 (MQ=255)
                                tataTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGaaa                                                                  <  1:8784174/101‑1 (MQ=255)
                                tataTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGaaa                                                                  >  1:8784967/1‑101 (MQ=255)
                                 atatCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAg                                                                 <  2:8784138/101‑1 (MQ=255)
                                 atatCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAg                                                                 <  2:8784482/101‑1 (MQ=255)
                                  tatCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGc                                                                <  2:8784347/101‑1 (MQ=255)
                                   atCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCa                                                               >  1:8784975/1‑101 (MQ=255)
                                      aGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCg                                                            >  1:8784981/1‑101 (MQ=255)
                                      aGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCg                                                            >  2:8784984/1‑101 (MQ=255)
                                        tACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAAc                                                                        >  2:8784989/1‑87 (MQ=255)
                                         aCCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGca                                                         >  1:8784990/1‑101 (MQ=255)
                                           cAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGcaca                                                       <  2:8784365/101‑1 (MQ=255)
                                            aaTCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGcacac                                                      >  2:8784997/1‑101 (MQ=255)
                                               ccTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCg                                                                                  >  2:8785004/1‑70 (MQ=255)
                                               ccTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGcacacaca                                                   >  1:8785003/1‑101 (MQ=255)
                                                cTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACAt                                                  >  2:8785009/1‑101 (MQ=255)
                                                cTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACAt                                                  <  2:8784937/101‑1 (MQ=255)
                                                  aCATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCaa                                                                          <  1:8784444/75‑1 (MQ=255)
                                                  aCATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATaa                                                >  2:8785016/1‑101 (MQ=255)
                                                  aCATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATaa                                                >  1:8785017/1‑101 (MQ=255)
                                                   cATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATaaa                                               <  2:8784630/101‑1 (MQ=255)
                                                   cATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATaaa                                               >  1:8785021/1‑101 (MQ=255)
                                                     tCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACa                                             <  1:8784461/101‑1 (MQ=255)
                                                       aGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATc                                           <  2:8784451/101‑1 (MQ=255)
                                                       aGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATc                                           <  1:8784319/101‑1 (MQ=255)
                                                       aGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATc                                           >  2:8785025/1‑101 (MQ=255)
                                                          gcCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGt                                        >  1:8785028/1‑101 (MQ=255)
                                                           ccATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTg                                       <  1:8784704/101‑1 (MQ=255)
                                                           ccATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTg                                       <  1:8784613/101‑1 (MQ=255)
                                                           ccATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTg                                       >  2:8785032/1‑101 (MQ=255)
                                                             atatTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGat                                     >  1:8785036/1‑101 (MQ=255)
                                                             atatTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGat                                     <  1:8784507/101‑1 (MQ=255)
                                                             atatTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGat                                     <  1:8784396/101‑1 (MQ=255)
                                                              tatTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGata                                    <  2:8784563/101‑1 (MQ=255)
                                                                ttCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGATAtt                                  >  2:8785041/1‑101 (MQ=255)
                                                                 tCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGaa                                                                   <  1:8784490/67‑1 (MQ=255)
                                                                  cccTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCa                                                               <  2:8784279/70‑1 (MQ=255)
                                                                  cccTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGATAtttt                                >  2:8785045/1‑101 (MQ=255)
                                                                   ccTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAg                                                                     <  2:8784105/63‑1 (MQ=255)
                                                                   ccTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCGTCGGCACACACATAAACATCGGTGATAttttt                               <  2:8784302/101‑1 (MQ=255)
                                                                   ccTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGATAttttt                               >  1:8785050/1‑101 (MQ=255)
                                                                   ccTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGATAtt                                  <  1:8784217/98‑1 (MQ=255)
                                                                    cTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGATATTTTTg                              <  1:8784385/101‑1 (MQ=255)
                                                                    cTGGCGATAACAGTCGGGCCGGGGGGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGATATTTTTg                              >  2:8785053/1‑101 (MQ=255)
                                                                      ggCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGATATTTTTGCt                            <  1:8784028/101‑1 (MQ=255)
                                                                       gCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGATATTTTTGCtt                           >  2:8785058/1‑101 (MQ=255)
                                                                        cGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGATATTTTTGCTTc                          <  2:8784189/101‑1 (MQ=255)
                                                                        cGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGATATTTTTGCTTc                          <  1:8784427/101‑1 (MQ=255)
                                                                        cGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGATATTTTTGCTTc                          >  1:8785060/1‑101 (MQ=255)
                                                                         gATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGATATTTTTGCTTCa                         <  1:8784456/101‑1 (MQ=255)
                                                                          aTAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGATATTTTTGCTTCAc                        >  2:8785064/1‑101 (MQ=255)
                                                                           tAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGAc                                                                                >  1:8785070/1‑44 (MQ=38)
                                                                             aCAGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGATATTTTTGCTTCACggg                     <  1:8784263/101‑1 (MQ=255)
                                                                               aGTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGATATTTTTGCTTCACGGGCa                   >  1:8785075/1‑101 (MQ=255)
                                                                                gTCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGATATTTTTGCTTCACGGGCaa                  >  2:8785077/1‑101 (MQ=255)
                                                                                 tCGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGATATTTTTGCTTCACGGGCAAt                 <  1:8784813/101‑1 (MQ=255)
                                                                                  cGGGCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGATATTTTTGCTTCACGGGCAATa                <  1:8784554/101‑1 (MQ=255)
                                                                                    ggCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGATATTTTTGCTTCACGGGCAATAGc              <  2:8784392/101‑1 (MQ=255)
                                                                                    ggCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGATATTTTTGCTTCACGGGCAATAGc              <  1:8784545/101‑1 (MQ=255)
                                                                                     gCCGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGATATTTTTGCTTCACGGGCAATAGCg             <  2:8784418/101‑1 (MQ=255)
                                                                                      ccGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGata                                    <  1:8784310/77‑1 (MQ=255)
                                                                                      ccGGGGTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGATATTTTTGCTTCACGGGCAATAGCGa            >  1:8785082/1‑101 (MQ=255)
                                                                                        ggggTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGATATTTTTGCTTCACGGGCAATAGCGAAt          <  2:8784659/101‑1 (MQ=255)
                                                                                        ggggTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGATATTTTTGCTTCACGGGCAATAGCGAAt          >  2:8785087/1‑101 (MQ=255)
                                                                                         gggTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGATATTTTTGCTTCACGGGCAATAGCGAATa         <  1:8784500/101‑1 (MQ=255)
                                                                                           gTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGATATTTTTGCTTCACGGGCAATAGCGAATaaa       <  2:8784712/101‑1 (MQ=255)
                                                                                           gTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGATATTTTTGCTTCACGGGCAATAGCGAATaaa       >  1:8785089/1‑101 (MQ=255)
                                                                                           gTGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGATATTTTTGCTTCACGGGCAATAGCGAATa         <  1:8784657/99‑1 (MQ=255)
                                                                                            tGCGGCGGCAAAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGATATTTTTGCTTCACGGGCAATAGCGAATaaaa      >  2:8785093/1‑101 (MQ=255)
                                                                                                     |                                                                                               
GAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGC‑AAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGATATTTTTGCTTCACGGGCAATAGCGAATAAAATAGG  >  NC_000913/4521459‑4521654

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