Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,648,144 (T)7→6 100% intergenic (‑356/‑384) dinQ ← / → arsR UV‑inducible membrane toxin, DinQ‑AgrB type I toxin‑antitoxin system/arsenical resistance operon transcriptional repressor; autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,648,1380T.95.4% 686.0 / 21.7 172intergenic (‑350/‑390)dinQ/arsRUV‑inducible membrane toxin, DinQ‑AgrB type I toxin‑antitoxin system/arsenical resistance operon transcriptional repressor; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (1/6);  new base . (88/76);  total (89/83)
Fisher's exact test for biased strand distribution p-value = 5.59e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 192 total aligned reads displayed.

AGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTC  >  NC_000913/3648040‑3648228
                                                                                                  |                                                                                          
aGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtt                                                                                           <  2:4042948/100‑1 (MQ=255)
  tATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttt                                                                                         <  2:2741510/100‑1 (MQ=255)
  tATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttt                                                                                         <  2:2295930/100‑1 (MQ=255)
    tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                       <  2:3841054/100‑1 (MQ=255)
    aggTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                       >  1:3250983/2‑100 (MQ=255)
     ggTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttta                                                                                      <  2:3475583/100‑2 (MQ=255)
      gTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttttt                                                                                      >  1:1857536/1‑99 (MQ=255)
        aTTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttttacat                                                                                   >  1:240965/1‑96 (MQ=255)
          tAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATcc                                                                                <  1:2710550/100‑1 (MQ=255)
           aGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCg                                                                               >  2:460273/1‑100 (MQ=255)
           aGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCg                                                                               >  2:4192027/1‑100 (MQ=255)
           aGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGATGGAAAGTAACAAACCGACGCAAAGTCGCTTTTTTACATCCg                                                                               <  2:2978943/100‑1 (MQ=21)
            gCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCg                                                                               >  1:341460/1‑99 (MQ=255)
               ttGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGAtt                                                                           <  2:1386965/100‑1 (MQ=255)
               ttGATGCGTTAGATGATCACTATCTTAGTCCAGCCAGAACCACCAACTCATGGTTGTAAGGTAAAAATCCGACACAAAGTGGGTTTTTTTAcg                                                                                   <  1:3517921/93‑2 (MQ=255)
                  aTGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCgg                                                                        <  1:3156551/100‑1 (MQ=255)
                  aTGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCgg                                                                        <  2:301722/100‑1 (MQ=255)
                    gCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGAc                                                                      <  2:2197700/100‑1 (MQ=255)
                    gCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGAc                                                                      <  2:1683067/100‑1 (MQ=255)
                     cGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACCTCCGGGTTCGGa                                                                       >  1:3908841/1‑98 (MQ=18)
                     cGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACa                                                                     >  1:1174518/1‑100 (MQ=255)
                     cGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACa                                                                     >  2:3172217/1‑100 (MQ=255)
                     cGTGAGATGATGGCTATCTCACTCCACTCAGAGCCACCAACTCAGTGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGAc                                                                      <  1:4074767/99‑1 (MQ=18)
                       ttAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAg                                                                   <  1:3664313/100‑1 (MQ=255)
                        tAGATGTTGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAgg                                                                  <  1:2169659/100‑1 (MQ=255)
                         aGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGc                                                                 <  2:1781762/100‑1 (MQ=255)
                         aGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGc                                                                 <  1:3765111/100‑1 (MQ=255)
                          gatgatGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCt                                                                >  2:1576839/1‑100 (MQ=255)
                          gatgatGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCt                                                                <  1:1334406/100‑1 (MQ=255)
                             gatgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTAGGTTTTTTACATCCGGATTCGGAAAAAGCTTa                                                              >  2:750321/1‑99 (MQ=18)
                               tgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAata                                                           <  2:374154/100‑1 (MQ=255)
                               tgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAata                                                           <  1:2371455/100‑1 (MQ=255)
                               tgGCAATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAATACCGACGCAAACTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAat                                                            <  2:640905/99‑1 (MQ=18)
                                 gCTATCGCACGGCAGTCAGAGCCACCAACTCAGGGTTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAatat                                                          <  1:926787/99‑1 (MQ=18)
                                  cTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGa                                                        <  1:753393/100‑1 (MQ=255)
                                    aTCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAANTCGNTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACg                                                      <  2:565764/100‑1 (MQ=255)
                                     tctcACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGGTTCGGACAACGCTTAATATGACg                                                      >  2:812646/1‑99 (MQ=18)
                                     tctcACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGa                                                     >  2:3539036/1‑100 (MQ=255)
                                     tctcACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGa                                                     <  2:2828783/100‑1 (MQ=255)
                                     tctcACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGa                                                     <  1:3627239/100‑1 (MQ=255)
                                      ctcACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGAt                                                    >  2:3413602/1‑100 (MQ=255)
                                       tcACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATg                                                   >  2:2371105/1‑100 (MQ=255)
                                        cACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATg                                                   >  2:2534180/1‑99 (MQ=255)
                                          cTCCAGTCAGAGCCACCAACTCAGTGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAcc                                                <  1:1397503/100‑1 (MQ=255)
                                          cTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAcc                                                >  1:1496870/1‑100 (MQ=255)
                                          cTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAcc                                                >  1:3196716/1‑100 (MQ=255)
                                          cTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAcc                                                >  2:152126/1‑100 (MQ=255)
                                            ccAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGGCCCa                                              >  2:3334123/1‑100 (MQ=255)
                                            ccAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCa                                              >  2:1624515/1‑100 (MQ=255)
                                             cAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAg                                             >  2:4275108/1‑100 (MQ=255)
                                             cAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAg                                             <  1:844978/100‑1 (MQ=255)
                                              aGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGt                                            >  2:3959656/1‑100 (MQ=255)
                                              aGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGt                                            >  2:270708/1‑100 (MQ=255)
                                               gTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGGCCCAGTg                                           <  2:2309010/100‑1 (MQ=255)
                                                 cAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATTTGACGATGACCCAGTGaa                                         >  2:3375763/1‑100 (MQ=255)
                                                 cAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGAAGATGACCCAGTGaa                                         >  2:2252282/1‑100 (MQ=255)
                                                   gagCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAg                                       <  1:375933/100‑1 (MQ=255)
                                                     gCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGta                                     <  1:1963665/100‑1 (MQ=255)
                                                       caccaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtata                                   <  1:4291187/100‑1 (MQ=255)
                                                       caccaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtata                                   <  2:500971/100‑1 (MQ=255)
                                                          caACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATaa                                >  2:3296538/1‑100 (MQ=255)
                                                           aaCTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACAGCCGGAATCGGCCAAACCTTAAAAAGAAGATGAACGAGAGAGAGTATACaa                                >  1:1786146/1‑99 (MQ=9)
                                                             cTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAGGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATc                             <  2:2207287/100‑1 (MQ=255)
                                                             cTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATc                             >  1:3407314/1‑100 (MQ=255)
                                                             cTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACACAGTGAAAGTATATAAATc                             >  1:2916324/1‑100 (MQ=255)
                                                               cAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGt                           >  2:3268859/1‑100 (MQ=255)
                                                               cAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGt                           >  2:113760/1‑100 (MQ=255)
                                                                 gggCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCa                         >  2:3885374/1‑100 (MQ=255)
                                                                 gggCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACAGCCGGAGTGGAACAAGGCTTAATATGGCGAAGCCTCAGGAAAAGTATATAAATCGACa                         >  1:3294855/1‑100 (MQ=16)
                                                                  ggCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCAc                        <  1:460272/100‑1 (MQ=255)
                                                                   gCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACt                       <  1:2825074/100‑1 (MQ=255)
                                                                      ggAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCg                    >  2:2336832/1‑100 (MQ=255)
                                                                        aaaGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCTGGAAACGGAGAAGGCTTAATATGAAGTTGACCTATGGGTAGTTTATAAATGGTGACGGCGat                  >  1:2140574/1‑100 (MQ=12)
                                                                        aaaGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGat                  >  1:3420960/1‑100 (MQ=255)
                                                                        aaaGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGat                  <  1:3539011/100‑1 (MQ=255)
                                                                         ttccgatctaCCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGata                 <  1:1018818/91‑1 (MQ=255)
                                                                         aaGTAAAAAACCGACGCAAAGTCGGTTTTTTACATTCGTAATCCGAAAAGGCTTAATATGACGATGACCCAGTAAAAGTATATAAATTGTCACTCCGat                  >  1:600076/1‑99 (MQ=25)
                                                                         aaGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGata                 <  2:1526499/100‑1 (MQ=255)
                                                                         aaGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGata                 <  2:3737591/100‑1 (MQ=255)
                                                                         aaGTAAAAAACCGACGCAAAGTCGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGata                 <  1:4210297/100‑1 (MQ=255)
                                                                            tAAAAAACCGACGCAACGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGatatat              <  2:248522/100‑1 (MQ=255)
                                                                            tAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGatatat              >  2:299131/1‑100 (MQ=255)
                                                                             aaaaaaCCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAACTCGTCAATGCGAtatata             >  2:1908722/1‑100 (MQ=39)
                                                                             aaaaaaCCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGAtatata             >  1:2277013/1‑100 (MQ=255)
                                                                             aaaaaaCCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGAtatata             >  2:2607528/1‑100 (MQ=255)
                                                                             aaaaaaCCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGAtatata             >  1:3821550/1‑100 (MQ=255)
                                                                             aaaaaaCCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGAtatata             >  2:1055907/1‑100 (MQ=255)
                                                                               aaaaCCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATAcc           <  2:4351419/100‑1 (MQ=255)
                                                                               aaaaCCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATAcc           >  2:2921519/1‑100 (MQ=255)
                                                                                 aaCCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGa         >  2:129753/1‑100 (MQ=255)
                                                                                  aCCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGataagatcggaa        >  2:1018820/1‑91 (MQ=255)
                                                                                   ccGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAg       <  2:3669685/100‑1 (MQ=255)
                                                                                   ccGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAg       <  1:3468120/100‑1 (MQ=255)
                                                                                   ccGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAg       <  2:1481164/100‑1 (MQ=255)
                                                                                   ccGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAg       <  1:1576833/100‑1 (MQ=255)
                                                                                     gACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAg       >  1:2545002/1‑98 (MQ=255)
                                                                                     gACGCAAAGTCGGTCTCTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAg       >  1:99318/1‑98 (MQ=255)
                                                                                       cGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGct   <  2:3231498/100‑1 (MQ=255)
                                                                                       cGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGct   >  2:1100424/1‑100 (MQ=255)
                                                                                        gCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctc  <  2:4407319/100‑1 (MQ=255)
                                                                                                  |                                                                                          
AGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTC  >  NC_000913/3648040‑3648228

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