Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 531,546 C→A 5.8% E226* (GAA→TAA)  allS ← allD operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913531,5460CA 5.8% 228.7 / 2.5 69E226* (GAA→TAA) allSallD operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (28/37);  new base A (2/2);  total (30/39)
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.65e-01

AATGGCGAAGGAGGGCGCATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCC  >  NC_000913/531416‑531685
                                                                                                                                  |                                                                                                                                           
aaTGGCGAAGGAGGGCGCATCGTTGGGATTACCCGGCTGACCAGTTTTTGATTATCGATCATTGACTGGCAAATCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCggg                                                                                                                                    >  2:923029/1‑140 (MQ=255)
aaTGGCGAAGGAGGGCGCATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCggg                                                                                                                                    <  2:808417/140‑1 (MQ=255)
aaTGGCGAAGGAGGGCGCATCGTTGGGAGTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCggg                                                                                                                                    <  2:808416/140‑1 (MQ=255)
  tGGCGAAGGAGGGCGCATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGaa                                                                                                                                  >  2:1125807/1‑140 (MQ=255)
  tGGCGAAGGAGGGCGCATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGaa                                                                                                                                  >  2:652506/1‑140 (MQ=255)
       aaGGAGGGCGCATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACaata                                                                                                                             >  2:693361/1‑140 (MQ=255)
        aggaggGCTCATCGTTGGGATTACACGGCTGACCATTTGTTGATTATAGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACaataa                                                                                                                            >  1:383838/1‑140 (MQ=255)
        aggaggGCGCATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGACCcataa                                                                                                                            >  2:140243/1‑140 (MQ=255)
        aggaggGCGCATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACaataa                                                                                                                            >  2:140355/1‑140 (MQ=255)
         ggaggGCGCATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACaataat                                                                                                                           <  2:1362838/140‑1 (MQ=255)
         ggaggGCGCATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATGTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACaataat                                                                                                                           <  2:297151/140‑1 (MQ=255)
         ggaggGCGCATCGTTGGGATTACCCGGCTGACCAGTTGTAGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGAAGGGCGGCGATTTTCGTTTCCATATCGGGAACaataat                                                                                                                           <  2:1427839/140‑1 (MQ=255)
             ggCGCATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGAAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATctct                                                                                                                       >  2:11822/1‑140 (MQ=255)
               cgcATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCttt                                                                                                                     <  2:1796648/140‑1 (MQ=255)
               cgcATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCttt                                                                                                                     <  2:1796568/140‑1 (MQ=255)
                    cGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCATAGCGATTTTGGCAAAAATCCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACATTAATCTCTTTTTGcc                                                                                                                <  2:1301534/140‑1 (MQ=255)
                     gTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGccc                                                                                                               <  1:56158/140‑1 (MQ=255)
                     gTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGccc                                                                                                               <  1:56203/140‑1 (MQ=255)
                     gTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGccc                                                                                                               <  1:56241/140‑1 (MQ=255)
                       tGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTg                                                                                                             >  2:159651/1‑140 (MQ=255)
                         ggATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCTTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGc                                                                                                           <  1:1468259/140‑1 (MQ=255)
                         ggATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGc                                                                                                           <  1:312503/140‑1 (MQ=255)
                              aCCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCg                                                                                                      >  2:1719014/1‑140 (MQ=255)
                              aCCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCg                                                                                                      >  2:1574387/1‑140 (MQ=255)
                                ccGGCTGACCAGTTGTTGATTATCTATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGcc                                                                                                    >  2:1561617/1‑140 (MQ=255)
                                   gCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAgg                                                                                                 >  1:347720/1‑140 (MQ=255)
                                   gCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAgg                                                                                                 >  1:346764/1‑140 (MQ=255)
                                   gCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCAAATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCCgg                                                                                                 >  1:89209/1‑140 (MQ=255)
                                    cTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGc                                                                                                >  1:964598/1‑140 (MQ=255)
                                    cTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGc                                                                                                >  2:677498/1‑140 (MQ=255)
                                    cTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGc                                                                                                >  2:1490439/1‑140 (MQ=255)
                                    cTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGc                                                                                                >  2:1490633/1‑140 (MQ=255)
                                        ccAGTTGTTGATTATCGCTCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGAcg                                                                                            <  1:125169/140‑1 (MQ=255)
                                        ccAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGAcg                                                                                            <  1:98841/140‑1 (MQ=255)
                                            ttgttgATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGtt                                                                                        <  1:528763/140‑1 (MQ=255)
                                                 gATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGtt                                                                                   <  2:1162936/140‑1 (MQ=255)
                                                   ttATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTaa                                                                                 <  1:974079/140‑1 (MQ=255)
                                                   ttACCGATCTTTGCCTGGCTAAGCGATTTTGGCATAAAACCAATGCCAACGCCCCCCAGATGGGCGGCGATTTTCGTTTCCATTTCGGGAACAATTATCTCTTTTTGCCCGGGCAATCGCCAGGCGACGCGTTTGGTTaa                                                                                 <  2:1156682/140‑1 (MQ=255)
                                                      tCGATCATTGACTGGCAAAGCGATTTTGGCAAAAATCCAATGCCAATGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATAGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGt                                                                              >  2:1875223/1‑140 (MQ=255)
                                                      tCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCACAGGCAATCGCCAGGCGACGCGTTTGGTTAAGGt                                                                              >  2:1718038/1‑140 (MQ=255)
                                                       cGACCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTg                                                                             <  1:1513879/140‑1 (MQ=255)
                                                         aTCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTACATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCg                                                                           >  1:768071/1‑140 (MQ=255)
                                                         aTCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTACATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCg                                                                           >  1:766579/1‑140 (MQ=255)
                                                            aTTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGgc                                                                        >  2:692066/1‑140 (MQ=255)
                                                             ttGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGgcg                                                                       <  2:913732/140‑1 (MQ=255)
                                                             ttGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGgcg                                                                       <  2:396941/140‑1 (MQ=255)
                                                             atGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGgcg                                                                       <  1:160715/139‑1 (MQ=255)
                                                               gACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCt                                                                     <  1:159651/140‑1 (MQ=255)
                                                               gACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCt                                                                     <  1:726649/140‑1 (MQ=255)
                                                                    gCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCtt                                                                <  1:1414750/140‑1 (MQ=255)
                                                                           gATATTGTCAAAAAACCCATGCCAACGCCCGCCAGATGGGCGGCGATTTACGTCTCCATATCGGGAACAATAATCACTTTTAGCCCTGGCAATCGCCAGGCGACGCGTTTGGGTAAGGTGCGGGCGCTGTCTTCAATAtt                                                         <  2:1564597/140‑1 (MQ=255)
                                                                                      aaaaCCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGaa                                              <  1:1336849/140‑1 (MQ=255)
                                                                                                   aCGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTg                                 >  2:527326/1‑140 (MQ=255)
                                                                                                    cgcccgacAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCTATATTGTCCGCCGGAATGCGCCGCAAGTgc                                >  2:1183324/1‑140 (MQ=255)
                                                                                                       ccgccAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGCTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCt                             <  1:1363870/140‑1 (MQ=255)
                                                                                                        cgccAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTTAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCtt                            <  2:1408179/140‑1 (MQ=255)
                                                                                                        cgccAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCACTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCtt                            <  2:1266387/140‑1 (MQ=255)
                                                                                                          ccAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCt                          >  1:165132/1‑140 (MQ=255)
                                                                                                          ccAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCt                          >  1:1628482/1‑140 (MQ=255)
                                                                                                             gATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGtt                       <  2:1814287/140‑1 (MQ=255)
                                                                                                                gggCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGc                    >  1:623744/1‑140 (MQ=255)
                                                                                                                gggCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCGTCTGTTAGc                    >  1:526556/1‑140 (MQ=255)
                                                                                                                  gcggcgATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCgg                  <  1:838091/140‑1 (MQ=255)
                                                                                                                  gcggcgATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCgg                  <  1:191532/140‑1 (MQ=255)
                                                                                                                       gATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCtt             <  2:844568/140‑1 (MQ=255)
                                                                                                                       gATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCtt             <  2:844887/140‑1 (MQ=255)
                                                                                                                       gATTTTCGTTTACATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCtt             <  2:766579/140‑1 (MQ=255)
                                                                                                                       gATTTTCGTTTACATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCtt             <  2:768071/140‑1 (MQ=255)
                                                                                                                                 tCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGc   <  2:964598/140‑1 (MQ=255)
                                                                                                                                  ccATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGcc  <  1:43833/140‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                           
AATGGCGAAGGAGGGCGCATCGTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCC  >  NC_000913/531416‑531685

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