Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,358,043 C→A 7.5% *268L (TGA→TTA)  rhmA ← 2‑keto‑3‑deoxy‑L‑rhamnonate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,358,0430CA 7.5% 270.8 / 7.4 80*268L (TGA→TTA) rhmA2‑keto‑3‑deoxy‑L‑rhamnonate aldolase
Reads supporting (aligned to +/- strand):  ref base C (36/38);  new base A (3/3);  total (39/41)
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.78e-01

TTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACGCCAGGCATTGCTGCGCCA  >  NC_000913/2357907‑2358181
                                                                                                                                        |                                                                                                                                          
ttGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAAt                                                                                                                                         >  1:695572/1‑140 (MQ=255)
 tGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATa                                                                                                                                        <  1:1020390/140‑1 (MQ=255)
 tGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATa                                                                                                                                        <  1:1020389/140‑1 (MQ=255)
   cTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAAc                                                                                                                                      <  2:302880/140‑1 (MQ=255)
        cgcATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTTGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCAAGGGCCTTTGATATCAATAACAACCt                                                                                                                                 >  2:641181/1‑140 (MQ=255)
        cgcATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCt                                                                                                                                 >  2:375739/1‑140 (MQ=255)
                    cAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGc                                                                                                                     <  2:983639/140‑1 (MQ=255)
                    cAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGc                                                                                                                     <  2:1149155/140‑1 (MQ=255)
                       ggCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCa                                                                                                                  >  2:799213/1‑140 (MQ=255)
                           ttGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTAGATATCAATAACTACCTTTTATGCGTGGCCCAtttt                                                                                                              >  2:198771/1‑140 (MQ=255)
                              aTCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACGGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACATTTTATGCGTGGCCCATTTTTGc                                                                                                           <  2:1127484/140‑1 (MQ=255)
                              aTCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGc                                                                                                           <  2:639057/140‑1 (MQ=255)
                              aTCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGc                                                                                                           <  2:907480/140‑1 (MQ=255)
                                         gTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAAc                                                                                                >  1:121053/1‑140 (MQ=255)
                                              cgcgGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGc                                                                                           >  2:891728/1‑140 (MQ=255)
                                                    cACCCAGCTTAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACTCAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCg                                                                                     <  2:1183625/140‑1 (MQ=255)
                                                     aCCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGt                                                                                    <  1:23808/140‑1 (MQ=255)
                                                         aGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGAt                                                                                <  2:1192954/140‑1 (MQ=255)
                                                         aGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGAt                                                                                <  2:1192922/140‑1 (MQ=255)
                                                         aGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGAt                                                                                <  2:1192952/140‑1 (MQ=255)
                                                           cTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATcc                                                                              >  1:1053373/1‑140 (MQ=255)
                                                              aGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCTAgg                                                                           >  1:224080/1‑140 (MQ=255)
                                                               gCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAggg                                                                          <  1:653276/140‑1 (MQ=255)
                                                                  caGCCGAACTGGCAGCGAATAACCACAAAAAATCACATGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCAATTATGCGTGGCCCAATTTTGCCTGATTAAAACATCGCCAGTCGTTGATCCAGGGCAt                                                                       <  2:1087344/139‑1 (MQ=255)
                                                                  aaGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATACCTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCAt                                                                       <  2:945296/140‑1 (MQ=255)
                                                                  aaGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCAt                                                                       <  1:901433/140‑1 (MQ=255)
                                                                  aaGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCAt                                                                       <  1:11298/140‑1 (MQ=255)
                                                                  aaGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCAt                                                                       <  2:1120815/140‑1 (MQ=255)
                                                                  aaGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCAt                                                                       <  1:1044868/140‑1 (MQ=255)
                                                                  aaGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCAt                                                                       <  1:1076080/140‑1 (MQ=255)
                                                                  aaGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCAt                                                                       <  1:194080/140‑1 (MQ=255)
                                                                  aaGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCAt                                                                       <  2:945202/140‑1 (MQ=255)
                                                                  aaGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACAACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCAt                                                                       <  2:1087489/140‑1 (MQ=255)
                                                                  aaGCCGATCAGGCAGCGAACCACCACTAAATCCCACACGCCCACAGCCGATCCCCCTGGGCCTTTGATCTCAATAACTACCTATTATGCGTGGCCCATTTTTGCCTGATTTAAACATCTCCAGACGTTGATCCAGGGCAt                                                                       <  1:1040136/140‑1 (MQ=255)
                                                                    gCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCg                                                                     <  1:446881/140‑1 (MQ=255)
                                                                    gCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCg                                                                     <  1:446801/140‑1 (MQ=255)
                                                                    gCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCg                                                                     <  1:1214861/140‑1 (MQ=255)
                                                                     acGAACTGGCAGCGAATAACCACTAATTAACAAAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGc                                                                    >  1:223091/2‑140 (MQ=255)
                                                                      cGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCt                                                                   <  2:121053/140‑1 (MQ=255)
                                                                        aTCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGt                                                                 >  1:276676/1‑140 (MQ=255)
                                                                        aTCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGt                                                                 >  1:1148042/1‑140 (MQ=255)
                                                                               aGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCAt                                                          <  2:28743/140‑1 (MQ=255)
                                                                               aGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCAt                                                          >  1:1202600/1‑140 (MQ=255)
                                                                                gCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATc                                                         <  2:185067/140‑1 (MQ=255)
                                                                                           cTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCGAGGGCATCGCTGTAGAGGATCGTGTCAACGcc                                              >  2:688387/1‑140 (MQ=255)
                                                                                           cTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGcc                                              >  2:563296/1‑140 (MQ=255)
                                                                                           cTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGcc                                              >  2:192878/1‑140 (MQ=255)
                                                                                             aaTTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCTTCGTGTCAACGCCaa                                            <  2:511843/140‑1 (MQ=255)
                                                                                                   cacaGGCCCACAGCCGATCCCCATGGGCCTTTGATATAAATAACTa                                                                                                                                    <  2:691833/46‑1 (MQ=255)
                                                                                                   cacaGGCCCACAGCCGATCCCCATGGGCCTTTGATATAAATAACTa                                                                                                                                    <  2:691823/46‑1 (MQ=255)
                                                                                                   cacaGGCCCACAGCCGATCCCCATGGGCCTTTGATATAAATAACTa                                                                                                                                    <  2:691795/46‑1 (MQ=255)
                                                                                                   cacaGGCCCACAGCCGATCCCCATGGGCCTTTGATATAAATAACTa                                                                                                                                    >  1:691795/1‑46 (MQ=255)
                                                                                                   cacaGGCCCACAGCCGATCCCCATGGGCCTTTGATATAAATAACTa                                                                                                                                    >  1:691823/1‑46 (MQ=255)
                                                                                                   cacaGGCCCACAGCCGATCCCCATGGGCCTTTGATATAAATAACTa                                                                                                                                    >  1:691833/1‑46 (MQ=255)
                                                                                                        gCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCa                                                                                         >  2:1131702/1‑84 (MQ=255)
                                                                                                        gCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCa                                                                                         <  1:1131702/84‑1 (MQ=255)
                                                                                                        gCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAg                                 >  1:166301/1‑140 (MQ=255)
                                                                                                        gCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAg                                 >  1:392419/1‑140 (MQ=255)
                                                                                                             caGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGc                            >  1:946902/1‑140 (MQ=255)
                                                                                                             caGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGc                            >  1:81692/1‑140 (MQ=255)
                                                                                                             caGCCGATCCCCATGGGCCTGTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGAAGAGCATCGTGTCAACGCCAACAGCGACAAAGTACGc                            >  2:867428/1‑140 (MQ=255)
                                                                                                                 cGATCCCCATGGGCCTATGATATCAATAACTACCTATTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTATAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTccc                        <  1:160447/140‑1 (MQ=255)
                                                                                                                  gaaccccATGGGCCTTTGATAACAATAACTACGTTTTATGCGTGGCCCATTTTAGCCTGATTAATACATCGCCAGTCGTTGATCCAGGGCATCGCTGAAGAGCATCGAGTCAACGCCAACAGCGACAAAGTTCGCTcccc                       >  2:527813/4‑140 (MQ=255)
                                                                                                                  gATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGTTTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTcccc                       >  2:1120362/1‑140 (MQ=255)
                                                                                                                  gATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGACAACAGCGACAAAGTTCGCTcccc                       >  2:723940/1‑140 (MQ=255)
                                                                                                                   aTCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCa                      >  1:681252/1‑140 (MQ=255)
                                                                                                                   aTCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCa                      >  1:557308/1‑140 (MQ=255)
                                                                                                                   aTCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCa                      >  2:1078067/1‑140 (MQ=255)
                                                                                                                     ccccATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACg                    <  2:392419/140‑1 (MQ=255)
                                                                                                                     ccccAAGGGCCATTGAAATCAAAAACTCCCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACg                    <  2:1205668/140‑1 (MQ=255)
                                                                                                                      cccATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACGc                   >  1:778157/1‑140 (MQ=255)
                                                                                                                          tGGGCCTTTGATCTCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACGCCAgg               >  2:461498/1‑140 (MQ=255)
                                                                                                                          tGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACGCCAgg               >  1:621753/1‑140 (MQ=255)
                                                                                                                            ggCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACGCCAGGCa             >  1:652216/1‑140 (MQ=255)
                                                                                                                                tttGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACGCCAGGCATtgc         >  2:816010/1‑140 (MQ=255)
                                                                                                                                tttGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACGCCAGGCATtgc         >  2:920141/1‑140 (MQ=255)
                                                                                                                                  tGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACGCCAGGCATtgctg       <  1:93643/140‑1 (MQ=255)
                                                                                                                                       tCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACGCCAGGCATTGCTGCGCCa  <  2:404084/140‑1 (MQ=255)
                                                                                                                                       tCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACGCCAGGCATTGCTGCGCCa  >  1:356775/1‑140 (MQ=255)
                                                                                                                                       tCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACGCCAGGCATTGCTGCGCCa  >  1:180972/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                          
TTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTGGCAGCGAATAACCACTAATTAACACAGGCCCACAGCCGATCCCCATGGGCCTTTGATATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACGCCAGGCATTGCTGCGCCA  >  NC_000913/2357907‑2358181

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: