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.5% 354.5 / 10.6 88intergenic (‑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 . (44/40);  total (48/40)
Fisher's exact test for biased strand distribution p-value = 1.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 106 total aligned reads displayed.

CAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCC  >  NC_000913/3648039‑3648233
                                                                                                   |                                                                                               
cAGTATGTTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGt                                                                                                 <  1:1131937/100‑1 (MQ=255)
cAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGt                                                                                                 <  1:1868516/100‑1 (MQ=255)
  gTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttt                                                                                               >  1:1518078/1‑100 (MQ=255)
   tATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttt                                                                                              <  1:1835773/100‑1 (MQ=255)
    aTGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttt                                                                                             <  1:2153315/100‑1 (MQ=255)
    aTGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttt                                                                                             >  2:1459477/1‑100 (MQ=255)
    aTGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttt                                                                                             <  2:1802972/100‑1 (MQ=255)
     tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCCGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                            <  1:2036150/100‑1 (MQ=255)
     tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                            <  1:1869080/100‑1 (MQ=255)
      ggTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttta                                                                                           <  2:1882510/100‑2 (MQ=255)
       gTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttttac                                                                                          <  2:149385/100‑3 (MQ=255)
       gTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACAGACGCAAAGTCGGttttttac                                                                                          >  1:2035221/1‑98 (MQ=255)
        tATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttttaca                                                                                         <  1:92142/100‑4 (MQ=255)
        tATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttttaca                                                                                         >  1:1602445/1‑97 (MQ=255)
         aTTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttttacat                                                                                        >  2:1367673/1‑96 (MQ=255)
          ttAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCACAGCCGGtagtgta                                                                                           >  1:484843/1‑90 (MQ=21)
          ttAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATc                                                                                      <  2:2596383/100‑1 (MQ=255)
          ttAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATc                                                                                      >  2:2548732/1‑100 (MQ=255)
           tAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttttactttc                                                                                      <  2:1437546/100‑7 (MQ=255)
           tAGCATTGAGGCGTTAGATAATGGCTATCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttttt                                                                                           <  1:2548721/95‑1 (MQ=255)
            aGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTc                                                                                      >  1:1621989/1‑99 (MQ=255)
            aGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATcc                                                                                     >  1:311381/1‑100 (MQ=255)
            aGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATcca                                                                                    >  2:243585/1‑99 (MQ=255)
             gCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCgg                                                                                   >  2:1378956/1‑100 (MQ=255)
                 tGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTc                                                                               <  1:242104/100‑1 (MQ=255)
                  gATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCg                                                                              <  2:651172/100‑1 (MQ=255)
                    tGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGa                                                                            <  1:993755/100‑1 (MQ=255)
                      cGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACa                                                                          >  2:1506845/1‑100 (MQ=255)
                      cGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACa                                                                          >  2:321611/1‑100 (MQ=255)
                      cGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACa                                                                          >  2:1336864/1‑100 (MQ=255)
                        ttAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAg                                                                        <  1:341556/100‑1 (MQ=255)
                         tAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAgg                                                                       <  1:321609/100‑1 (MQ=255)
                          aGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGc                                                                      <  2:353915/100‑1 (MQ=255)
                             tgatgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTa                                                                   >  1:2177007/1‑100 (MQ=255)
                             tgatgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCAGATTCGGACAAGGCTTa                                                                   <  1:243583/100‑1 (MQ=255)
                               atgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAat                                                                 <  1:2113039/100‑1 (MQ=255)
                                tgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAata                                                                >  1:1995032/1‑100 (MQ=255)
                                  gCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATg                                                              >  2:1551255/1‑100 (MQ=255)
                                    tATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGAc                                                            <  1:228490/100‑1 (MQ=255)
                                    tATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGAc                                                            <  2:1204030/100‑1 (MQ=255)
                                     aTCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACg                                                           >  1:2877428/1‑100 (MQ=255)
                                     aTCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACg                                                           <  2:795239/100‑1 (MQ=255)
                                      tctcACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGa                                                          >  2:2791449/1‑100 (MQ=255)
                                      tctcACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGa                                                          >  2:2397218/1‑100 (MQ=255)
                                        tcACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATg                                                        >  2:2364948/1‑100 (MQ=255)
                                         cACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGa                                                       >  2:2521018/1‑100 (MQ=255)
                                           cTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAcc                                                     >  1:1157401/1‑100 (MQ=255)
                                           cTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAcc                                                     <  1:1502481/100‑1 (MQ=255)
                                           cTCCAGTCAGAGCCACCAACTCAGGACTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAcc                                                     >  1:123272/1‑100 (MQ=255)
                                             ccAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAccc                                                    >  2:2019201/1‑99 (MQ=255)
                                             ccAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCa                                                   >  2:1665602/1‑100 (MQ=255)
                                              cAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAg                                                  <  1:614286/100‑1 (MQ=255)
                                               aGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGt                                                 <  2:638923/100‑1 (MQ=255)
                                               aGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGt                                                 <  1:1378954/100‑1 (MQ=255)
                                                    gagCCACAAACTCAGGTCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAg                                            <  1:1425831/100‑1 (MQ=25)
                                                      gCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGta                                          >  2:522522/1‑100 (MQ=255)
                                                      gCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGta                                          <  1:1336862/100‑1 (MQ=255)
                                                      ctcaccaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGta                                          <  1:2202670/98‑1 (MQ=255)
                                                       ccaccaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCAGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGt                                           <  1:2841856/98‑1 (MQ=255)
                                                        caccaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGATGATGACCCAGTGAAAGtata                                        <  1:1784019/100‑1 (MQ=255)
                                                        caccaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtata                                        <  1:68562/100‑1 (MQ=255)
                                                        caccaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtaag                                        >  2:2202678/1‑98 (MQ=255)
                                                         accaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtatat                                       >  2:2058544/1‑100 (MQ=255)
                                                          ccaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtatata                                      <  1:960169/100‑1 (MQ=255)
                                                          ccaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtatata                                      <  1:607645/100‑1 (MQ=255)
                                                           caACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATaa                                     <  1:349286/100‑1 (MQ=255)
                                                               tCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCg                                 >  2:588645/1‑100 (MQ=255)
                                                               tCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCg                                 >  2:497306/1‑100 (MQ=255)
                                                                cAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATTCGGAATCGCACAAGGCTTAATATGACTATGAACCAATAAAATtataca                                      >  1:1145977/1‑92 (MQ=18)
                                                                cAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGt                                >  2:585188/1‑100 (MQ=255)
                                                                  gggCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCa                              <  1:2747686/100‑1 (MQ=255)
                                                                  gggCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGATTAATATGACTATGACCCAGTGAAAGTATATaaa                                    >  1:55575/1‑94 (MQ=25)
                                                                     cTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTg                           >  2:1922533/1‑100 (MQ=255)
                                                                      tGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGc                          <  1:169294/100‑1 (MQ=255)
                                                                       ggAAAGTCAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGc                          >  2:1528736/1‑100 (MQ=255)
                                                                       ggAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGAATCGGACAAGGCTTAATATGAAGATGACCCAGTGAAAGTATATAAATCGTCACTGCg                         >  1:2242137/1‑100 (MQ=38)
                                                                        ngatgTAAAAAACCGACTTTGCGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGa                        >  2:2130550/5‑100 (MQ=37)
                                                                        gAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGc                          <  1:700147/98‑1 (MQ=255)
                                                                         aaaGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGat                       <  1:1390041/100‑1 (MQ=255)
                                                                            tctaaaaaCCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGatata                    <  1:2310962/97‑1 (MQ=255)
                                                                             tAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGatatat                   >  2:781881/1‑100 (MQ=255)
                                                                             tAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGGGAAAGTATATAAATCGTCACTGCGatatat                   <  2:1932388/100‑1 (MQ=255)
                                                                              aaaaaaCCGACTTTGCGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGCtatata                  >  1:706745/1‑100 (MQ=37)
                                                                              aaaaaaCCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGAtatata                  >  2:829885/1‑100 (MQ=255)
                                                                              aaaaaaCCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGAtatata                  >  2:573346/1‑100 (MQ=255)
                                                                               naaaaCCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTACTATGACGATGACCCAGTAAAAGTATCTACATCGTAAATGCGaaaaaaa                  >  1:2401114/2‑93 (MQ=25)
                                                                               aaaaaCCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGatataaga                 >  2:2310972/1‑97 (MQ=255)
                                                                                aaaaCCGACGCAAAGTCGGTTTTTTACATTCGGATTCTGACAATTCTTTATATGACGTTGACCCAGTGAAAGTATATAAATCGTCACTGCGatatat                   >  2:1482840/1‑97 (MQ=37)
                                                                                aaaaCCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATAcc                <  2:2562487/100‑1 (MQ=255)
                                                                                  aaCCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGa              <  2:335835/100‑1 (MQ=255)
                                                                                    ccGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAg            >  2:204966/1‑100 (MQ=255)
                                                                                     cGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAg            >  1:1670912/1‑99 (MQ=255)
                                                                                     cGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGt           >  1:1670725/1‑100 (MQ=255)
                                                                                      gACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTg          <  2:1973420/100‑1 (MQ=255)
                                                                                      gACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTg          <  1:1985711/100‑1 (MQ=255)
                                                                                          cAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctcc      >  2:2662833/1‑100 (MQ=255)
                                                                                           aaaGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctccc     >  1:1608948/1‑100 (MQ=255)
                                                                                           aaaGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctccc     <  2:1273802/100‑1 (MQ=255)
                                                                                            aaGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctccct    >  2:2873406/1‑100 (MQ=255)
                                                                                              gTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctccctcc  >  1:2248656/1‑100 (MQ=255)
                                                                                                   |                                                                                               
CAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCC  >  NC_000913/3648039‑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: