Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,337,024 G→A 100% S3L (TCG→TTG)  adiC ← arginine:agmatine antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,337,0240GA99.1% 353.7 / ‑4.6 105S3L (TCG→TTG) adiCarginine:agmatine antiporter
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (55/49);  total (56/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

Only 100 of 105 total aligned reads displayed.

AGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGG  >  NC_000913/4336928‑4337120
                                                                                                |                                                                                                
aGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAg                                                                                               <  1:1513168/100‑1 (MQ=255)
 ggTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGa                                                                                              <  2:840349/100‑1 (MQ=255)
  gTTTGCAGGTAACAGAACAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGAc                                                                                             <  1:616133/100‑1 (MQ=255)
  gTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGAc                                                                                             <  2:1196142/100‑1 (MQ=255)
     tGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATa                                                                                          >  1:1442873/1‑100 (MQ=255)
        agttAACGGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGAGATTACGCCCACTTTGTGAGCATCAGCATCCAAAGACATAAtt                                                                                       <  1:1387902/97‑1 (MQ=255)
        aGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAAtt                                                                                       <  1:876960/100‑1 (MQ=255)
         ggTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTa                                                                                      >  2:404218/1‑100 (MQ=255)
           tAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTaaa                                                                                    >  1:1790554/1‑100 (MQ=255)
            aaCAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAAc                                                                                   >  2:398358/1‑100 (MQ=255)
             tctgAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACt                                                                                  <  2:615104/97‑1 (MQ=255)
              cAGAAAAACATCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTc                                                                                 >  2:1812798/1‑100 (MQ=255)
              cAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTc                                                                                 >  2:562002/1‑100 (MQ=255)
              cAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTc                                                                                 >  2:1085609/1‑100 (MQ=255)
              cAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTc                                                                                 >  2:2383814/1‑100 (MQ=255)
               aGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTcc                                                                                <  1:915236/100‑1 (MQ=255)
                gAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACtaga                                                                               >  1:615101/1‑97 (MQ=255)
                gAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCt                                                                               >  1:770005/1‑100 (MQ=255)
                  aaaaCACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGc                                                                             >  2:2544646/1‑100 (MQ=255)
                   aaaCACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATTAGCATCCAAAGACATAATTAAACTCCTGc                                                                             >  1:1276089/1‑99 (MQ=255)
                   aaaCACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCg                                                                            >  2:934878/1‑100 (MQ=255)
                    aaCACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGa                                                                           <  2:484019/100‑1 (MQ=255)
                       acCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAgg                                                                        <  2:569486/100‑1 (MQ=255)
                        ccTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAACATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGc                                                                       >  1:1967112/1‑100 (MQ=255)
                          ttaCCCCATAATATTCCCTGACACCAACAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAgg                                                                        <  1:1349408/95‑1 (MQ=255)
                           gACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAg                                                                    <  1:1229822/100‑1 (MQ=255)
                              cccATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCtt                                                                 <  2:2434085/100‑1 (MQ=255)
                               ccATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTc                                                                >  2:1641678/1‑100 (MQ=255)
                                cATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCg                                                               <  2:694495/100‑1 (MQ=255)
                                  taataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTg                                                             <  2:173788/100‑1 (MQ=255)
                                  taataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTg                                                             <  2:173689/100‑1 (MQ=255)
                                  taataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTg                                                             <  2:1640687/100‑1 (MQ=255)
                                  taataTCCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTg                                                             <  2:1977288/100‑1 (MQ=255)
                                   aataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGa                                                            <  1:1967400/100‑1 (MQ=255)
                                    atatTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGAc                                                           <  1:1040513/100‑1 (MQ=255)
                                      atTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACaa                                                         >  2:1693427/1‑100 (MQ=255)
                                       atCCCCGTCACCCTGATGGTGTCGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAAt                                                        >  2:762868/2‑100 (MQ=37)
                                        tCCCCGTCACCCTGATGGTGTCGGGGATTAAGCCCACTTTGTGCGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATa                                                       >  1:1260250/1‑100 (MQ=37)
                                         ccccGTCACCCTGATGGTGTCGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATACGTAAACTCCTGCGAAGGCGAGCCGCGTGACAATa                                                       >  1:2130980/1‑99 (MQ=25)
                                         ccccGTCACCCTGATGGTGTCGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTCAACTCCTGCGAAGGCGAGCATCGTGACAAtca                                                      >  1:1057591/1‑98 (MQ=37)
                                         ccccGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATaa                                                      >  2:802254/1‑100 (MQ=255)
                                         ccccGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATaa                                                      >  2:701639/1‑100 (MQ=255)
                                         ccccGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATaa                                                      >  2:2302408/1‑100 (MQ=255)
                                         ccccGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATaa                                                      >  2:2099249/1‑100 (MQ=255)
                                         ccccGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATaa                                                      >  2:1998938/1‑100 (MQ=255)
                                          cccGTCACCCTGATGGTGTCGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATaaa                                                     >  1:625842/1‑100 (MQ=255)
                                          cccGACACCATCAGGGTGACGGGGATTCAGCCCACTTTGTGCGCATCAGCATCCAACGACATAATTAACCTCCTGCGAAGGCGAGCTTCGTGACAATaa                                                      <  1:1156970/99‑1 (MQ=255)
                                           ccGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAAt                                                    <  2:1693641/100‑1 (MQ=255)
                                             gACACCAGCAGGGTGACTGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAGCTCCTGCGAAGGCGAGCTTCGTGACAATAAATTa                                                  <  1:2583664/100‑1 (MQ=255)
                                               cacCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTAc                                                 >  1:2228101/1‑99 (MQ=255)
                                                 ccATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtc                                              >  2:923724/1‑100 (MQ=255)
                                                 ccATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtc                                              <  2:1597339/100‑1 (MQ=255)
                                                 ccATCAGCGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtc                                              >  2:1362739/1‑100 (MQ=255)
                                                  cATCAGGGTGACGGGGATTAAGCCCACTTTGTAAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtca                                             <  1:585873/100‑1 (MQ=255)
                                                    tCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcata                                           <  2:2375595/100‑1 (MQ=255)
                                                       gggTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcatatca                                        <  1:2337583/100‑1 (MQ=255)
                                                         gTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcatatcata                                      >  1:2327160/1‑100 (MQ=255)
                                                         gTGACGGGGATTAAGCCCACTTTGTGAGCAGCAGCATGCAAAGTCATCATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcatatcata                                      <  1:1747899/100‑1 (MQ=255)
                                                           gACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACg                                    <  2:102321/100‑1 (MQ=255)
                                                             cGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGcc                                  >  1:2379664/1‑100 (MQ=255)
                                                             cGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGcc                                  >  2:2087092/1‑100 (MQ=255)
                                                             cGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGcc                                  >  2:1720214/1‑100 (MQ=255)
                                                                ggATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGc                               >  2:1333076/1‑100 (MQ=255)
                                                                    tAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTc                           >  1:457620/1‑100 (MQ=255)
                                                                    tAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTc                           <  1:1280165/100‑1 (MQ=255)
                                                                     aaGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCg                          >  1:90116/1‑100 (MQ=255)
                                                                     aaGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCg                          <  2:431431/100‑1 (MQ=255)
                                                                      aGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGa                         >  1:376224/1‑100 (MQ=255)
                                                                      aGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGa                         >  2:552021/1‑100 (MQ=255)
                                                                       gCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAATCTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGAt                        <  1:594030/100‑1 (MQ=255)
                                                                       gCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGAt                        <  2:777904/100‑1 (MQ=255)
                                                                       gCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGAt                        >  2:1122621/1‑100 (MQ=255)
                                                                       gCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGAt                        >  2:312662/1‑100 (MQ=255)
                                                                        cccACTTTGTGAGCATGAGCAACCAAAGACATAATTAAACTCCTGCGAAGGAGTGGTTCGTGCAAATAAATTACGTCATATCATACGACTGGATTCGAt                        >  1:493019/1‑99 (MQ=255)
                                                                        cccACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATa                       >  2:787035/1‑100 (MQ=255)
                                                                         ccACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAg                      <  2:1244045/100‑1 (MQ=255)
                                                                          cACTTTGTGAGCATCAGCATCCAAAGACATAATTAACCGCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGt                     <  2:1034481/100‑1 (MQ=255)
                                                                          cACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGt                     <  2:376226/100‑1 (MQ=255)
                                                                            cTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTg                   <  2:32991/100‑1 (MQ=255)
                                                                            cTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTg                   <  2:702580/100‑1 (MQ=255)
                                                                            cTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTg                   <  1:118911/100‑1 (MQ=255)
                                                                             tttcgcagcatcagcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTTACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTg                   <  1:2093794/93‑1 (MQ=255)
                                                                              ttGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGag                 >  1:702577/1‑99 (MQ=255)
                                                                               tgtgAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAat                <  2:945284/100‑1 (MQ=255)
                                                                               tgtgAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAat                <  1:313168/100‑1 (MQ=255)
                                                                               tgtgAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAat                <  2:1790563/100‑1 (MQ=255)
                                                                               tgtgAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAat                <  1:1270277/100‑1 (MQ=255)
                                                                                gtgAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAatt               >  2:610158/1‑100 (MQ=255)
                                                                                gtgAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAatt               >  2:2021326/1‑100 (MQ=255)
                                                                                      atcagcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATc         <  1:1297751/100‑1 (MQ=255)
                                                                                      atcagcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATc         <  1:1145387/100‑1 (MQ=255)
                                                                                       tcagcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCg        >  1:897584/1‑100 (MQ=255)
                                                                                       tcagcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCg        >  2:1859538/1‑100 (MQ=255)
                                                                                       tcagcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCg        <  1:210459/100‑1 (MQ=255)
                                                                                        cagcatcCAAAGACATAATTAAACTCCAGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGa       <  1:531773/100‑1 (MQ=255)
                                                                                         agcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGAc      >  1:1295665/1‑100 (MQ=255)
                                                                                         agcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGAc      >  2:269816/1‑100 (MQ=255)
                                                                                          gcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGAcc     >  1:1075822/1‑100 (MQ=255)
                                                                                           catcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCg    >  1:1600633/1‑100 (MQ=255)
                                                                                             tcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCggg  >  2:2355725/1‑100 (MQ=255)
                                                                                                |                                                                                                
AGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGG  >  NC_000913/4336928‑4337120

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