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.96.5% 483.5 / 9.3 114intergenic (‑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 (2/2);  new base . (53/57);  total (55/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 131 total aligned reads displayed.

AGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCC  >  NC_000913/3648040‑3648233
                                                                                                  |                                                                                               
aGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtt                                                                                                >  1:3469978/1‑100 (MQ=255)
  tATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttt                                                                                              <  2:2956622/100‑1 (MQ=255)
  tATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttt                                                                                              <  1:576988/100‑1 (MQ=255)
   aTGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttt                                                                                             <  2:1536726/100‑1 (MQ=255)
    tGGTATTATCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTTGAAAGTAAAAAACCGACTTTGCGTCGGtttttt                                                                                            <  1:490339/100‑1 (MQ=21)
    tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                            >  2:63852/1‑100 (MQ=255)
    tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                            <  2:2520711/100‑1 (MQ=255)
    tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                            <  2:2909848/100‑1 (MQ=255)
     ggTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttta                                                                                           >  2:777089/1‑99 (MQ=255)
     ggTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttta                                                                                           <  2:544574/100‑2 (MQ=255)
     ggTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttta                                                                                           >  2:1799169/1‑99 (MQ=255)
      gTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttttac                                                                                          >  2:3764027/1‑98 (MQ=255)
       tATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttttt                                                                                           >  1:615696/1‑98 (MQ=255)
       tATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGCCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttttaca                                                                                         >  1:2667163/1‑97 (MQ=255)
        aTTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttttacat                                                                                        >  2:3221870/1‑96 (MQ=255)
          tAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATcc                                                                                     <  2:2636171/100‑1 (MQ=255)
           aGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATcc                                                                                     >  2:265517/1‑99 (MQ=255)
           aGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCg                                                                                    >  2:2548354/1‑100 (MQ=255)
           aGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCg                                                                                    <  2:93450/100‑1 (MQ=255)
           aGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCg                                                                                    >  2:3780436/1‑100 (MQ=255)
            gCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCgg                                                                                   >  2:1596314/1‑100 (MQ=255)
            gCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCgg                                                                                   <  1:2868890/100‑1 (MQ=255)
              aTTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGAt                                                                                 <  2:299914/100‑1 (MQ=255)
              aTTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGAt                                                                                 <  2:2171609/100‑1 (MQ=255)
               ttGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGAtt                                                                                >  2:1315025/1‑100 (MQ=255)
               ttGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGAtt                                                                                <  1:1200333/100‑1 (MQ=255)
                 gATGCGTTAGATGATGGCTATCTCATGCCAGTCAGAGCTACCAACTCAGTGCTGGAAACTAAAAAACCTACGCCAATTCGG‑TTTTTTACATCCGGATTCg                                                                              <  1:3468830/100‑1 (MQ=255)
                 gATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCTAAGTCGG‑TTTTTTACATCCGGATTCg                                                                              <  2:2310189/100‑1 (MQ=255)
                 gATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCg                                                                              <  2:488590/100‑1 (MQ=255)
                 gATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCg                                                                              >  2:3621901/1‑100 (MQ=255)
                 gATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCg                                                                              >  2:2727622/1‑100 (MQ=255)
                   tGCGTTAGATGATGGCGATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGa                                                                            <  1:1369203/100‑1 (MQ=255)
                      gTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTAACTTc                                                                                      >  1:3260243/1‑88 (MQ=255)
                       ttAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAg                                                                        <  2:2449950/100‑1 (MQ=255)
                       ttAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAg                                                                        <  2:3017763/100‑1 (MQ=255)
                            tgatgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCttn                                                                   <  1:3725702/100‑2 (MQ=255)
                            tgatgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTa                                                                   >  2:1889124/1‑100 (MQ=255)
                             gatgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACAGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTaa                                                                  >  1:3109219/1‑100 (MQ=255)
                              atgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAat                                                                 <  2:2216174/100‑1 (MQ=255)
                              atgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAat                                                                 <  2:795837/100‑1 (MQ=255)
                               tgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTAAATTCGGATTCCAAAAAGGCATAat                                                                 >  1:385679/1‑99 (MQ=14)
                                ggCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAata                                                                >  2:3123163/1‑99 (MQ=255)
                                 gCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATg                                                              <  2:1938827/100‑1 (MQ=255)
                                  cTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGa                                                             >  1:1091292/1‑100 (MQ=255)
                                  cTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGa                                                             >  1:1376493/1‑100 (MQ=255)
                                   tATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGa                                                             <  1:311570/100‑1 (MQ=255)
                                   tATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGat                                                            <  1:353599/100‑2 (MQ=255)
                                   tATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGAc                                                            <  2:2350396/100‑1 (MQ=255)
                                      ctcACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGAt                                                         <  2:560229/100‑1 (MQ=255)
                                      ctcACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑GTTTTTAGATCCGGATTCGGACAAGGCTTAATATGACGAt                                                         >  1:2663976/1‑100 (MQ=21)
                                       tgaccccAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATg                                                        <  2:2983080/95‑1 (MQ=255)
                                        cACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGa                                                       >  2:810753/1‑100 (MQ=255)
                                        cACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGa                                                       >  2:2431837/1‑100 (MQ=255)
                                            ccAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGCCCCa                                                   >  2:1409071/1‑100 (MQ=255)
                                             cAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAACCCTTAATATGACTATGcccc                                                    >  1:2012446/1‑99 (MQ=18)
                                             cAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAg                                                  >  2:89165/1‑100 (MQ=255)
                                               gTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑ATTTTTACATACGGATTCGGACAAGGCTTAATATGACGATGACCCAGTg                                                >  1:1106368/1‑100 (MQ=25)
                                                tCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGa                                               <  2:690816/100‑1 (MQ=255)
                                                   gagCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAg                                            <  2:3726478/100‑1 (MQ=255)
                                                   gagCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAg                                            <  2:2417800/100‑1 (MQ=255)
                                                   gagCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAg                                            <  2:2981498/100‑1 (MQ=255)
                                                   gagCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAg                                            <  1:26101/100‑1 (MQ=255)
                                                     gCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGta                                          >  2:2914454/1‑100 (MQ=255)
                                                     gCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGta                                          >  2:2810131/1‑100 (MQ=255)
                                                      ccaccaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtat                                         <  1:3621873/100‑1 (MQ=255)
                                                       caccaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtata                                        <  2:3797758/100‑1 (MQ=255)
                                                        accaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtatat                                       <  2:3480651/100‑1 (MQ=255)
                                                        accaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtatat                                       <  2:1836390/100‑1 (MQ=255)
                                                           aaCTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATaaa                                    <  1:2760767/100‑1 (MQ=255)
                                                           aaCTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATaa                                     <  1:1075619/99‑1 (MQ=255)
                                                            aCTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACCTTCCTAATCGGAAAAGCCTTAATATGAAGAAGCCCCAGCAAAAGTATaaaa                                     >  1:271185/1‑98 (MQ=9)
                                                            aCTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTACTATGACGATGACCCAGTGAAAGTATATaaa                                    >  2:1192877/1‑99 (MQ=255)
                                                            aCTCAGGGCTGGAAAATAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAAt                                   <  2:2759773/100‑1 (MQ=255)
                                                             cTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATc                                  <  2:1506260/100‑1 (MQ=255)
                                                               cAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGt                                >  1:1537410/1‑100 (MQ=255)
                                                               cAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGt                                >  2:1328614/1‑100 (MQ=255)
                                                                aGGGCTGGATAGTTATTAATCGACGCAAAGCCGG‑TTTTTTATATCCGGATGCAGACAAGTCTCATTATGACGATGACCCAGTGAAAGTATATAAATCGTc                               <  1:996611/100‑1 (MQ=255)
                                                                 gggCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCa                              <  1:1799796/100‑1 (MQ=255)
                                                                   gCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGATTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACt                            >  2:1742663/1‑100 (MQ=255)
                                                                    cTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTg                           <  1:2836231/100‑1 (MQ=255)
                                                                      ggAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCg                         >  2:3449051/1‑100 (MQ=255)
                                                                      ggAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCg                         >  2:1670837/1‑100 (MQ=255)
                                                                       ggatgTAAAATACCGACTCTGCGTCGG‑TTTTTTACATCCGGATTCGGTCAAGGCTTAATATGACGATGAACCAGTGAAAGTATATAAATTGTGGCTGCGa                        >  1:2805067/5‑100 (MQ=255)
                                                                        aaaGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGat                       <  2:3569875/100‑1 (MQ=255)
                                                                        aaaGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGat                       >  2:366497/1‑100 (MQ=255)
                                                                        aaaGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGat                       <  2:1376498/100‑1 (MQ=255)
                                                                             aaaaaaCCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGAtatata                  >  2:3131907/1‑100 (MQ=255)
                                                                             aaaaaaCCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGAtatata                  <  2:2344119/100‑1 (MQ=255)
                                                                             aaaaaaCCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGAtatata                  >  1:858880/1‑100 (MQ=255)
                                                                               aaaaCCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATAcc                >  1:2889885/1‑100 (MQ=255)
                                                                               aaaaCCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATAcc                >  2:1374763/1‑100 (MQ=255)
                                                                                  aCCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGAGGATGACCCAGGGAAAGTATATAAATCGTCACTGCGATATATACCg               >  2:935615/1‑98 (MQ=255)
                                                                                   ccGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAg            <  2:2081197/100‑1 (MQ=255)
                                                                                     gACGCAAAGGCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGt           <  2:3613077/100‑1 (MQ=255)
                                                                                       cGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGct        >  2:666196/1‑100 (MQ=255)
                                                                                        gCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAATGCTTAATATGACGATTACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctc       <  1:3331649/100‑1 (MQ=255)
                                                                                          ataGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctccc     <  1:870547/98‑1 (MQ=255)
                                                                                          aaaGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctccc     <  2:1369813/100‑1 (MQ=255)
                                                                                           aaGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACTGAAGTGctccct    >  1:378176/1‑100 (MQ=255)
                                                                                           aaGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctccc     <  1:3214179/99‑1 (MQ=255)
                                                                                                  |                                                                                               
AGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCC  >  NC_000913/3648040‑3648233

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: ‑