Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,490,141 C→T *333* (TAG→TAA idnR ← transcriptional repressor, 5‑gluconate‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,490,1410CT98.7% 244.9 / NA 76*333* (TAG→TAAidnRtranscriptional repressor, 5‑gluconate‑binding
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (39/36);  minor base G (1/0);  total (40/36)
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

CGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGTTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATCAAAACGGGGAGTAATTACGCTGGCG  >  NC_000913/4490007‑4490276
                                                                                                                                      |                                                                                                                                       
ttagccattcccggcagtaaaaacagttctgtgtccggcgtgcgggcaattaacaggggttcactcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTaaa                                                                                                                                     >  2:521417‑M2/121‑139 (MQ=255)
ttagccattcccggcagtaaaaacagttctgtgtccggcgtgcgggcaattaacaggggttcactcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTa                                                                                                                                       >  1:139938‑M2/121‑137 (MQ=255)
     cattcccggcagtaaaaacagttctgtgtccggcgtgcgggcaattaacaggggttcactcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCgtg                                                                                                                                <  2:386479‑M2/24‑1 (MQ=255)
         cccggcagtaaaaacagttctgtgtccggcgtgcgggcaattaacaggggttcactcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGc                                                                                                                            <  1:1613576‑M2/28‑1 (MQ=255)
             gcagtaaaaacagttctgtgtccggcgtgcgggcaattaacaggggttcactcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTg                                                                                                                        <  1:215037‑M2/32‑1 (MQ=255)
                           tctgtgtccggcgtgcgggcaattaacaggggttcactcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATAcc                                                                                                          >  1:668144‑M2/94‑139 (MQ=255)
                           tctgtgtccggcgtgcgggcaattaacaggggttcactcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATAcc                                                                                                          >  1:381605‑M2/94‑139 (MQ=255)
                            ctgtgtccggcgtgcgggcaattaacaggggttcactcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATAccc                                                                                                         <  1:521417‑M2/47‑1 (MQ=255)
                               agtccgccgtcagggcaattaacaggggttcactcatgagctatcgtccagtttcctgcctggagtataggcaacgtccggcaaagagaCAAGCGAAACTAACTTAAAGTGTGTTGCCGTGAGAAATCTGTTACCCAAg                                                                                                      <  1:1613245‑M2/50‑1 (MQ=255)
                                    ggcgtgcgggcaattaacaggggttcactcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGt                                                                                                    <  2:1013780‑M2/52‑1 (MQ=255)
                                    ggcgtgcgggcaattaacaggggttcactcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACGTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGt                                                                                                    >  1:1013780‑M2/85‑136 (MQ=255)
                                      cgtgcgggcaattaacaggggttcactcatgagatatcctccagtttcctgcctggagtataggcgacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGcc                                                                                                                           <  2:231351‑M2/29‑1 (MQ=255)
                                      cgtgcgggcaattaacaggggttcactcatgagatatcctccagtttcctgcctggagtataggcgacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGcc                                                                                                                           >  1:231351‑M2/83‑111 (MQ=255)
                                           gggcaattaacaggggttcactcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGtgtt                                                                                          <  2:381605‑M2/62‑1 (MQ=255)
                                                    acaggggttcactcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTgttgtggttgtg                                                                                 <  2:715782‑M2/71‑1 (MQ=255)
                                                       ggggttcactcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATc                                                                              >  1:835671‑M2/66‑139 (MQ=255)
                                                           ttcactcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGc                                                                                                                            <  1:481144‑M2/28‑1 (MQ=255)
                                                           ttcactcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGc                                                                                                                            >  2:481144‑M2/62‑89 (MQ=255)
                                                             cactcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTAtt                                                                        >  2:527573‑M2/60‑139 (MQ=255)
                                                                  atgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTaa                                                                   >  2:860232‑M2/55‑139 (MQ=255)
                                                                          tcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGt                                                                                                    <  2:1453900‑M2/52‑1 (MQ=255)
                                                                          tcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGt                                                                                                    >  1:1453900‑M2/47‑98 (MQ=255)
                                                                          tcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTgttgtggtt                                                                                    <  1:874218‑M2/68‑1 (MQ=255)
                                                                          tcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTgttgtggtt                                                                                    >  2:874218‑M2/47‑114 (MQ=255)
                                                                          tcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCt                                                            >  1:635122‑M2/47‑138 (MQ=255)
                                                                          tcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCt                                                            <  2:635122‑M2/92‑1 (MQ=255)
                                                                            ctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAg                                                         >  2:550933‑M2/45‑139 (MQ=255)
                                                                            ctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAg                                                         <  2:835671‑M2/95‑1 (MQ=255)
                                                                             tccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGc                                                        >  2:405086‑M2/44‑139 (MQ=255)
                                                                              ccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCa                                                       >  1:502164‑M2/43‑139 (MQ=255)
                                                                               cagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCaa                                                      >  2:880064‑M2/42‑139 (MQ=255)
                                                                                       tgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTggg                                              >  2:526050‑M2/34‑139 (MQ=255)
                                                                                              agtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCAt                                       >  2:1499695‑M2/27‑139 (MQ=255)
                                                                                              agtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCAt                                       <  1:550933‑M2/113‑1 (MQ=255)
                                                                                                   aggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCaa                                                                                                 >  1:1211781‑M2/22‑76 (MQ=255)
                                                                                                   aggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCaa                                                                                                 <  2:1211781‑M2/55‑1 (MQ=255)
                                                                                                   aggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAg                                                         <  2:1555811‑M2/95‑1 (MQ=255)
                                                                                                   aggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAg                                                         >  1:1555811‑M2/22‑116 (MQ=255)
                                                                                                   aggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGAc                                  >  1:743534‑M2/22‑139 (MQ=255)
                                                                                                    ggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGAcc                                 <  1:1499695‑M2/119‑1 (MQ=255)
                                                                                                        acgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGtt                                                                           <  2:237714‑M2/77‑1 (MQ=255)
                                                                                                        acgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGtt                                                                           >  1:237714‑M2/17‑93 (MQ=255)
                                                                                                        acgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTg                                                              >  2:514486‑M2/17‑106 (MQ=255)
                                                                                                        acgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTg                                                              <  1:514486‑M2/90‑1 (MQ=255)
                                                                                                        acgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGAc                                  >  1:806484‑M2/17‑134 (MQ=255)
                                                                                                        acgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGAc                                  <  2:806484‑M2/118‑1 (MQ=255)
                                                                                                        acgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAata                             >  2:290291‑M2/17‑139 (MQ=255)
                                                                                                        acgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAata                             >  2:1382171‑M2/17‑139 (MQ=255)
                                                                                                        acgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAata                             >  2:1092276‑M2/17‑139 (MQ=255)
                                                                                                        acgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAata                             >  1:597977‑M2/17‑139 (MQ=255)
                                                                                                         cgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAgtgt                                                                                           <  2:813867‑M2/61‑1 (MQ=255)
                                                                                                         cgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAgtgt                                                                                           >  1:813867‑M2/16‑76 (MQ=255)
                                                                                                         cgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGc                                                             >  1:1091733‑M2/16‑106 (MQ=255)
                                                                                                         cgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGc                                                             <  2:1091733‑M2/91‑1 (MQ=255)
                                                                                                          gtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGggcg                                            >  2:371266‑M2/15‑122 (MQ=255)
                                                                                                          gtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGggcg                                            <  1:371266‑M2/108‑1 (MQ=255)
                                                                                                           tccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATCa                          <  1:527573‑M2/126‑1 (MQ=255)
                                                                                                            ccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGt                                                                                                    <  1:1153255‑M2/52‑1 (MQ=255)
                                                                                                            ccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGt                                                                                                    >  2:1153255‑M2/13‑64 (MQ=255)
                                                                                                                caaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGTAACATTTGGGCGGCCATCCGACCAATATCAAAACg                     <  1:860232‑M2/131‑1 (MQ=255)
                                                                                                                  aagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATCAAAACggg                   <  2:743534‑M2/133‑1 (MQ=255)
                                                                                                                   agagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATCAAAACgggg                  >  2:947005‑M2/6‑139 (MQ=255)
                                                                                                                   agagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATCAAAACgggg                  <  2:1462879‑M2/134‑1 (MQ=255)
                                                                                                                   agagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATCAAAACgggg                  <  2:597977‑M2/134‑1 (MQ=255)
                                                                                                                    gagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATCAAAACGGGGa                 <  1:1344646‑M2/135‑1 (MQ=255)
                                                                                                                     agaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATCAAAACGGGGAg                >  1:1369100‑M2/4‑139 (MQ=255)
                                                                                                                         aaGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATCAAAACGGGGAGTAAt            >  1:1555052/1‑139 (MQ=255)
                                                                                                                              aaaaTAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATCAAAACg                     <  2:581914/125‑1 (MQ=255)
                                                                                                                              aaaaTAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATCAAAACg                     >  1:581914/1‑125 (MQ=255)
                                                                                                                               aaaTAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATCAAAACg                     <  1:1527248/124‑1 (MQ=255)
                                                                                                                               aaaTAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATCAAAACg                     >  2:1527248/1‑124 (MQ=255)
                                                                                                                                aaTAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCaa                                                                                                 <  1:436210/47‑1 (MQ=255)
                                                                                                                                aaTAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCaa                                                                                                 >  2:436210/1‑47 (MQ=255)
                                                                                                                                 aTAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATCAAAACGGGGAGTAATTACGCTgg    >  1:260601/1‑139 (MQ=255)
                                                                                                                                  tAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATCAAAACGGGGAGTAATTACGCTGGc   <  2:1369100/139‑1 (MQ=255)
                                                                                                                                   aaCTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATCAAAACGGGGAGTAATTACGCTGGCg  <  2:260601/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                       
CGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGTTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATCAAAACGGGGAGTAATTACGCTGGCG  >  NC_000913/4490007‑4490276

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: