Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,648,042 T→C intergenic (‑254/‑486) 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,0420TC96.7% 316.7 / ‑0.8 93intergenic (‑254/‑486)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/2);  new base C (45/44);  total (47/46)
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.65e-01

ATGCGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGT  >  NC_000913/3647946‑3648138
                                                                                                |                                                                                                
aTGCGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATg                                                                                               >  1:781237/1‑100 (MQ=255)
 tGCGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATgg                                                                                              <  1:463817/100‑1 (MQ=255)
    gcATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTAt                                                                                           >  1:1743866/1‑100 (MQ=255)
    gcATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCACGGTAt                                                                                           >  2:1211499/1‑100 (MQ=255)
      aTTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTa                                                                                         <  2:396565/100‑1 (MQ=255)
      aTTGCCAGAGTGAAAATGCCTCAGAGCTAACCGGATATTCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTa                                                                                         <  1:1307301/100‑1 (MQ=255)
       ttGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAg                                                                                        <  1:452879/100‑1 (MQ=255)
        tGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGc                                                                                       <  1:87847/100‑1 (MQ=255)
        tGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGc                                                                                       <  2:1070919/100‑1 (MQ=255)
         gCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCa                                                                                      >  2:1524706/1‑100 (MQ=255)
         gCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCa                                                                                      <  1:433313/100‑1 (MQ=255)
         gCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCa                                                                                      >  2:1661407/1‑100 (MQ=255)
         gCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCa                                                                                      <  2:787689/100‑1 (MQ=255)
           cAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCAtt                                                                                    >  1:991591/1‑100 (MQ=255)
           cAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGAATGGTATTAGCAtt                                                                                    >  2:380461/1‑100 (MQ=255)
           cAGAGTGAAAATGCCTCAGAGCTAACTGGAGAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATTGTATTAGCAtt                                                                                    <  1:839193/100‑1 (MQ=255)
               gTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATg                                                                                <  2:557756/100‑1 (MQ=255)
                tGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGc                                                                               <  2:645284/100‑1 (MQ=255)
                    cgatcttCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTa                                                                           <  2:599267/94‑1 (MQ=255)
                    aaTGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTa                                                                           <  2:26852/100‑1 (MQ=255)
                    aaTGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTa                                                                           <  2:205166/100‑1 (MQ=255)
                    aaTGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTa                                                                           <  1:1137048/100‑1 (MQ=255)
                       gCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAgat                                                                        <  1:1503846/100‑1 (MQ=255)
                        ccTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAgatg                                                                       >  2:1395850/1‑100 (MQ=255)
                        ccTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAgatg                                                                       >  2:1395957/1‑100 (MQ=255)
                          tCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCCTGGTATTAGCATTGATGCGTTAgatgat                                                                     >  1:1684504/1‑100 (MQ=255)
                          tCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTaagatcg                                                                     >  1:599265/1‑94 (MQ=255)
                           cAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGatgatg                                                                    >  2:948285/1‑100 (MQ=255)
                           cAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGatgatg                                                                    <  1:415462/100‑1 (MQ=255)
                           cAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGatgatg                                                                    >  1:113989/1‑100 (MQ=255)
                           cAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGatgatg                                                                    <  1:1118200/100‑1 (MQ=255)
                           cAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGatgatg                                                                    <  1:147454/100‑1 (MQ=255)
                            agagCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATgg                                                                   <  1:1701624/100‑1 (MQ=255)
                            agagCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATgg                                                                   <  2:124347/100‑1 (MQ=255)
                            agagCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATgg                                                                   <  2:255252/100‑1 (MQ=255)
                            agagCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATgg                                                                   <  1:259699/100‑1 (MQ=255)
                               tcTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTa                                                                <  2:558536/99‑1 (MQ=255)
                               gCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTa                                                                <  1:1139603/100‑1 (MQ=255)
                                 tAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTAtc                                                              >  2:248598/1‑100 (MQ=255)
                                 tAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTGGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTAtc                                                              <  1:1885378/100‑1 (MQ=255)
                                  aaCTGTAGAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTAtc                                                              <  1:942000/99‑1 (MQ=255)
                                  aaCTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTaaga                                                             >  1:558534/1‑97 (MQ=255)
                                  aaCTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTAtct                                                             >  2:1784495/1‑100 (MQ=255)
                                  aaCTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTAtct                                                             <  2:1477594/100‑1 (MQ=255)
                                   aCTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTAtctc                                                            <  2:1693475/100‑1 (MQ=255)
                                     tGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGTATTGATGCGTTAGATGATGGCTATCTCAc                                                          >  2:511279/1‑100 (MQ=255)
                                     tGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCAc                                                          <  1:39471/100‑1 (MQ=255)
                                     tGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCAc                                                          <  1:1593647/100‑1 (MQ=255)
                                     tGGATAATCATACAGTACATGCAGGTCATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCAc                                                          >  2:1186749/1‑100 (MQ=255)
                                      ggATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACt                                                         >  1:609425/1‑100 (MQ=255)
                                      ggATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACt                                                         >  1:1361536/1‑100 (MQ=255)
                                        aTAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTcc                                                       >  2:1703236/1‑100 (MQ=255)
                                           aTCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGt                                                    >  1:1223516/1‑100 (MQ=255)
                                            tCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTc                                                   >  1:90559/1‑100 (MQ=255)
                                             cATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCa                                                  >  1:1517062/1‑100 (MQ=255)
                                                aCAGTACATGCAGGTTATAAACCCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCagag                                               <  2:1074211/100‑1 (MQ=255)
                                                 cAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGc                                              >  1:301239/1‑100 (MQ=255)
                                                  aGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGcc                                             >  2:1936339/1‑100 (MQ=255)
                                                        tGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAAc                                       >  2:1098659/1‑100 (MQ=255)
                                                         gCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACt                                      >  2:1056304/1‑100 (MQ=255)
                                                          cAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTc                                     >  1:1852122/1‑100 (MQ=255)
                                                           aGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCa                                    <  1:1108152/100‑1 (MQ=255)
                                                             gTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAgg                                  <  2:1743877/100‑1 (MQ=255)
                                                                 taAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTg                              >  2:628937/1‑100 (MQ=255)
                                                                 taAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTg                              >  2:215050/1‑100 (MQ=255)
                                                                 taAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTg                              >  1:514108/1‑100 (MQ=255)
                                                                  aaaaCCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTTATGCTTTAGATGATGGCTATCTCACTCCAGTCATAGCCACCAACTCAGGGCTgg                             >  1:646791/1‑100 (MQ=38)
                                                                  aaaaCCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTgg                             <  1:1524698/100‑1 (MQ=255)
                                                                   aaaCCAGCACGTGCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGa                            >  2:1278442/1‑100 (MQ=255)
                                                                   aaaCCAGCACGTCCTTTCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGa                            <  1:1872640/100‑1 (MQ=255)
                                                                   aaaCCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGa                            >  1:379075/1‑100 (MQ=255)
                                                                   aaaCCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGa                            >  1:1932905/1‑100 (MQ=255)
                                                                     aCCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGaaa                          >  1:1708808/1‑100 (MQ=255)
                                                                      ccAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAg                         <  2:631736/100‑1 (MQ=255)
                                                                      ccAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAg                         <  1:830283/100‑1 (MQ=255)
                                                                       cAGCACGTCCTTGCATTAGTTTCAGCATGGTATTAGCATTGATGCGTTCGTTGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGaaan                         <  1:425269/99‑2 (MQ=37)
                                                                        aGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTa                       >  2:1103523/1‑100 (MQ=255)
                                                                         gCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTaa                      <  1:1355729/100‑1 (MQ=255)
                                                                               ccTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAAc                 <  2:380442/99‑1 (MQ=255)
                                                                               ccTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAAcc                >  1:132326/1‑100 (MQ=255)
                                                                               ccTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAAcc                >  2:44491/1‑100 (MQ=255)
                                                                               ccTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATATCACTCCAGTCAGAGCCACCAACACAGGGCTGGAAAGTAAAAAAcc                <  1:1061979/100‑1 (MQ=255)
                                                                                  tGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGAc             >  2:575944/1‑100 (MQ=255)
                                                                                   gCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACg            >  1:694123/1‑100 (MQ=255)
                                                                                   gCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCACCTCAGGGCTGGAAAGTAAAAAACCGACg            <  2:1359520/100‑1 (MQ=255)
                                                                                   gCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACg            >  1:927292/1‑100 (MQ=255)
                                                                                    cAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGc           >  1:980160/1‑100 (MQ=255)
                                                                                      aTAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCaa         <  2:465869/100‑1 (MQ=255)
                                                                                      aTAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCaa         <  2:1361544/100‑1 (MQ=255)
                                                                                        aGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAg       <  1:1703225/100‑1 (MQ=255)
                                                                                         gTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGt      >  1:1074056/1‑100 (MQ=255)
                                                                                          tttCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAgac     >  1:1183468/1‑98 (MQ=255)
                                                                                             cAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGt  >  1:1496175/1‑100 (MQ=255)
                                                                                                |                                                                                                
ATGCGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGT  >  NC_000913/3647946‑3648138

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