Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,648,144 (T)7→6 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.96.8% 540.4 / 10.9 124intergenic (‑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 (4/0);  new base . (57/63);  total (61/63)
Fisher's exact test for biased strand distribution p-value = 5.56e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

Only 100 of 139 total aligned reads displayed.

GTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGC  >  NC_000913/3648041‑3648235
                                                                                                 |                                                                                                 
gTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttt                                                                                                 <  2:557974/100‑1 (MQ=255)
gTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttt                                                                                                 >  1:1498019/1‑100 (MQ=255)
   tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCATTCCAGTCAGAGCCACCAACCCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                              <  2:1629826/100‑1 (MQ=255)
   tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAGACCGACGCAAAGTCGGtttttt                                                                                              <  2:1735498/100‑1 (MQ=255)
   tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACTTTGCGTCGGtttttt                                                                                              <  1:1958855/100‑1 (MQ=25)
   tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                              <  1:1463971/100‑1 (MQ=255)
   tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                              <  2:189245/100‑1 (MQ=255)
   tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                              <  1:1275545/100‑1 (MQ=255)
   tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                              <  2:591229/100‑1 (MQ=255)
    ggTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttta                                                                                             <  2:1412195/100‑2 (MQ=255)
      tATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGatttttaca                                                                                           >  1:2786889/1‑97 (MQ=255)
       aTTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttttacat                                                                                          >  2:1250994/1‑96 (MQ=255)
       aTTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttttacat                                                                                          <  2:1769607/100‑5 (MQ=255)
        ttAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATc                                                                                        <  1:928089/100‑1 (MQ=255)
        ttAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATc                                                                                        <  2:2197704/100‑1 (MQ=255)
         tAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATcc                                                                                       <  2:1497106/100‑1 (MQ=255)
         tAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATcc                                                                                       >  2:530934/1‑100 (MQ=255)
         tAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATcc                                                                                       >  2:217655/1‑100 (MQ=255)
          aGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCg                                                                                      >  2:2660182/1‑100 (MQ=255)
          aGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCg                                                                                      >  1:35801/1‑100 (MQ=255)
          aGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCg                                                                                      >  1:2906527/1‑100 (MQ=255)
          aGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCg                                                                                      >  2:1881850/1‑100 (MQ=255)
           gCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCgg                                                                                     >  2:657286/1‑100 (MQ=255)
           gCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCgg                                                                                     >  2:2020826/1‑100 (MQ=255)
            cATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGa                                                                                    >  2:2959970/1‑100 (MQ=255)
                gATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCg                                                                                <  1:1341268/100‑1 (MQ=255)
                 aTGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTc                                                                                 <  2:1013811/98‑1 (MQ=255)
                 aTGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCgg                                                                               <  1:2161751/100‑1 (MQ=255)
                  tGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGa                                                                              >  2:1940512/1‑100 (MQ=255)
                  tGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGa                                                                              <  1:976052/100‑1 (MQ=255)
                      ttAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAg                                                                          >  2:236877/1‑100 (MQ=255)
                       tAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAgg                                                                         <  1:1632148/100‑1 (MQ=255)
                       tAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAgg                                                                         <  1:718680/100‑1 (MQ=255)
                        aGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGc                                                                        >  1:2933357/1‑100 (MQ=255)
                         gatgatGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCt                                                                       >  2:1036864/1‑100 (MQ=255)
                           tgatgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCtt                                                                      <  2:1296687/99‑1 (MQ=255)
                           tgatgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTa                                                                     >  1:1481905/1‑100 (MQ=255)
                             atgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAat                                                                   <  1:1129243/100‑1 (MQ=255)
                                gCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCTGATTCGGACAAGGCTTAata                                                                  <  2:1818394/98‑1 (MQ=255)
                                gCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAatatn                                                                <  2:2518984/100‑2 (MQ=255)
                                 cTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGa                                                               <  2:2329410/100‑1 (MQ=255)
                                  tATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGAc                                                              >  2:1995146/1‑100 (MQ=255)
                                  tATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGAc                                                              <  2:321615/100‑1 (MQ=255)
                                   aTCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGAc                                                              <  2:492678/99‑1 (MQ=255)
                                   aTCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGAc                                                              <  1:2350673/99‑1 (MQ=255)
                                      tcACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATg                                                          <  1:1852098/100‑1 (MQ=255)
                                      tcACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATg                                                          >  2:1625952/1‑100 (MQ=255)
                                      tcACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATg                                                          >  2:2231752/1‑100 (MQ=255)
                                         cTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAcc                                                       >  2:2950764/1‑100 (MQ=255)
                                         cTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAcc                                                       >  1:1314034/1‑100 (MQ=255)
                                          tCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAccc                                                      >  1:1531146/1‑100 (MQ=255)
                                          tCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAccc                                                      >  1:765578/1‑100 (MQ=255)
                                            cAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAg                                                    <  1:399010/100‑1 (MQ=255)
                                              gTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTg                                                  >  2:2590394/1‑100 (MQ=255)
                                                cAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGaa                                                <  2:2505568/100‑1 (MQ=255)
                                                cAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGaa                                                >  2:68012/1‑100 (MQ=255)
                                                 agagCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGaaa                                               >  1:1969002/1‑100 (MQ=255)
                                                  gagCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAg                                              <  1:830938/100‑1 (MQ=255)
                                                   agCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGt                                             <  2:1498010/100‑1 (MQ=255)
                                                    gCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGta                                            >  1:784615/1‑100 (MQ=255)
                                                    gCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGta                                            >  1:515697/1‑100 (MQ=255)
                                                      caccaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtata                                          >  1:595939/1‑100 (MQ=255)
                                                      caccaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtata                                          <  1:1759996/100‑1 (MQ=255)
                                                       accaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtatat                                         <  2:2978135/100‑1 (MQ=255)
                                                       accaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtatat                                         >  1:2436298/1‑100 (MQ=255)
                                                        ccaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtatata                                        <  2:2943392/100‑1 (MQ=255)
                                                         caACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtatata                                        >  1:2856183/1‑100 (MQ=255)
                                                          aaCTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATaaa                                      <  1:1625961/100‑1 (MQ=255)
                                                            cTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAAt                                     >  2:720580/1‑100 (MQ=255)
                                                            cTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATc                                    >  1:850097/1‑100 (MQ=255)
                                                            cTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATc                                    >  1:622613/1‑100 (MQ=255)
                                                            cTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATc                                    <  2:2463969/100‑1 (MQ=255)
                                                            cTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGGTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATc                                    >  1:2093188/1‑100 (MQ=255)
                                                             tCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCg                                   <  1:2463844/100‑1 (MQ=255)
                                                             tCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCg                                   <  2:227523/100‑1 (MQ=255)
                                                               aGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTc                                 <  2:2163158/100‑1 (MQ=255)
                                                                 ggCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCAc                               <  1:2713559/100‑1 (MQ=255)
                                                                    tGGAAAGTAAAAAACCGACAGATAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGGATATAAATCGTCAc                               <  2:229353/97‑1 (MQ=38)
                                                                     ggAAATTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCg                           <  2:527467/100‑1 (MQ=255)
                                                                     ggAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCg                           >  1:1331600/1‑100 (MQ=255)
                                                                     ggAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCg                           <  1:1251000/100‑1 (MQ=255)
                                                                      gAAAGTAAAAAACCGACGCAAAGTCGGTCTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGa                          <  2:1453086/100‑1 (MQ=255)
                                                                       aaaGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGat                         <  2:765575/100‑1 (MQ=255)
                                                                       aaaGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGat                         <  2:2913286/100‑1 (MQ=255)
                                                                           tAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGatatat                     >  1:897117/1‑100 (MQ=255)
                                                                            aaaaaaCCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGAtatata                    <  2:673150/100‑1 (MQ=255)
                                                                              aaaaCCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATAcc                  <  2:1302177/100‑1 (MQ=255)
                                                                               aaaCCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCg                 >  1:1308056/1‑100 (MQ=255)
                                                                                 aCCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGaa               <  1:1179647/100‑1 (MQ=255)
                                                                                 aCCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGaa               >  1:1493727/1‑100 (MQ=255)
                                                                                  ccGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAg              <  1:1881861/100‑1 (MQ=255)
                                                                                  ccGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAg              >  1:1959235/1‑100 (MQ=255)
                                                                                  ccGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCAATATATACCGAAg              <  2:35800/100‑1 (MQ=255)
                                                                                    gACGCAAAGTCGGTTTTTTACATCCCTGTTCCAGCAAGCCTTAACGTGACGATGTCCCGGTGGAAGTATATAAATCGTCACTGCGATATATACCGAAg              >  2:2776641/1‑98 (MQ=255)
                                                                                     aCGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGc           <  1:236879/100‑1 (MQ=255)
                                                                                       gCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctc         <  2:2303002/100‑1 (MQ=255)
                                                                                        cAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctcc        >  1:2091601/1‑100 (MQ=255)
                                                                                          aaGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctccct      <  1:68013/100‑1 (MQ=255)
                                                                                            gTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctccctcc    >  2:2101463/1‑100 (MQ=255)
                                                                                               ggTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGc  >  2:1511446/1‑100 (MQ=255)
                                                                                                 |                                                                                                 
GTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGC  >  NC_000913/3648041‑3648235

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