Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,003,385 Δ1 bp coding (186/681 nt) ycbF → putative chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,003,3840A.88.8% 707.3 / 75.4 206coding (185/681 nt)ycbFputative chaperone
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (119/64);  minor base G (12/11);  total (131/75)
Fisher's exact test for biased strand distribution p-value = 2.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.01e-02

Only 100 of 235 total aligned reads displayed.

AGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAA  >  NC_000913/1003292‑1003473
                                                                                             |                                                                                         
 gaATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGag                                                                                          >  2:4775400/1‑92 (MQ=255)
   aTCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAGc                                                                                        >  1:3481190/1‑93 (MQ=255)
      gATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCCACCGCCCCACGCGTTGGGCCGGGGGGGAGCGTca                                                                                    >  1:1265852/1‑92 (MQ=255)
      gATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCg                                                                                    >  2:6532288/1‑93 (MQ=255)
       atatCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGa                                                                                   <  1:1301201/93‑1 (MQ=255)
       atatCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGa                                                                                   >  2:9195902/1‑93 (MQ=255)
        tatCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGAc                                                                                  >  1:5865793/1‑93 (MQ=255)
         atCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACa                                                                                 >  2:6591525/1‑93 (MQ=255)
         atCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACa                                                                                 >  1:6064901/1‑93 (MQ=255)
         atCTATTTAAGGGACTAGTACGAGTCAGGAAGCGTGGCTGATTAGCAGTAAATACACCGTCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACa                                                                                 <  2:5013940/93‑1 (MQ=255)
          tCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGGCGTCGACAg                                                                                >  2:3749487/1‑93 (MQ=255)
           cTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGt                                                                               >  2:9648603/1‑93 (MQ=255)
            tATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGGAGCGTCGACAGt                                                                               >  1:9840881/1‑93 (MQ=255)
            tATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTg                                                                              >  1:579434/1‑93 (MQ=255)
            tATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTg                                                                              <  1:5237066/93‑1 (MQ=255)
            tATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTg                                                                              >  1:2650358/1‑93 (MQ=255)
            tATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTg                                                                              >  2:2187289/1‑93 (MQ=255)
            tATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTg                                                                              >  1:4201741/1‑93 (MQ=255)
             tttttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGc                                                                             <  1:4828423/92‑1 (MQ=255)
              ttttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGcc                                                                            >  1:4340453/1‑93 (MQ=255)
                ttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAg                                                                          >  2:9626112/1‑93 (MQ=255)
                ttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAg                                                                          >  1:3662093/1‑93 (MQ=255)
                  tcTGACCAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCa                                                                        <  1:1065218/92‑1 (MQ=255)
                  aCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCa                                                                        >  2:3358436/1‑93 (MQ=255)
                     gACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCa                                                                     <  1:4974177/93‑1 (MQ=255)
                     gACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCa                                                                     <  2:3423119/93‑1 (MQ=255)
                     gACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCa                                                                     >  1:2238474/1‑93 (MQ=255)
                        tAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTa                                                                  >  1:5964082/1‑93 (MQ=255)
                         aaTACCAGTCAGGAATCCTGTCTTATAAACAGTAAAATCACCCGCCCAACGGGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTAc                                                                 <  2:8452388/93‑1 (MQ=255)
                          aTACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACt                                                                >  2:9663329/1‑93 (MQ=255)
                           tACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTg                                                               >  1:9575780/1‑93 (MQ=255)
                           aaCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTg                                                               <  2:6417741/92‑1 (MQ=255)
                             ccAGTGAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGc                                                             <  2:7823055/93‑1 (MQ=255)
                               aGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGccc                                                           >  1:3083211/1‑92 (MQ=255)
                                gTCAGGAATCCTGGCTTATAAACAGTAAAATCAACCGCCCAACGGGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGACGc                                                          <  1:6254923/93‑1 (MQ=255)
                                 tCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCt                                                         >  1:7348410/1‑93 (MQ=255)
                                 tCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCt                                                         >  1:1919354/1‑93 (MQ=255)
                                   agaaaTCCGGGTTTATAAAAATAAAAACCACCCGCCCAACGCTTGGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc                                                       <  2:9476113/90‑1 (MQ=255)
                                   aGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc                                                       <  1:8470737/93‑1 (MQ=255)
                                    ggAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCg                                                      >  1:6488188/1‑93 (MQ=255)
                                     gAATCCTGGCTTATTGACACTAAAATCAACCGCCCAACCCGATGGGCCGGGGGGGCGCCTCCGCCGAGCCAGCACCCTTACCGGGCGCCCCGc                                                     >  2:2111335/1‑93 (MQ=255)
                                       aTCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGG‑‑AGCCTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCAc                                                  >  2:4488878/1‑93 (MQ=255)
                                        tCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCAc                                                  <  2:6879093/93‑1 (MQ=255)
                                        tCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCAc                                                  >  2:2784858/1‑93 (MQ=255)
                                         ccTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACt                                                 >  1:475187/1‑93 (MQ=255)
                                         ccTGGCTTATTAACAGAAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACt                                                 <  1:8372651/93‑1 (MQ=255)
                                          cTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACtt                                                >  1:3291910/1‑93 (MQ=255)
                                            ggCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTAt                                              >  2:9099418/1‑93 (MQ=255)
                                            ggCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTAt                                              >  2:4777289/1‑93 (MQ=255)
                                            ggCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTAt                                              >  1:4929602/1‑93 (MQ=255)
                                               ttattaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATtct                                           >  1:9199603/1‑93 (MQ=255)
                                                tattaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATtctc                                          >  2:9610668/1‑93 (MQ=255)
                                                 attaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTcc                                         >  1:8796054/1‑93 (MQ=255)
                                                  ttaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCt                                        >  2:5245894/1‑93 (MQ=255)
                                                  ttaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCt                                        >  2:3980093/1‑93 (MQ=255)
                                                  ttaACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCt                                        >  2:6319752/1‑93 (MQ=255)
                                                    aaCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGa                                      >  1:9915344/1‑93 (MQ=255)
                                                    aaCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGGGCCAGCACCATACCTGGCCGCTCCCCCCCTTATtctcttta                                      >  2:6092860/1‑89 (MQ=255)
                                                       aGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGc                                   >  1:2506038/1‑93 (MQ=255)
                                                         aaaaaTCAACCGCCCAACGGGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGccc                                 <  2:3080475/92‑1 (MQ=255)
                                                          aaaaTCACCCGCCCAAGGCGTTGGGCCGGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGccc                                 <  1:5336383/93‑1 (MQ=255)
                                                          aaaaTCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCg                                <  1:5004741/93‑1 (MQ=255)
                                                            aaTCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGt                              >  2:8065891/1‑93 (MQ=255)
                                                             accAACCCCCCAGCGGTTGGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTa                             <  1:2337087/91‑1 (MQ=255)
                                                             aTCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTa                             >  2:6474207/1‑93 (MQ=255)
                                                             aTCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTa                             >  2:6869085/1‑93 (MQ=255)
                                                             aTCAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTa                             >  2:8307329/1‑93 (MQ=255)
                                                             aTCAACCCCCCAACGGCGTGGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGt                              <  2:4337854/93‑1 (MQ=255)
                                                              acAACCGCCCACCGCGTTGGGCCGGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTa                             <  1:1467968/92‑1 (MQ=255)
                                                               cAACCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACg                           >  1:8123036/1‑93 (MQ=255)
                                                                aaCCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGa                          >  1:6546445/1‑93 (MQ=255)
                                                                aaCCGCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGCACGa                          >  1:3340188/1‑93 (MQ=255)
                                                                 aCCGCCCACCGCGTGGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGAc                         <  2:1447215/93‑1 (MQ=255)
                                                                  ccccccACCGCGTTGGCCCGGGGGGG‑‑GCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAg                       <  2:2866999/90‑1 (MQ=255)
                                                                  ccGCCCACCGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACa                        <  1:2658072/93‑1 (MQ=255)
                                                                   cGCCTAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAg                       <  1:1668510/93‑1 (MQ=255)
                                                                    gCCCAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAgg                      >  2:3579502/1‑93 (MQ=255)
                                                                     cccAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGc                     >  2:6101856/1‑93 (MQ=255)
                                                                     cacAACGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGATCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGc                     <  2:5277240/91‑1 (MQ=255)
                                                                      ccACCGGGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCa                    <  2:8422991/93‑1 (MQ=255)
                                                                         aCGCGTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGt                 >  1:9004479/1‑93 (MQ=255)
                                                                            gcgTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTg               >  2:4598438/1‑93 (MQ=255)
                                                                            gcgTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTg               <  1:3588448/93‑1 (MQ=255)
                                                                              gTTGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTgcg             >  2:1751079/1‑93 (MQ=255)
                                                                              gTTGGGCCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTgc              <  2:2475726/93‑1 (MQ=255)
                                                                               ttGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTgcgc            <  2:10014209/93‑1 (MQ=255)
                                                                                tGGGCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCt           <  2:4988894/93‑1 (MQ=255)
                                                                                ggggCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCt           <  2:3021462/92‑1 (MQ=255)
                                                                                 gggCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCtt          >  1:7023998/1‑93 (MQ=255)
                                                                                  ggCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTg         <  1:4887026/93‑1 (MQ=255)
                                                                                   gCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGc        <  2:1329072/93‑1 (MQ=255)
                                                                                   gCCGGGGGGGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGc        >  2:3163394/1‑93 (MQ=255)
                                                                                     gcggggggGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCt       <  1:2834143/91‑1 (MQ=255)
                                                                                       ggggggGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTga     <  2:492155/93‑1 (MQ=255)
                                                                                       ggggggGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTga     >  2:7280146/1‑93 (MQ=255)
                                                                                       ggggggGAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTga     <  2:2039221/93‑1 (MQ=255)
                                                                                       cggggggggCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTga     <  1:9710324/92‑1 (MQ=255)
                                                                                        ggggggAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTgag    >  2:710450/1‑93 (MQ=255)
                                                                                          ggggAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGaa  >  1:945328/1‑93 (MQ=255)
                                                                                          ggggAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGaa  <  2:9773666/93‑1 (MQ=255)
                                                                                             |                                                                                         
AGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAA  >  NC_000913/1003292‑1003473

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