Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,374,263 T→C Q50R (CAA→CGA)  sspB ← stringent starvation protein B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,374,2630TC99.5% 686.7 / ‑4.3 215Q50R (CAA→CGA) sspBstringent starvation protein B
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (117/97);  minor base A (0/1);  total (117/98)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01

Only 100 of 222 total aligned reads displayed.

GGAATGCCACCAAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTTGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCAGCAAC  >  NC_000913/3374173‑3374352
                                                                                          |                                                                                         
 gAATGCCACCAAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGcc                                                                                        >  2:623152/1‑93 (MQ=255)
 gAATGCCACCAAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGcc                                                                                        >  1:1002176/1‑93 (MQ=255)
    tGCCACCAAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGGACGCGGCGCAATGTTGAGTACGATTCGACCgc                                                                                     >  2:662829/1‑92 (MQ=255)
      ccaccaAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCa                                                                                   <  1:3724013/93‑1 (MQ=255)
      ccaccaAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCa                                                                                   <  1:4374389/93‑1 (MQ=255)
       caccaAAGCGCGCGTTAAAGCGCACCTCATCATTCGCTAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCAc                                                                                  <  1:3965118/93‑1 (MQ=255)
       caccaAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCAc                                                                                  <  2:680799/93‑1 (MQ=255)
       caccaAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCAc                                                                                  >  1:2305325/1‑93 (MQ=255)
        cccaAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCAcg                                                                                 <  1:607863/92‑1 (MQ=255)
         ccaAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCAcgc                                                                                >  1:1070135/1‑93 (MQ=255)
         ccaAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCAcgc                                                                                <  1:2689037/93‑1 (MQ=255)
         ccaAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCAcgc                                                                                <  2:1406036/93‑1 (MQ=255)
           aaaGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACgcgc                                                                              >  1:4413583/1‑93 (MQ=255)
            aaGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCa                                                                             <  2:417558/93‑1 (MQ=255)
            aaGCGCGCGTTAAAGCGCACATCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCa                                                                             <  2:2600147/93‑1 (MQ=255)
              gcgcgcGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCata                                                                           >  2:2957229/1‑93 (MQ=255)
               cgcgcgTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCatat                                                                          <  2:1369845/93‑1 (MQ=255)
               cgcgcgTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCatat                                                                          <  2:4213660/93‑1 (MQ=255)
               agcgcgTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGCCACGCGCatat                                                                          >  2:4595038/2‑93 (MQ=255)
                 cgcgTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGGCAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTc                                                                        <  1:2636895/93‑1 (MQ=255)
                 cgcgTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATGGCCGACAGCACGCGGCGCAATGTTGAGTACGATGCGCCCGTCACGCGCATATTc                                                                        <  2:2783981/93‑1 (MQ=255)
                    gTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCAt                                                                     >  1:672627/1‑93 (MQ=255)
                    gTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCAt                                                                     >  2:360754/1‑93 (MQ=255)
                    gTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCAt                                                                     >  1:796552/1‑93 (MQ=255)
                    gTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCAt                                                                     >  1:3237128/1‑93 (MQ=255)
                     ttAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATa                                                                    >  2:1699962/1‑93 (MQ=255)
                      tAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAg                                                                   <  1:1338387/93‑1 (MQ=255)
                      tAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAg                                                                   <  2:1997003/93‑1 (MQ=255)
                       aaaGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAgg                                                                  >  2:2211034/1‑93 (MQ=255)
                       aaaGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCAAGCGCATATTCCATacg                                                                  >  2:4351464/1‑91 (MQ=255)
                        aaGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGa                                                                 <  2:882427/93‑1 (MQ=255)
                          gcgcACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAAc                                                               >  1:2121206/1‑93 (MQ=255)
                            gcACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCt                                                             >  1:2310316/1‑93 (MQ=255)
                            gcACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCt                                                             >  1:3873097/1‑93 (MQ=255)
                             cACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATATGAACCTg                                                            >  1:1728801/1‑93 (MQ=255)
                              aCCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGc                                                           <  1:2476119/93‑1 (MQ=255)
                              aCCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGc                                                           >  2:2431868/1‑93 (MQ=255)
                               ccTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGTAACCTGCa                                                          >  2:2664168/1‑93 (MQ=255)
                               ccTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCa                                                          >  1:543822/1‑93 (MQ=255)
                                cTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCAc                                                         >  1:1544463/1‑93 (MQ=255)
                                cTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCAc                                                         >  1:2888783/1‑93 (MQ=255)
                                cTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCAc                                                         >  1:1193829/1‑93 (MQ=255)
                                  catcatTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGc                                                       <  1:1538994/93‑1 (MQ=255)
                                   atcatTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGcc                                                      <  1:1869127/93‑1 (MQ=255)
                                    tcatTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCa                                                     <  1:2727804/93‑1 (MQ=255)
                                    tcatTCGCCAGTTCCACATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCa                                                     <  1:2983846/93‑1 (MQ=255)
                                      atTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAgg                                                   >  1:1649158/1‑93 (MQ=255)
                                      atTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAgg                                                   >  2:912078/1‑93 (MQ=255)
                                       ttCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAggg                                                  >  2:1674017/1‑93 (MQ=255)
                                        tCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGa                                                 <  1:499229/93‑1 (MQ=255)
                                         cGCCAGTTCCAGATTGCCGACAGCGCGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAg                                                <  2:3090051/93‑1 (MQ=255)
                                          gCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGc                                               <  1:2426571/93‑1 (MQ=255)
                                          gCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGc                                               >  1:4050038/1‑93 (MQ=255)
                                          cccAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGc                                               <  2:3920608/92‑1 (MQ=255)
                                           ccAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCg                                              >  2:2711414/1‑93 (MQ=255)
                                            cAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGt                                             >  2:1424177/1‑93 (MQ=255)
                                              gTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTca                                           >  1:4326088/1‑93 (MQ=255)
                                               ttCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTcac                                          >  2:3712062/1‑93 (MQ=255)
                                               ttCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTcac                                          >  2:2903069/1‑93 (MQ=255)
                                               ttCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTcac                                          <  1:3634221/93‑1 (MQ=255)
                                               ttCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTcac                                          >  2:4051570/1‑93 (MQ=255)
                                                tCCAGATTGCCGACAGCACTCGGCGCAATGTTGAGAGAGATAGGCCTGTCGCGCGCAGATGCCAGAGAAACCAGCGCGCCAGTGAGCGTcaca                                         <  2:3963503/93‑1 (MQ=255)
                                                tCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTcaca                                         >  1:3508363/1‑93 (MQ=255)
                                                    gATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATcca                                     >  2:3463404/1‑93 (MQ=255)
                                                     aTTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATccac                                    <  1:2537726/93‑1 (MQ=255)
                                                      ttGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATccacc                                   <  2:4028934/93‑1 (MQ=255)
                                                      ttGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATccacc                                   >  2:718087/1‑93 (MQ=255)
                                                        gCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATccaccac                                 <  1:4587694/93‑1 (MQ=255)
                                                          cGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATccaccacca                               >  1:4093394/1‑93 (MQ=255)
                                                            aCAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAgg                             >  2:4391195/1‑93 (MQ=255)
                                                              aGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTg                           >  1:4616475/1‑93 (MQ=255)
                                                              aGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTg                           >  2:1290551/1‑93 (MQ=255)
                                                                 aCGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTgcgg                        >  1:4448701/1‑93 (MQ=255)
                                                                  cgcgGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTgcggc                       >  1:851569/1‑93 (MQ=255)
                                                                   gcggcgCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTgcggcg                      <  2:650722/93‑1 (MQ=255)
                                                                    cggcgCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGt                     >  2:1041251/1‑93 (MQ=255)
                                                                     ggcgCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTg                    <  1:68827/93‑1 (MQ=255)
                                                                     ggcgCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTg                    >  1:462844/1‑93 (MQ=255)
                                                                      gcgcAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGa                   >  1:3443800/1‑93 (MQ=255)
                                                                        gcAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGc                 <  1:3360120/93‑1 (MQ=255)
                                                                        gcAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGc                 >  1:1924843/1‑93 (MQ=255)
                                                                         cAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCt                <  1:118667/93‑1 (MQ=255)
                                                                          aaTGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTg               >  2:2127920/1‑93 (MQ=255)
                                                                           aTGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTgg              >  1:2248494/1‑93 (MQ=255)
                                                                            tGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGt             <  2:547180/93‑1 (MQ=255)
                                                                            tGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGt             >  2:1800196/1‑93 (MQ=255)
                                                                             gTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGtt            >  2:4104820/1‑93 (MQ=255)
                                                                              ttGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTa           <  2:2087389/93‑1 (MQ=255)
                                                                                 aGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATcc        >  1:2128174/1‑93 (MQ=255)
                                                                                  gTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCCGGGATCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATcct       >  2:3148408/1‑92 (MQ=255)
                                                                                  gTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCa       <  2:2326900/93‑1 (MQ=255)
                                                                                  gTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCa       <  1:1507613/93‑1 (MQ=255)
                                                                                  gTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCa       <  2:1176512/93‑1 (MQ=255)
                                                                                  gTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCa       >  2:3295409/1‑93 (MQ=255)
                                                                                   tACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCAg      >  2:3206972/1‑93 (MQ=255)
                                                                                   tACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCAg      <  1:340956/93‑1 (MQ=255)
                                                                                      gATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCAGCaa   >  2:3409151/1‑93 (MQ=255)
                                                                                       aTTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCAGCAAc  <  2:2241314/93‑1 (MQ=255)
                                                                                       aTTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCAGCAAc  >  1:4137052/1‑93 (MQ=255)
                                                                                       aTTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCAGCAAc  <  1:3065651/93‑1 (MQ=255)
                                                                                          |                                                                                         
GGAATGCCACCAAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTTGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCAGCAAC  >  NC_000913/3374173‑3374352

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: