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,0240GA97.8% 300.1 / ‑3.0 91S3L (TCG→TTG) adiCarginine:agmatine antiporter
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (44/45);  total (44/47)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGG  >  NC_000913/4336929‑4337120
                                                                                               |                                                                                                
ggTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGa                                                                                              <  2:2330205/100‑1 (MQ=255)
      cAGGTAACAGAAATACACCTGACCCCATAATAGTCCCCGACACGATCAGGGTGACGGGGAGTAAGCCAACTTTGTGAGCATCAGCATCCAAAGACATAAt                                                                                        >  2:1858756/1‑100 (MQ=255)
      cAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAAt                                                                                        <  2:2469909/100‑1 (MQ=255)
      cAGGTAACAGAAAAACACCTGACCCCATAATATCCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACATTGTGAGCATCAGCATCCAAAGACATAAt                                                                                        >  2:1594088/1‑100 (MQ=255)
       aGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAAtt                                                                                       >  1:1877562/1‑100 (MQ=255)
        ggTAACATAAATACACCTGACCGCATAATATTCCCCGTCACCATCAGAGTGACGGGTATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTa                                                                                      <  2:867307/100‑1 (MQ=255)
        ggTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTa                                                                                      <  2:1208854/100‑1 (MQ=255)
        ggTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCAGCCAAAGACATAATTa                                                                                      >  1:2617509/1‑100 (MQ=255)
         gTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTaa                                                                                     >  1:1256876/1‑100 (MQ=255)
            aCAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACt                                                                                  <  1:1857298/100‑1 (MQ=255)
             cAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTc                                                                                 >  2:1116380/1‑100 (MQ=255)
               gAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCt                                                                               <  2:825795/100‑1 (MQ=255)
                aaaaaCACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTg                                                                              <  1:1911926/100‑1 (MQ=255)
                 aaaaCACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGc                                                                             >  2:2315072/1‑100 (MQ=255)
                  aaaCACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCg                                                                            <  1:695416/100‑1 (MQ=255)
                   aaCACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGTCCACTTTGTGATCATCAGCATCCAAAGACATAATTAAACTCCTGCGa                                                                           <  2:701149/100‑1 (MQ=255)
                    cttccgatctccATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGaa                                                                          <  1:61951/90‑1 (MQ=255)
                    acacCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGaa                                                                          >  1:1004995/1‑100 (MQ=255)
                        cTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCg                                                                      <  1:1458110/100‑1 (MQ=255)
                        cTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCg                                                                      <  1:126967/100‑1 (MQ=255)
                        cTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCg                                                                      >  2:1178028/1‑100 (MQ=255)
                         tGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGa                                                                     <  2:1169910/100‑1 (MQ=255)
                          gACCCCATAATATTCCCCGACACTATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAg                                                                    <  1:1650821/100‑1 (MQ=255)
                           aCCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGc                                                                   <  2:1079990/100‑1 (MQ=255)
                            ccccATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCt                                                                  >  1:297703/1‑100 (MQ=255)
                             cccATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCtt                                                                 >  2:767709/1‑100 (MQ=255)
                              ccATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGaaagatcggaag                                                                >  2:61951/1‑90 (MQ=255)
                               cATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCg                                                               <  2:2644649/100‑1 (MQ=255)
                                agaataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGt                                                              <  2:1607776/98‑1 (MQ=255)
                                 taataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTg                                                             <  2:1137793/100‑1 (MQ=255)
                                 taataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTg                                                             >  1:1161457/1‑100 (MQ=255)
                                 taataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTg                                                             <  1:1207937/100‑1 (MQ=255)
                                  cttaatCCCCGTCACCCTGATGGTGTCGGGGATTAAGCCCACTTTGTGAGCATCAGCATCAAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGa                                                            >  1:672038/6‑100 (MQ=25)
                                  aataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGa                                                            >  2:1210013/1‑100 (MQ=255)
                                    tatTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACa                                                          >  1:174554/1‑100 (MQ=255)
                                     atTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACaa                                                         <  2:4398/100‑1 (MQ=255)
                                      ttCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAAt                                                        >  1:2212320/1‑100 (MQ=255)
                                      ttCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGACGAGCTTCGTGACaa                                                         >  1:1543006/1‑99 (MQ=255)
                                      atCCCCGTCACCCTGATGGTGTCGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAAt                                                        >  2:1579993/2‑100 (MQ=37)
                                        ccccGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATaa                                                      >  2:184886/1‑100 (MQ=255)
                                         cccGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATaaa                                                     >  1:1237156/1‑100 (MQ=255)
                                          ccGACACCATCAGGGTGACGGGGATTAAGCCCACGTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGTTTCGTGACAATAAAt                                                    >  2:656757/1‑100 (MQ=255)
                                           cGACACCATCAGGGTGACGGGGATTAAGGCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAAt                                                    <  1:2642632/99‑1 (MQ=255)
                                           cGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAAtt                                                   <  1:2503227/100‑1 (MQ=255)
                                                ccATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtc                                              >  2:393511/1‑100 (MQ=255)
                                                 cATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtca                                             <  1:1633211/100‑1 (MQ=255)
                                                  aTCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcat                                            <  2:2617478/100‑1 (MQ=255)
                                                    cAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcatat                                          >  1:1122558/1‑100 (MQ=255)
                                                      gggTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACGCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcatat                                          <  2:1905490/98‑1 (MQ=255)
                                                      gggTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcatatca                                        <  1:726466/100‑1 (MQ=255)
                                                      gggTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcatatca                                        <  1:2064502/100‑1 (MQ=255)
                                                        gTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcatatcata                                      >  1:1211354/1‑100 (MQ=255)
                                                           aCGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGc                                   >  2:957745/1‑100 (MQ=255)
                                                            cGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGcc                                  >  1:1134405/1‑100 (MQ=255)
                                                            cGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGcc                                  >  1:2462470/1‑100 (MQ=255)
                                                            cGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGcc                                  >  1:1300599/1‑100 (MQ=255)
                                                            cGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGcc                                  >  1:431712/1‑100 (MQ=255)
                                                              gggATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTg                                <  1:2141494/100‑1 (MQ=255)
                                                                 aTTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCAt                             <  2:1953030/100‑1 (MQ=255)
                                                                   tAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTc                           <  1:2293156/100‑1 (MQ=255)
                                                                    aaGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCg                          >  2:1530040/1‑100 (MQ=255)
                                                                    aaGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCg                          >  2:1645650/1‑100 (MQ=255)
                                                                     aGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGa                         >  1:2010677/1‑100 (MQ=255)
                                                                      gCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAATTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGAt                        <  1:999585/100‑1 (MQ=255)
                                                                      gCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGAt                        <  2:2189836/100‑1 (MQ=255)
                                                                      gCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGAt                        <  2:1975770/100‑1 (MQ=255)
                                                                       cccACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATa                       >  1:2479500/1‑100 (MQ=255)
                                                                        cctcTTTGTGAGCATCAGCTTCCATAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAg                      <  1:1707047/97‑1 (MQ=255)
                                                                           cTTTGTGAGCATCAGCATCCAAAGACATAATTAGACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTg                   <  2:1256864/100‑1 (MQ=255)
                                                                           cTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTg                   <  2:2604299/100‑1 (MQ=255)
                                                                           cTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTg                   <  1:1024185/100‑1 (MQ=255)
                                                                            tttGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGa                  <  1:1712187/100‑1 (MQ=255)
                                                                             ttGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGaa                 <  1:2633331/100‑1 (MQ=255)
                                                                             ttGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGaa                 <  2:1142952/100‑1 (MQ=255)
                                                                              tgtgAGCATCAGCATCCAAAGACATAATTATACTCCTGCGAAGGTGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAat                <  1:2315101/100‑1 (MQ=255)
                                                                              tgtgAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAat                >  1:1926689/1‑100 (MQ=255)
                                                                                tgAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAatta              >  1:1196818/1‑100 (MQ=255)
                                                                                 gAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAattat             <  1:462513/100‑1 (MQ=255)
                                                                                  agcatcagcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAattatt            >  2:647804/1‑100 (MQ=255)
                                                                                   gcatcagcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAAttatta           >  2:887600/1‑100 (MQ=255)
                                                                                   gcatcagcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAAttatta           >  2:1503305/1‑100 (MQ=255)
                                                                                    catcagcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATtattat          <  2:1211342/100‑1 (MQ=255)
                                                                                      tcagcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCg        <  1:1116390/100‑1 (MQ=255)
                                                                                      tcagcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCg        >  1:509431/1‑100 (MQ=255)
                                                                                      tcagcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCg        <  2:1161447/100‑1 (MQ=255)
                                                                                         gcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGAcc     >  2:1281258/1‑100 (MQ=255)
                                                                                          catcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCg    >  1:586147/1‑100 (MQ=255)
                                                                                           atcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCgg   <  2:2075991/100‑1 (MQ=255)
                                                                                            tcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCggg  <  1:143739/100‑1 (MQ=255)
                                                                                            tcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCggg  >  1:1975242/1‑100 (MQ=255)
                                                                                            tcCAAAGACATAATTAAACTACTGCGAAGGCGTGCTTCGTAACAATAAATTACGTAATATCATACGCCCGCATTCGCTAGCTGAATTATTATCGAc      >  2:865527/1‑96 (MQ=255)
                                                                                               |                                                                                                
GGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGG  >  NC_000913/4336929‑4337120

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: