Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,066,293 T→C 5.0% D294G (GAC→GGC)  ygfI ← putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,066,2930TC 5.0% 732.4 / 16.4 218D294G (GAC→GGC) ygfIputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (103/104);  new base C (2/9);  total (105/113)
Fisher's exact test for biased strand distribution p-value = 6.07e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.97e-01

Only 100 of 227 total aligned reads displayed.

CGTGCCAATGTGCAGCCATGTTGATCACAGCGAACACAGCGTCAAAGTGATCATCACCGAACAAGGGATCGCCGATCTGCGCGGTCTTTCCCCGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCACCGTCAGAAAATGTCCGCGCTCCACGGCGCTGATTTTGGCCATAAAACTGTAGGGACGATAAATAATGGTTGGATTGATCCCACAAGGTCGCATATTCGCATCAATCATTGCTTC  >  NC_000913/3066022‑3066561
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                             
 gtGCCAATGTGCAGCCATGTTGATCACAGCGAACACAGCGTCAAAGTGATCATCACCGAACAAGGGATCGCCGATCTGCGCGGTCTTTCCCCGCTTCGACGCGCCCGCACTATCATTGATAATTTTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTCCACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACg                                                                                                                                                                                                                                                                           >  4:176838/1‑275 (MQ=255)
      aaTGTGCCGCCAGGTTGATCACAGCGAACACATCGTCAAACTGATCATCACCGAACAAGGGATCGCCGTTCTGCGCGGTCTTTCCCCGGTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCtt                                                                                                                                                                                                                                                                     <  4:70804/275‑1 (MQ=255)
                                  cacaGCGTCAAAGTGATCATCACCGAACAAGGGATCGCCGATCTGCGCGGTCTTTCCCCGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCtt                                                                                                                                                                                                                                                                     <  2:279046/247‑1 (MQ=255)
                                                           aaCAAGGGATCGCCGATCTGCGCGGTCTTTCCCCGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATaa                                                                                                                                                                                                                >  2:339523/1‑275 (MQ=255)
                                                                                         ccccGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGt                                                                                                                                                                                                                                                                               >  3:54895/1‑183 (MQ=255)
                                                                                                                           ttGTGCACATCTTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAGTTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGt                                                                                                                                                                                                                                                                               >  7:961065/1‑149 (MQ=255)
                                                                                                                           ttGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGt                                                                                                                                                                                                                                                                               >  5:317394/1‑149 (MQ=255)
                                                                                                                           ttGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGt                                                                                                                                                                                                                                                                               >  5:1377297/1‑149 (MQ=255)
                                                                                                                              tgCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCAc                                                                                                                                                                                                                                                                            <  6:574130/149‑1 (MQ=255)
                                                                                                                              tgCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGAtt                                                                                                                                                                                                                      >  2:85785/1‑202 (MQ=255)
                                                                                                                                cacaTCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGAt                                                                                                                                                                                                                       >  4:273485/1‑199 (MQ=255)
                                                                                                                                 acaTCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAt                                                                                                                                                                                                                                                                         >  5:717832/1‑149 (MQ=255)
                                                                                                                                  caTCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAAt                                                                                                                                                                                                                  <  4:85785/202‑1 (MQ=255)
                                                                                                                                    tCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAt                                                                                                                                                                                                                                                                         <  6:409630/146‑1 (MQ=255)
                                                                                                                                     ccTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTc                                                                                                                                                                                                                                                                    >  6:343813/1‑149 (MQ=255)
                                                                                                                                      cTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCa                                                                                                                                                                                                                                                                   <  6:1005346/149‑1 (MQ=255)
                                                                                                                                        aTGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGc                                                                                                                                                                                                                                                                 <  7:1073277/149‑1 (MQ=255)
                                                                                                                                        aTGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTcc                                                                                                                                                                                                       <  1:68319/207‑1 (MQ=255)
                                                                                                                                            aTCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAata                                                                                                                                                                                                                                                             >  6:560037/1‑149 (MQ=255)
                                                                                                                                             tCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAatat                                                                                                                                                                                                                                                            >  6:677088/1‑149 (MQ=255)
                                                                                                                                              cGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATg                                                                                                                                                                                                                                                           <  7:568405/149‑1 (MQ=255)
                                                                                                                                               gggATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATTTTCAGCAATCTGCTGAAGCATCGTGGAGGTATAACGATGATTTCTGATTGAATAAAGATATTACTGCATATACCTCGGCCTGTGCGGTTCCTCAAAATATTTCAGAGTTTCCGGATTCAACAGTAGCCAGGCACAAttc                                                                                                                            >  1:40669/1‑273 (MQ=255)
                                                                                                                                                ggATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCatattgct                                                                                                                                                                                                                                                         >  8:540182/1‑142 (MQ=255)
                                                                                                                                                ggATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCt                                                                                                                                                                                                                                                         >  6:1284195/1‑149 (MQ=255)
                                                                                                                                                   ttATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGaa                                                                                                                                                                                                                                                      <  5:1284195/149‑1 (MQ=255)
                                                                                                                                                   ttATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGaa                                                                                                                                                                                                                                                      <  7:560266/149‑1 (MQ=255)
                                                                                                                                                    tATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCa                                                                                                                                  <  2:176982/264‑1 (MQ=255)
                                                                                                                                                     aTCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGc                                                                                                                                                                                                                                                    <  8:408461/149‑1 (MQ=255)
                                                                                                                                                       cTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCAt                                                                                                                                                                                                                                                  <  5:1168822/149‑1 (MQ=255)
                                                                                                                                                        tGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGCATATTGCTGAAgatc                                                                                                                                                                                                                                                  <  8:1236593/149‑4 (MQ=255)
                                                                                                                                                           aTCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCt                                                                                                                                                                                                                                                                      >  5:41800/1‑125 (MQ=255)
                                                                                                                                                           aTCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTCAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCt                                                                                                                                                                                                                                                                      <  6:41800/125‑1 (MQ=255)
                                                                                                                                                            tCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTgg                                                                                                                                                                                                                                             >  8:836389/1‑149 (MQ=255)
                                                                                                                                                            tCGCTATCTGGAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGaa                                                                                                                                                                                                                                                      >  8:767770/1‑139 (MQ=255)
                                                                                                                                                                  tCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGCATATTGCTGAAGATCGTgac                                                                                                                                                                                                                                            <  5:590847/144‑3 (MQ=38)
                                                                                                                                                                     ggAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAAc                                                                                                                                                                                                                                    >  7:1428257/1‑149 (MQ=255)
                                                                                                                                                                      gAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCa                                                                                                                                                                                                                                                                   >  5:393006/1‑117 (MQ=255)
                                                                                                                                                                        aaaTGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACgat                                                                                                                                                                                                                                 <  8:370014/149‑1 (MQ=255)
                                                                                                                                                                         aaTGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACgatg                                                                                                                                                                                                                                >  7:87051/1‑149 (MQ=255)
                                                                                                                                                                           tGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGc                                                                                                                                                                                                                                                    >  7:157118/1‑127 (MQ=255)
                                                                                                                                                                           tGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGc                                                                                                                                                                                                                                                    <  8:157118/127‑1 (MQ=255)
                                                                                                                                                                           tGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACgatgat                                                                                                                                                                                                                              <  8:495107/149‑1 (MQ=255)
                                                                                                                                                                             cgcCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGAttt                                                                                                                                                                                                                            <  7:502010/149‑1 (MQ=255)
                                                                                                                                                                              gcCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTc                                                                                                                                                                                                                           <  6:992738/149‑1 (MQ=255)
                                                                                                                                                                              gcCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTc                                                                                                                                                                                                                           <  7:336565/149‑1 (MQ=255)
                                                                                                                                                                               ccTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCt                                                                                                                                                                                                                          <  5:1211299/149‑1 (MQ=255)
                                                                                                                                                                                  ggCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGAt                                                                                                                                                                                                                       <  8:250471/149‑1 (MQ=255)
                                                                                                                                                                                  ggCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGAt                                                                                                                                                                                                                       <  8:824079/149‑1 (MQ=255)
                                                                                                                                                                                          tatTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAATCATCGTGGCGGTATAACGATGATTTCTGATTGAATaaa                                                                                                                                                                                                               <  8:1201969/149‑1 (MQ=255)
                                                                                                                                                                                           atTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAat                                                                                                                                                                                                                                                              <  8:595994/101‑1 (MQ=255)
                                                                                                                                                                                           atTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAat                                                                                                                                                                                                                                                              >  7:595994/1‑101 (MQ=255)
                                                                                                                                                                                           atTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAg                                                                                                                                                                                                              >  5:471985/1‑149 (MQ=255)
                                                                                                                                                                                            ttCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGa                                                                                                                                                                                                             >  7:914529/1‑149 (MQ=255)
                                                                                                                                                                                              caccacGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGata                                                                                                                                                                                                           >  8:1439724/1‑149 (MQ=255)
                                                                                                                                                                                               accacGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGatat                                                                                                                                                                                                          >  5:804431/1‑149 (MQ=255)
                                                                                                                                                                                               accacGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGatat                                                                                                                                                                                                          >  7:688143/1‑149 (MQ=255)
                                                                                                                                                                                                 cacGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTc                                                                                                                                                                                                        <  6:1205378/149‑1 (MQ=255)
                                                                                                                                                                                                     aTCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATaaa                                                                                                                                                                                                               <  6:729095/138‑1 (MQ=255)
                                                                                                                                                                                                      tCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCa                                                                                                                                                                                                   >  8:703476/1‑149 (MQ=255)
                                                                                                                                                                                                        ttAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCata                                                                                                                                                                                                 <  5:1337577/149‑1 (MQ=255)
                                                                                                                                                                                                          tgCCACTTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATTCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCACCGTCAGAAAATGTCCGCGCTCCACGGCGCTGATTTTGGCCATAAAACTGTa                                                                 <  1:13573/274‑1 (MQ=255)
                                                                                                                                                                                                            ccACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACa                                                                                                                                                                                             >  5:500917/1‑149 (MQ=255)
                                                                                                                                                                                                            ccACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACa                                                                                                                                                                                             >  8:1119268/1‑149 (MQ=255)
                                                                                                                                                                                                                gTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGCATATTGCTGAAGATCGTGACgggacga                                                                                                                                                                                                                                     >  7:758228/1‑100 (MQ=18)
                                                                                                                                                                                                                  tttCGAATTACACCGTAATTTAATTGCAACCGGCTGTATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCACCGTCAGAAAATGTCCGCGCTCCACGGCGCTGATTTTGGCCATAAAACTGTAGGGACGat                                                         <  4:224777/274‑1 (MQ=255)
                                                                                                                                                                                                                  cTTCGACTTACACCGTAATTTAATTGAAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGcc                                                                                                                                                                                       <  5:1450815/149‑1 (MQ=255)
                                                                                                                                                                                                                   ttCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAg                                                                                                                                                                                                              >  6:513733/1‑125 (MQ=255)
                                                                                                                                                                                                                   ttCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCt                                                                                                                                                                                      >  6:611187/1‑149 (MQ=255)
                                                                                                                                                                                                                   ttCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCt                                                                                                                                                                                      >  6:461774/1‑149 (MQ=255)
                                                                                                                                                                                                                         ttACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCgg                                                                                                                                                                                >  5:331072/1‑149 (MQ=255)
                                                                                                                                                                                                                          tACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGt                                                                                                                                                                               <  8:1026305/149‑1 (MQ=255)
                                                                                                                                                                                                                            cacCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGc                                                                                                                                                                             <  7:1247008/149‑1 (MQ=255)
                                                                                                                                                                                                                             acCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGTTCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCa                                                                                                                                                                            >  5:281091/1‑149 (MQ=255)
                                                                                                                                                                                                                             acCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCa                                                                                                                                                                            <  8:87051/149‑1 (MQ=255)
                                                                                                                                                                                                                                  aatttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCACCGTCAGAAAATGTCCGCGCTCCACGGCGCTGATTTTGGCCATAAAACTGTAGGGACGATAAATAATGGTTGGAtt                                         >  2:47953/1‑275 (MQ=255)
                                                                                                                                                                                                                                   atttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaaaa                                                                                                                                                                      <  7:836389/149‑1 (MQ=255)
                                                                                                                                                                                                                                         ttGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTc                                                                                                                                                                <  5:420451/149‑1 (MQ=255)
                                                                                                                                                                                                                                             aaCCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCg                                                                                                                                                            >  6:1355088/1‑149 (MQ=255)
                                                                                                                                                                                                                                                   cTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCg                                                                                                                                                      >  7:203694/1‑149 (MQ=255)
                                                                                                                                                                                                                                                        tGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTc                                                                                                                                                 <  7:1131720/149‑1 (MQ=255)
                                                                                                                                                                                                                                                            gggTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCa                                                                                                                                             <  5:677088/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                  aTGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGcc                                                                                                                                       <  7:1186697/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                  aTGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGNTGCCGGATTCACCATTCGCCAGGCAAAATGCGGGACCACCGTCAGAAAATGTCCGCGCTCCACGGCGATNATTTTGGCCATAAAACTGAAGGGACGATAAAAAAGGGTTGGGTTGNTNCGAAAAGGGCAGATATTTGcgacaaaac         >  3:38341/1‑267 (MQ=255)
                                                                                                                                                                                                                                                                   tGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCa                                                                                                                                      <  6:1390620/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                    gACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAg                                                                                                                                     >  7:786683/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    gACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAg                                                                                                                                     >  8:74831/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                       tCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCACCGTCAGAAAATGTCCGCGCTCCACGGCGCTGATTTTGGCCATAAAACTGTAGGGACGATAAATAATGGTTGGATTGATCCCACAAGGTCGCATAtt                    >  4:314587/1‑259 (MQ=255)
                                                                                                                                                                                                                                                                         gTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCCCCGTCCGAAAATGTCCGCTCTCCACGGCGCTGATTTTGGCCATACAAATGTAGGGCAGATAAATAATGGTTGGATTGATCCCCCAAGGTCGCATATTCGAATCAATCATTGCttt  >  2:248907/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                          tCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCaaaa                                                                                                                               >  7:476842/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                           cccGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATaaa                                                                                                                                                                                                               <  5:1058480/68‑1 (MQ=255)
                                                                                                                                                                                                                                                                           cccGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATaaa                                                                                                                                                                                                               >  6:1058480/1‑68 (MQ=255)
                                                                                                                                                                                                                                                                           cccGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAAt                                                                                                                              >  5:1015083/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                             cGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCAc                                                                                                                                               <  5:155330/130‑1 (MQ=255)
                                                                                                                                                                                                                                                                             cGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCAc                                                                                                                                               >  6:155330/1‑130 (MQ=255)
                                                                                                                                                                                                                                                                             cGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGc                                                                                                                            >  6:1102613/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                              gTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCg                                                                                                                           <  6:281651/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                              gTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCg                                                                                                                           <  8:1155174/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                              gTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGCTGCCGGATTCACCAGTCGCCAGGCAAAATGCg                                                                                                                           >  6:586857/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                               tCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATc                                                                                                                                                                                           >  6:1306297/1‑84 (MQ=255)
                                                                                                                                                                                                                                                                               tCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCgg                                                                                                                          <  8:229568/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                             
CGTGCCAATGTGCAGCCATGTTGATCACAGCGAACACAGCGTCAAAGTGATCATCACCGAACAAGGGATCGCCGATCTGCGCGGTCTTTCCCCGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCACCGTCAGAAAATGTCCGCGCTCCACGGCGCTGATTTTGGCCATAAAACTGTAGGGACGATAAATAATGGTTGGATTGATCCCACAAGGTCGCATATTCGCATCAATCATTGCTTC  >  NC_000913/3066022‑3066561

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