Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,337,024 G→A 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,0240GA96.6% 388.7 / 0.2 117S3L (TCG→TTG) adiCarginine:agmatine antiporter
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (55/58);  total (57/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 117 total aligned reads displayed.

AGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCG  >  NC_000913/4336928‑4337118
                                                                                                |                                                                                              
aGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAATCCCACTTTGTGAGCATCAGCATCCAAAg                                                                                             >  2:1631442/1‑100 (MQ=255)
aGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAg                                                                                             <  1:1212195/100‑1 (MQ=255)
 ggTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGa                                                                                            >  2:3196921/1‑100 (MQ=255)
    ttGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACAt                                                                                         <  1:2178369/100‑1 (MQ=255)
     tGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATa                                                                                        >  2:558243/1‑100 (MQ=255)
      gCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATaa                                                                                       >  1:1167557/1‑100 (MQ=255)
      gCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATaa                                                                                       >  1:3153420/1‑100 (MQ=255)
      gCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATaa                                                                                       >  2:668726/1‑100 (MQ=255)
       cAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAAt                                                                                      <  2:1639952/100‑1 (MQ=255)
         ggTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTa                                                                                    >  1:2282960/1‑100 (MQ=255)
          gTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTaa                                                                                   <  1:1851416/100‑1 (MQ=255)
          gTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTaa                                                                                   <  1:1572538/100‑1 (MQ=255)
           tAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTaaa                                                                                  <  2:2443436/100‑1 (MQ=255)
            aaCAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAAc                                                                                 <  2:1244796/100‑1 (MQ=255)
            aaCAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAAc                                                                                 <  2:1120033/100‑1 (MQ=255)
             aCAGAAAAACACCGGACCCCATAATATTCCCCGACACGATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTaacc                                                                                 <  2:285670/99‑3 (MQ=255)
              cAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTc                                                                               >  2:1608549/1‑100 (MQ=255)
              cAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTc                                                                               <  1:1631446/100‑1 (MQ=255)
              cAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTc                                                                               >  2:642554/1‑100 (MQ=255)
                gAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCt                                                                             >  1:3340617/1‑100 (MQ=255)
                   aaaCACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCg                                                                          <  2:1602133/100‑1 (MQ=255)
                   aaaCACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCg                                                                          <  2:1329533/100‑1 (MQ=255)
                      cacCTGACCCCATAATATTCCCCGACACCCTCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATTAGCATCCAAAGACATAATTAAACTCCTGCGAAg                                                                       <  2:2022068/100‑1 (MQ=255)
                        ccTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGc                                                                     >  2:3218182/1‑100 (MQ=255)
                        ccTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGc                                                                     >  2:2854952/1‑100 (MQ=255)
                        ccTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGc                                                                     <  1:2577022/100‑1 (MQ=255)
                           gACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAg                                                                  >  1:2600636/1‑100 (MQ=255)
                            aCCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGc                                                                 <  2:2097174/100‑1 (MQ=255)
                             ccccATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCt                                                                <  2:2294207/100‑1 (MQ=255)
                             ccccATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCt                                                                >  1:3074254/1‑100 (MQ=255)
                             ccccATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCt                                                                >  1:3074222/1‑100 (MQ=255)
                               ccATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTc                                                              >  1:63755/1‑100 (MQ=255)
                                 ataataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGt                                                            <  1:1564800/100‑1 (MQ=255)
                                 ataataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGt                                                            >  1:130820/1‑100 (MQ=255)
                                  taataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTg                                                           <  1:2453644/100‑1 (MQ=255)
                                     tatTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACa                                                        >  1:1938142/1‑100 (MQ=255)
                                       ttCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAAt                                                      <  1:642555/100‑1 (MQ=255)
                                        tCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATa                                                     <  1:2053401/100‑1 (MQ=255)
                                        tCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATa                                                     >  1:391722/1‑100 (MQ=255)
                                        tCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATa                                                     >  1:1500599/1‑100 (MQ=255)
                                          cccGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATaaa                                                   >  1:1287207/1‑100 (MQ=255)
                                            cGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAAtt                                                 <  1:1958726/100‑1 (MQ=255)
                                             gACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTa                                                <  1:2303735/100‑1 (MQ=255)
                                               cacCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACg                                              <  2:3340598/100‑1 (MQ=255)
                                               cacCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACg                                              >  2:1733436/1‑100 (MQ=255)
                                                 ccATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtc                                            <  1:98080/100‑1 (MQ=255)
                                                 ccATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtc                                            >  1:1342465/1‑100 (MQ=255)
                                                 ccATCAGGGTGACGGAGATAAAGCCCACTAAGAGAACATCAACATCCGTAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGt                                             <  1:283535/99‑1 (MQ=255)
                                                  cATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtca                                           <  2:2572584/100‑1 (MQ=255)
                                                  cATCAGGGCGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtca                                           >  1:1074118/1‑100 (MQ=255)
                                                   aTCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcat                                          <  1:1347702/100‑1 (MQ=255)
                                                   aTCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcat                                          <  1:668727/100‑1 (MQ=255)
                                                      aGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcatatc                                       >  1:1921573/1‑100 (MQ=255)
                                                       gggTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcatatcn                                      <  2:2949391/100‑2 (MQ=255)
                                                       gggTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcatatca                                      <  2:3204517/100‑1 (MQ=255)
                                                       gggTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcatatca                                      <  2:1388944/100‑1 (MQ=255)
                                                       gggTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcatatca                                      >  1:2486295/1‑100 (MQ=255)
                                                       gggTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcatatca                                      >  2:1057765/1‑100 (MQ=255)
                                                        ggTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTACACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcatatcat                                     <  1:290593/100‑1 (MQ=255)
                                                         gTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcatatcata                                    >  1:2877965/1‑100 (MQ=255)
                                                          tGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATAc                                   >  1:1077314/1‑100 (MQ=255)
                                                           gACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATAc                                   >  1:1896138/1‑99 (MQ=255)
                                                             cGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGCGCTTCGTGACAATAAATTACGTCATATCATAcccc                                >  1:952884/1‑100 (MQ=255)
                                                             cGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGcc                                >  1:2232379/1‑100 (MQ=255)
                                                             cGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGcc                                >  1:1402019/1‑100 (MQ=255)
                                                             cGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGcc                                <  1:889611/100‑1 (MQ=255)
                                                                ggATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGc                             >  2:1407406/1‑100 (MQ=255)
                                                                ggATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGc                             <  1:766266/100‑1 (MQ=255)
                                                                 gATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCa                            <  2:2877950/100‑1 (MQ=255)
                                                                 gATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCa                            <  1:3058398/100‑1 (MQ=255)
                                                                  aTTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCAt                           <  1:3196939/100‑1 (MQ=255)
                                                                    tAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTc                         <  2:1724185/100‑1 (MQ=255)
                                                                     aaGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCg                        <  1:2304053/100‑1 (MQ=255)
                                                                      aGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGa                       >  1:3351871/1‑100 (MQ=255)
                                                                      aGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGa                       <  2:2967354/100‑1 (MQ=255)
                                                                      aGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGa                       <  2:1495463/100‑1 (MQ=255)
                                                                       gCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGAt                      <  1:666769/100‑1 (MQ=255)
                                                                        cccACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATa                     >  2:1858302/1‑100 (MQ=255)
                                                                         ccACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAg                    <  2:431476/100‑1 (MQ=255)
                                                                         ccACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAg                    >  1:574663/1‑100 (MQ=255)
                                                                         ccACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAg                    >  1:2199407/1‑100 (MQ=255)
                                                                          cACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGt                   >  2:1329399/1‑100 (MQ=255)
                                                                          cACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGt                   >  1:1347833/1‑100 (MQ=255)
                                                                          cACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGt                   >  2:2770369/1‑100 (MQ=255)
                                                                            cTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGttn                 <  2:778558/100‑2 (MQ=255)
                                                                            cTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTg                 <  2:1342461/100‑1 (MQ=255)
                                                                            cTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTg                 <  2:1167554/100‑1 (MQ=255)
                                                                              ttGTGAGCATCAGCGTCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGaa               <  2:2282949/100‑1 (MQ=255)
                                                                                gtgAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAatt             >  1:3261256/1‑100 (MQ=255)
                                                                                gtgAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAatt             <  1:3267356/100‑1 (MQ=255)
                                                                                gtgAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAatt             <  1:892437/100‑1 (MQ=255)
                                                                                    gcatcagcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAAttatta         <  2:3074205/100‑1 (MQ=255)
                                                                                    gcatcagcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAAttatta         <  2:3074237/100‑1 (MQ=255)
                                                                                     catcagcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATtattat        >  1:2514368/1‑100 (MQ=255)
                                                                                     catcagcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATtattat        >  1:1394807/1‑100 (MQ=255)
                                                                                      atcagcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATc       >  1:1789457/1‑100 (MQ=255)
                                                                                      atcagcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACACCTGCATTCGATAGTTGAATTATTATc       <  2:2517539/100‑1 (MQ=255)
                                                                                       tcagcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCg      >  2:2498792/1‑100 (MQ=255)
                                                                                           catcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCg  >  1:1800639/1‑100 (MQ=255)
                                                                                           catcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCg  <  1:1524962/100‑1 (MQ=255)
                                                                                                |                                                                                              
AGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCG  >  NC_000913/4336928‑4337118

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: