Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,141,153 A→G 6.2% intergenic (+167/+29) flgL → / ← rne flagellar hook‑filament junction protein/endoribonuclease; RNA‑binding protein;RNA degradosome binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,141,1530AG 6.2% 314.8 / 8.4 97intergenic (+167/+29)flgL/rneflagellar hook‑filament junction protein/endoribonuclease; RNA‑binding protein;RNA degradosome binding protein
Reads supporting (aligned to +/- strand):  ref base A (50/41);  new base G (2/4);  total (52/45)
Fisher's exact test for biased strand distribution p-value = 4.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.44e-01

CATGAAACTGGGTATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCTGCCAGTCACTGTGACGC  >  NC_000913/1141005‑1141299
                                                                                                                                                    |                                                                                                                                                  
cATGAAACTGGGTATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTa                                                                                                                                                    >  1:808520/1‑149 (MQ=255)
     aaCTGGGTATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAgg                                                                                                                                               <  3:267213/149‑1 (MQ=255)
        tGGGTATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGTAACTGCCAGGGCt                                                                                                                                            <  2:847740/149‑1 (MQ=255)
         gggTATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCtt                                                                                                                                           <  2:297549/149‑1 (MQ=255)
          ggTATGTTTTGTCTGCCTGCTCTGGGATCGCCGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGTAACTGCCAGGGCttt                                                                                                                                          <  4:200818/149‑2 (MQ=255)
           gTATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGTAACTGCCAGGGCtttt                                                                                                                                         <  3:574012/149‑3 (MQ=255)
           gTATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGa                                                                                                                                         >  1:645628/1‑149 (MQ=255)
           gTATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGa                                                                                                                                         >  1:420370/1‑149 (MQ=255)
            tATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGAt                                                                                                                                        >  2:219931/1‑149 (MQ=255)
                ttttGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACt                                                                                                                                    >  1:109177/1‑149 (MQ=255)
                   tGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTg                                                                                                                                 <  2:142381/149‑1 (MQ=255)
                    gTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGa                                                                                                                                <  3:631083/149‑1 (MQ=255)
                     tCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAg                                                                                                                               <  1:354094/149‑1 (MQ=255)
                       tgcctgcTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCt                                                                                                                             >  1:746217/1‑149 (MQ=255)
                       tgcctgcTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCt                                                                                                                             >  3:843811/1‑149 (MQ=255)
                         cctgcTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTaa                                                                                                                           >  2:217650/1‑149 (MQ=255)
                           tgcTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGTAACTGCCAGGGCtttttta                                                                                                                                      <  3:291560/136‑6 (MQ=255)
                           tgcTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGTAACTGCCAGGGCtttttta                                                                                                                                      >  4:291560/1‑131 (MQ=255)
                           tgcTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAatt                                                                                                                         >  3:538744/1‑149 (MQ=255)
                                tGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTAc                                                                                                                    <  3:52612/149‑1 (MQ=255)
                                    aTCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCaa                                                                                                                >  3:610589/1‑149 (MQ=255)
                                      cGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAACTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACa                                                                                                              <  2:880835/149‑1 (MQ=255)
                                       gCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAg                                                                                                             >  3:418482/1‑149 (MQ=255)
                                            ggcgggcATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGa                                                                                                                                <  3:237072/125‑1 (MQ=255)
                                            ggcgggcATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGa                                                                                                                                >  4:237072/1‑125 (MQ=255)
                                                ggcATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGAc                                                                                                    >  4:767500/1‑149 (MQ=255)
                                                 gcATTTTTTTGCGTATTTTGCATTGTTGGTGAGGAAGGCTACCTTTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATGATTACTCAACAGGTTGCGGAcg                                                                                                   >  4:510579/1‑149 (MQ=255)
                                                 gcATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGAcg                                                                                                   >  2:195691/1‑149 (MQ=255)
                                                   aTTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGAcgcg                                                                                                 >  3:652744/1‑149 (MQ=255)
                                                   aTTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGAcgcg                                                                                                 >  1:576783/1‑149 (MQ=255)
                                                          tGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGc                                                                                          <  2:813244/149‑1 (MQ=255)
                                                            ccTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCgg                                                                                        <  2:727912/149‑1 (MQ=255)
                                                            ccTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCgg                                                                                        >  4:459767/1‑149 (MQ=255)
                                                             cTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAg                                                                                              <  1:86265/142‑1 (MQ=255)
                                                             cTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAg                                                                                              >  2:86265/1‑142 (MQ=255)
                                                             cTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGc                                                                                       >  2:405269/1‑149 (MQ=255)
                                                             cTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGc                                                                                       >  2:498254/1‑149 (MQ=255)
                                                              tATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCa                                                                                      >  3:857421/1‑149 (MQ=255)
                                                               aTTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCag                                                                                     >  4:92074/1‑149 (MQ=255)
                                                                ttttGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCaga                                                                                    <  1:405269/149‑1 (MQ=255)
                                                                 tttGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCagag                                                                                   <  3:90618/149‑1 (MQ=255)
                                                                 tttGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCagag                                                                                   <  1:658824/149‑1 (MQ=255)
                                                                    gCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCa                                                                                >  2:240852/1‑149 (MQ=255)
                                                                      aTTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAAttatta                                                                                                                     >  1:171311/1‑110 (MQ=255)
                                                                      aTTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAAttatta                                                                                                                     <  2:171311/110‑1 (MQ=255)
                                                                      aTTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCatg                                                                              >  4:679889/1‑149 (MQ=255)
                                                                       ttgttgGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGc                                                                                                                                             >  1:372683/1‑85 (MQ=255)
                                                                       ttgttgGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGc                                                                                                                                             <  2:372683/85‑1 (MQ=255)
                                                                       ttgttgGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCatga                                                                             <  4:418482/149‑1 (MQ=255)
                                                                            ggTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATgtgt                                                                        <  3:846335/149‑1 (MQ=255)
                                                                             gTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGtt                                                                       <  1:498254/149‑1 (MQ=255)
                                                                                aGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGcc                                                                    <  4:652744/149‑1 (MQ=255)
                                                                                 gCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCg                                                                   >  2:725489/1‑149 (MQ=255)
                                                                                                 aTGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGc                                                   >  1:246152/1‑149 (MQ=255)
                                                                                                   gAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCAc                                                 <  1:725489/149‑1 (MQ=255)
                                                                                                     aTTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCt                                               <  3:510579/149‑1 (MQ=255)
                                                                                                      ttttAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCtt                                              >  3:668505/1‑149 (MQ=255)
                                                                                                      ttttAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTCGCCGTATGACCACCTGCGGCACCtt                                              >  2:357967/1‑149 (MQ=255)
                                                                                                         tAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTAc                                           <  2:348784/149‑1 (MQ=255)
                                                                                                             atGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGAcgc                                                                                                  >  1:350510/1‑90 (MQ=255)
                                                                                                             atGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGAcgc                                                                                                  <  2:350510/90‑1 (MQ=255)
                                                                                                               gTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGa                                     >  2:368778/1‑149 (MQ=255)
                                                                                                               gTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGa                                     >  1:620334/1‑149 (MQ=255)
                                                                                                                ttGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGaa                                    <  2:576783/149‑1 (MQ=255)
                                                                                                                ttGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGaa                                    <  4:668505/149‑1 (MQ=255)
                                                                                                                ttGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTCGCCGTATGACCACCTGCGGCACCTTTACCTTCGaa                                    <  1:357967/149‑1 (MQ=255)
                                                                                                                   aTTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGc                                 <  2:246152/149‑1 (MQ=255)
                                                                                                                    ttCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCa                                >  4:296118/1‑149 (MQ=255)
                                                                                                                    ttCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCa                                >  2:583434/1‑149 (MQ=255)
                                                                                                                     tCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCaa                               >  2:39207/1‑149 (MQ=255)
                                                                                                                     tCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCaa                               >  4:129140/1‑149 (MQ=255)
                                                                                                                       aaaGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTc                                                                                                                  >  2:400675/1‑64 (MQ=255)
                                                                                                                       aaaGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTc                                                                                                                  <  1:400675/64‑1 (MQ=255)
                                                                                                                       aaaGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCaaaa                             <  1:195691/149‑1 (MQ=255)
                                                                                                                       aaaGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCaaaa                             <  2:746217/149‑1 (MQ=255)
                                                                                                                       aaaGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCaaaa                             <  4:843811/149‑1 (MQ=255)
                                                                                                                          gATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTa                          >  1:350295/1‑149 (MQ=255)
                                                                                                                          gATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTa                          <  4:368132/149‑1 (MQ=255)
                                                                                                                          gATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTa                          <  3:459767/149‑1 (MQ=255)
                                                                                                                           aTGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAg                         >  2:618051/1‑149 (MQ=255)
                                                                                                                           aTGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAg                         <  2:592453/149‑1 (MQ=255)
                                                                                                                               aaTAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGgcg                     >  1:610760/1‑149 (MQ=255)
                                                                                                                                aTAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCatgat                                                                            >  4:589328/1‑93 (MQ=255)
                                                                                                                                aTAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCatgat                                                                            <  3:589328/93‑1 (MQ=255)
                                                                                                                                aTAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGgcgc                    >  3:526939/1‑149 (MQ=255)
                                                                                                                                aTAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGgcgc                    >  2:308887/1‑149 (MQ=255)
                                                                                                                                aTAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGgcgc                    <  2:620334/149‑1 (MQ=255)
                                                                                                                                   aaaaaGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCTGc                 <  2:574995/149‑1 (MQ=255)
                                                                                                                                    tcaaGCCCTGGTAACTGCCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCTGcc                >  2:40613/3‑149 (MQ=255)
                                                                                                                                     aaaGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGa                                     <  1:464772/127‑1 (MQ=255)
                                                                                                                                     aaaGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGa                                     >  2:464772/1‑127 (MQ=255)
                                                                                                                                        gCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCTGCCAGTc            >  1:715923/1‑149 (MQ=255)
                                                                                                                                            tGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCTGCCAGTCACtg        >  3:611947/1‑149 (MQ=255)
                                                                                                                                            tGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCTGCCAGTCACtg        <  3:440339/149‑1 (MQ=255)
                                                                                                                                             ggCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCTGCCAGTCACtgt       <  1:217650/149‑1 (MQ=255)
                                                                                                                                                aGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCTGCCAGTCACTGTGAc    >  3:459746/1‑149 (MQ=255)
                                                                                                                                                aGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCTGCCAGTCACTGTGAc    >  4:568150/1‑149 (MQ=255)
                                                                                                                                                 gTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTACCTTCGAAGGCAAAAGTAGGGCGCTGCCAGTCACTGTGAcgc  >  2:749416/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CATGAAACTGGGTATGTTTTGTCTGCCTGCTCTGGGATCGCTGGGGCGGGCATTTTTTTGCCTATTTTGCATTGTTGGTTAGCAAGGATGCCATTCGATGAATTTTAATATGTTGATTCAAAGATGAAATAAAAAAGCCCTGGCAGTTACCAGGGCTTGATTACTTTGAGCTAATTATTACTCAACAGGTTGCGGACGCGCAGGAGCGGCAGAGGCATGATGTGTTGCCGTATGACCACCTGCGGCACCTTTACCTTCGAAGGCAAAAGTAGGGCGCTGCCAGTCACTGTGACGC  >  NC_000913/1141005‑1141299

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