Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,472,437 T→C 7.1% N33D (AAT→GAT)  pyrL ← pyrBI operon leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,472,4370TC 7.1% 673.8 / 27.8 212N33D (AAT→GAT) pyrLpyrBI operon leader peptide
Reads supporting (aligned to +/- strand):  ref base T (99/97);  new base C (8/7);  total (107/105)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.27e-01

Only 100 of 214 total aligned reads displayed.

TTTCGAAAGAGAGGCGGGTACGGGTAGAGGCTTCGAAGAAACAGCTGGCAATGACTTTGTGCTTCAACAGCTCTGGTTGCGGGTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGA  >  NC_000913/4472212‑4472706
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                             
                           agGCTTCGAAGAAACAGCTGGCAATGACTTTGTGCTTCAACAGCTCTGGTTGCGGGTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGGGTATTTTTTCAGACGCGGTAAGaaaaaa                                                                                                                                                                                                   >  1:48610/1‑275 (MQ=255)
                               ttCGAAGAAACAGCTGGCAATGACTTTGTGCTTCAACAGCTCTGGTTGCGGGTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGa                                                                                                                                                                                                                                                                             >  4:367478/1‑197 (MQ=255)
                                 cGAAGAAACAGCTGGCAATGACTTTGTGCTTCAACAGCTCTGGTTGCGGGTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTCGGACAAAAAAAGCACCTCGATTGAGTGGCTGGGCATGGGTGATCAACGGGAAGAAAAACGGCAAGCTAGCGTCTTTGTTTAGACGCAGTATGAAAAAATCTCGa                                                                                                                                                                                             >  1:8739/1‑275 (MQ=255)
                                   aagaaACAGCTGGCACTGTCTTTTTGCTGCAACAGCTCTTGTTGCGCGTTTGCTTTCTGTTGTGCCTGTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAAGGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGaa                                                                                                                                                                                            <  2:167292/274‑1 (MQ=255)
                                             tGGCAATGACTTTGTGCTTCAACAGCTCTGGTTGCGGGTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAAc                                                                                                                                                                                  <  2:88228/274‑1 (MQ=255)
                                                                    aGCTCTTGTTGCGGGTTTGCTTTCAGTTTAGCCGCTGTCGCCAGCACCAGATTAAGGTCATGGCGACTAAGGTCGTTTATGGAAATGATATGTTGCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTc                                                                                                                                                          <  2:197996/275‑1 (MQ=255)
                                                                                gggTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAg                                                                                                                                                                                                                                                                            <  8:1139295/149‑1 (MQ=255)
                                                                                  gTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAggg                                                                                                                                                                                                                                                                          >  8:944658/1‑149 (MQ=255)
                                                                                    ttGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCAATCGAGGGGc                                                                                                                                                                                                                                                                        <  5:1057787/149‑1 (MQ=255)
                                                                                      gCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTg                                                                                                                                                                                                                                                                      <  6:990719/149‑1 (MQ=255)
                                                                                       cTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTgg                                                                                                                                                                                                                                                                     >  7:584770/1‑149 (MQ=255)
                                                                                       cTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTgg                                                                                                                                                                                                                                                                     <  7:495929/149‑1 (MQ=255)
                                                                                        tttCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTggg                                                                                                                                                                                                                                                                    >  7:1277918/1‑149 (MQ=255)
                                                                                        tttCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCAATCGAGGGGc                                                                                                                                                                                                                                                                        <  7:570755/145‑1 (MQ=255)
                                                                                        tttCAGTTTAGCCGCTGTCGCCAGCACCAGACTAGGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTggg                                                                                                                                                                                                                                                                    >  6:1121627/1‑149 (MQ=255)
                                                                                           cAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAAt                                                                                                                                                                                                                                                                 <  5:407930/149‑1 (MQ=255)
                                                                                             gTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATgg                                                                                                                                                                                                                                                               <  6:1283010/149‑1 (MQ=255)
                                                                                               ttCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGt                                                                                                                                                                                                                                                             <  8:1065853/149‑1 (MQ=255)
                                                                                                 cgccgcTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGa                                                                                                                                                                                                                                                           <  6:364121/149‑1 (MQ=255)
                                                                                                   ccgcTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATc                                                                                                                                                                                                                                                         <  5:114703/149‑1 (MQ=255)
                                                                                                       tGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACg                                                                                                                                                                                                                                                     >  5:916541/1‑149 (MQ=255)
                                                                                                          tgcaaGCACCAGATTAAGGACATCACAACTAAGTTCGTTTATGCTAATGAAATGTGTCTAAAATAACCGATTAGCTATTTTTTATCCTCTGTCGCATGGGAAAAAAAAAGCCCCTCGATTTAGGGGCTGGGAATGGGTGATGAACGGga                                                                                                                                                                                                                                                  <  6:800043/145‑1 (MQ=255)
                                                                                                          cGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGga                                                                                                                                                                                                                                                  <  6:1075194/149‑1 (MQ=255)
                                                                                                           gCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGgaa                                                                                                                                                                                                                                                 >  7:157534/1‑149 (MQ=255)
                                                                                                           gCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGgaa                                                                                                                                                                                                                                                 <  7:712009/149‑1 (MQ=255)
                                                                                                           gCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATCAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGgaa                                                                                                                                                                                                                                                 <  8:1048625/149‑1 (MQ=255)
                                                                                                            ccAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTCTCTGATATAACGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGgaag                                                                                                                                                                                                                                                <  7:1295724/149‑1 (MQ=255)
                                                                                                                cacaAGATTAAGGTCAACGCGACAAAGGTCGTTGTTGGAAATGATAGGGTTCGGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCGCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCAGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGc                                                                                                               <  2:44083/274‑1 (MQ=255)
                                                                                                                 aCCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGaaaaa                                                                                                                                                                                                                                           >  6:906033/1‑149 (MQ=255)
                                                                                                                 aCCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGaaaaa                                                                                                                                                                                                                                           >  5:654692/1‑149 (MQ=255)
                                                                                                                     gATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGAt                                                                                                                                                                                                                                                          <  7:892284/130‑1 (MQ=255)
                                                                                                                     gATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCacggtctcttatac                                                                                                          >  1:11149/1‑261 (MQ=255)
                                                                                                                      aTTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGt                                                                                                                                                                                                                                                             >  5:1250698/1‑126 (MQ=255)
                                                                                                                      aTTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGt                                                                                                                                                                                                                                                             <  6:1250698/126‑1 (MQ=255)
                                                                                                                        tAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAgg                                                                                                                                                                                                                                    <  6:1235963/149‑1 (MQ=255)
                                                                                                                           ggTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCa                                                                                                                                                                                                                                 >  7:409528/1‑149 (MQ=255)
                                                                                                                             tCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGc                                                                                                                                                                                                                               >  5:719735/1‑149 (MQ=255)
                                                                                                                                tCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGc                                                                                                                                                                                                                                       >  8:634350/1‑138 (MQ=255)
                                                                                                                                tCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTc                                                                                                                                                                                                                            >  8:5511/1‑149 (MQ=255)
                                                                                                                                    gACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCtttt                                                                                                                                                                                                                        <  7:427818/149‑1 (MQ=255)
                                                                                                                                           gTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCt                                                                                                                                                                                                                                                                       >  7:33402/1‑95 (MQ=255)
                                                                                                                                                     gAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGt                                                                                                                                                                                                                                                             >  5:518507/1‑95 (MQ=255)
                                                                                                                                                     gAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGAc                                                                                                                                                                                                       >  6:699301/1‑149 (MQ=255)
                                                                                                                                                      aaaTGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATg                                                                                                                                                                                                                                                                <  6:1258510/91‑1 (MQ=255)
                                                                                                                                                         tGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCa                                                                                                                                                                                                                                 >  8:166691/1‑119 (MQ=255)
                                                                                                                                                         tGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACaaaa                                                                                                                                                                                                   <  7:187132/149‑1 (MQ=255)
                                                                                                                                                          gATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGgaagaa                                                                                                                                                                                                                                              <  8:900118/105‑1 (MQ=255)
                                                                                                                                                          gATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAAGGGGTGATCAACGGGAAGAAAAACGGAAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAAt                                                                                                                                                                                                  >  8:391342/1‑149 (MQ=255)
                                                                                                                                                           atatGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATg                                                                                                                                                                                                 >  6:202373/1‑149 (MQ=255)
                                                                                                                                                             atGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTc                                                                                                                                                                                               <  5:1272471/149‑1 (MQ=255)
                                                                                                                                                               gTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGa                                                                                                                                                                                             >  8:1212718/1‑149 (MQ=255)
                                                                                                                                                                tttCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGaa                                                                                                                                                                                            >  7:870602/1‑149 (MQ=255)
                                                                                                                                                                 ttCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAac                                                                                                                                                                                           >  5:1477466/1‑149 (MQ=255)
                                                                                                                                                                  tcgtatatAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTGTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATg                                                            <  2:287813/271‑1 (MQ=255)
                                                                                                                                                                    tGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAacaca                                                                                                                                                                                        >  8:1261287/1‑149 (MQ=255)
                                                                                                                                                                     gATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTGTTTTTTCAGACGCGGTAAGACAAAATGTCGAacacac                                                                                                                                                                                       <  5:1214048/149‑1 (MQ=255)
                                                                                                                                                                     gATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAacacac                                                                                                                                                                                       <  8:1277918/149‑1 (MQ=255)
                                                                                                                                                                       tataGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTg                                                                                                                                                                                     <  7:1408844/149‑1 (MQ=255)
                                                                                                                                                                         taGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGaa                                                                                                                                                                                   <  8:594218/149‑1 (MQ=255)
                                                                                                                                                                          aGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAAc                                                                                                                                                                                  <  6:1191317/149‑1 (MQ=255)
                                                                                                                                                                           gCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCtttt                                                                                                                                                                                                                        >  6:256769/1‑110 (MQ=255)
                                                                                                                                                                           gCGGATTAGCCATCTGTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCtttt                                                                                                                                                                                                                        <  5:256769/110‑1 (MQ=255)
                                                                                                                                                                              gATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATa                                                                                                                                                                              <  6:528229/149‑1 (MQ=255)
                                                                                                                                                                              gATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATa                                                                                                                                                                              <  7:1090526/149‑1 (MQ=255)
                                                                                                                                                                              gATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATa                                                                                                                                                                              <  5:347399/149‑1 (MQ=255)
                                                                                                                                                                               aTTGGCTATGTTTTATCCCCTGGAGCCTGNNAAAAAAAAAGCCCCTCGATGGGGGGGCTGGNAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACANACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGa                                               <  3:48610/275‑1 (MQ=255)
                                                                                                                                                                                    ccATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACAtcct                                                                                                                                                                        >  8:279472/1‑149 (MQ=255)
                                                                                                                                                                                    ccATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACAtcct                                                                                                                                                                        <  5:1080224/149‑1 (MQ=255)
                                                                                                                                                                                      aTCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACAtcctcc                                                                                                                                                                      >  8:557500/1‑149 (MQ=255)
                                                                                                                                                                                         ttttATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGt                                                                                                                                                                                                >  7:1172144/1‑120 (MQ=255)
                                                                                                                                                                                         ttttATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGt                                                                                                                                                                                                <  8:1172144/120‑1 (MQ=255)
                                                                                                                                                                                         ttttATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCgg                                                                                                                                                                   >  7:820669/1‑149 (MQ=255)
                                                                                                                                                                                         ttttATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCgg                                                                                                                                                                   >  8:1363809/1‑149 (MQ=255)
                                                                                                                                                                                             aTCTCCTGCCGCCTGGGCAAAAAAAAGCCCCTCGATTGGGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCaaa                                                                                                                                                               <  5:688217/149‑1 (MQ=255)
                                                                                                                                                                                                   tGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTcc                                                                                                                                                         <  8:870602/149‑1 (MQ=255)
                                                                                                                                                                                                         cTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAac                                                                                                                                                                                           <  6:475554/109‑1 (MQ=255)
                                                                                                                                                                                                         cTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAac                                                                                                                                                                                           >  5:475554/1‑109 (MQ=255)
                                                                                                                                                                                                            ggCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCgg                  >  4:146857/1‑275 (MQ=255)
                                                                                                                                                                                                             gCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATAc                                                                                                                                               >  6:1268262/1‑149 (MQ=255)
                                                                                                                                                                                                               aaaaaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACtt                                                                                                                                             >  5:1324660/1‑148 (MQ=255)
                                                                                                                                                                                                               aaaaaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTc                                                                                                                                             >  5:1287237/1‑149 (MQ=255)
                                                                                                                                                                                                               aaaaaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTc                                                                                                                                             >  5:988770/1‑149 (MQ=255)
                                                                                                                                                                                                               aaaaaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTc                                                                                                                                             >  5:789872/1‑149 (MQ=255)
                                                                                                                                                                                                                 aaaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCAt                                                                                                                                           >  5:198931/1‑149 (MQ=255)
                                                                                                                                                                                                                 aaaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCAt                                                                                                                                           >  5:560605/1‑149 (MQ=255)
                                                                                                                                                                                                                  tcccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATc                                                                                                                                          >  6:405369/5‑149 (MQ=255)
                                                                                                                                                                                                                    ccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGaa                                                                                                                                                                                            >  6:646955/3‑97 (MQ=255)
                                                                                                                                                                                                                    ccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCTGCATTATACTCATCGt                                                                                                                                        >  7:931632/3‑149 (MQ=255)
                                                                                                                                                                                                                    ccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGt                                                                                                                                        >  8:441656/3‑149 (MQ=255)
                                                                                                                                                                                                                     caGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTc                                                                                                                                       >  5:28308/2‑149 (MQ=255)
                                                                                                                                                                                                                     aaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTc                                                                                                                                       >  5:749003/1‑149 (MQ=255)
                                                                                                                                                                                                                       gCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAg                                                                                                                                     <  7:299994/149‑1 (MQ=255)
                                                                                                                                                                                                                         cccTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGaa                                                                                                                                   >  6:668194/1‑149 (MQ=255)
                                                                                                                                                                                                                           cTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAgc                                                                                                                                 >  8:145923/1‑149 (MQ=255)
                                                                                                                                                                                                                           cTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAgc                                                                                                                                 >  6:1161338/1‑149 (MQ=255)
                                                                                                                                                                                                                           cTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAgc                                                                                                                                 >  5:942664/1‑149 (MQ=255)
                                                                                                                                                                                                                           cTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAgc                                                                                                                                 >  5:881131/1‑149 (MQ=255)
                                                                                                                                                                                                                            tCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAgcg                                                                                                                                <  8:409528/149‑1 (MQ=255)
                                                                                                                                                                                                                             cGATTGAGGGGCTGGGCATGGGTGATCAAAGGGAAGAAAAACGGCAGGGCAGCGTCTCTTTTCAGACTCGGTAAGACAAAGCGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAgcgg                                                                                                                               <  6:1049617/149‑1 (MQ=255)
                                                                                                                                                                                                                             cGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCCTTATACTCATCGTCAGACGCGGCGCAAGCATGTGATGCAATATTTTGTCAGCGCAAACGGTGTATTTGAANTAAAAGCCAAGGTATATGCATCTTTTTTTGATTTGATTTTGAAGGGGAACTTTTTTTCGTTATCaacactttcaaa  >  3:67889/1‑264 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                             
TTTCGAAAGAGAGGCGGGTACGGGTAGAGGCTTCGAAGAAACAGCTGGCAATGACTTTGTGCTTCAACAGCTCTGGTTGCGGGTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGA  >  NC_000913/4472212‑4472706

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: