Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 562,105 C→T 10.2% N803N (AAC→AAT sfmD → putative outer membrane export usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913562,1050CT10.2% 538.4 / 37.7 176N803N (AAC→AATsfmDputative outer membrane export usher protein
Reads supporting (aligned to +/- strand):  ref base C (74/84);  new base T (5/13);  total (79/97)
Fisher's exact test for biased strand distribution p-value = 1.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 179 total aligned reads displayed.

ATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACA  >  NC_000913/561956‑562251
                                                                                                                                                     |                                                                                                                                                  
aTCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAAc                                                                                                                                                    >  2:2696430/1‑150 (MQ=255)
 tCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAAcg                                                                                                                                                   >  1:1537824/1‑149 (MQ=255)
    aTGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTa                                                                                                                                                <  1:132661/150‑1 (MQ=255)
       ttGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCg                                                                                                                                             <  1:1091100/150‑1 (MQ=255)
       ttGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCg                                                                                                                                             <  1:2516497/150‑1 (MQ=255)
           aaaaaaTATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCggcaa                                                                                                                                         >  2:268984/1‑147 (MQ=255)
            aaaaaTATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACAGGAATa                                                                                                                                        >  1:290805/1‑150 (MQ=255)
              aaaTATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAAc                                                                                                                                      <  2:1864854/150‑1 (MQ=255)
              aaaTATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAAc                                                                                                                                      <  1:856538/150‑1 (MQ=255)
              aaaTATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAAc                                                                                                                                      <  1:2228933/150‑1 (MQ=255)
                atatTAGCAGCGTTGTGCCGACGCACGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCa                                                                                                                                    <  2:2052138/150‑1 (MQ=255)
                atatTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCa                                                                                                                                    <  1:2196237/150‑1 (MQ=255)
                atatTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCa                                                                                                                                    <  1:881197/150‑1 (MQ=255)
                atatTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCa                                                                                                                                    <  1:881157/150‑1 (MQ=255)
                      gcagcGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTAtgg                                                                                                                              <  2:2319216/150‑1 (MQ=255)
                      gcagcGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTAtgg                                                                                                                              <  1:403821/150‑1 (MQ=255)
                       cagcGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATAGAGATACCCGCATAGCCGTGCGGGCGCTCATTACCCTTACCCAGGGCGGAAAACCGGTGCCGGTTGGATCACGCGTCCGGGAAAACAGTACCCGACTAACCAGAAgggt                                                                                                                             >  2:1629369/1‑150 (MQ=255)
                        agcGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTAtggtg                                                                                                                            <  2:1019668/150‑1 (MQ=255)
                        agcGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTAtggtg                                                                                                                            >  2:2497107/1‑150 (MQ=255)
                              gtgCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGTGtgat                                                                                                                      <  1:2437914/150‑1 (MQ=255)
                                 ccGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGAc                                                                                                                   >  2:684481/1‑150 (MQ=255)
                                 ccGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGAc                                                                                                                   >  1:1894646/1‑150 (MQ=255)
                                  cGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACg                                                                                                                  <  1:559348/150‑1 (MQ=255)
                                  cGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACg                                                                                                                  >  2:1527257/1‑150 (MQ=255)
                                  cGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACg                                                                                                                  <  1:949122/150‑1 (MQ=255)
                                   gACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTCCGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACgg                                                                                                                 >  2:1722526/1‑150 (MQ=255)
                                      gCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAATAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCa                                                                                                              <  2:1898641/150‑1 (MQ=255)
                                       cAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAATAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCaa                                                                                                             <  2:2482168/150‑1 (MQ=255)
                                       cAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCaa                                                                                                             >  2:1288991/1‑150 (MQ=255)
                                          ggCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGtt                                                                                                          <  1:2086097/150‑1 (MQ=255)
                                            cgcgTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGttta                                                                                                        <  2:181900/150‑1 (MQ=255)
                                             gcgTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGtttat                                                                                                       >  1:557450/1‑150 (MQ=255)
                                             gcgTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGgttat                                                                                                       >  2:1936396/1‑150 (MQ=255)
                                               gTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGtttattt                                                                                                     >  1:2191947/1‑150 (MQ=255)
                                                       gTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAgtggtg                                                                                             >  1:1989686/1‑150 (MQ=255)
                                                        tGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTgc                                                                                            >  2:2042588/1‑150 (MQ=255)
                                                            aaTTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCa                                                                                        >  1:882590/1‑150 (MQ=255)
                                                             aTTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCAt                                                                                       >  2:168567/1‑150 (MQ=255)
                                                                  gATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCt                                                                                  <  1:2589280/150‑1 (MQ=255)
                                                                     aCCCGCATAGGCGTGGGGGCGCTCAGTACCGGTACCCCGGGCGGAAAACCGGTGCCGTTTGGATCACGGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGt                                                                               <  1:268984/150‑1 (MQ=255)
                                                                     aCCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGt                                                                               <  2:1682392/150‑1 (MQ=255)
                                                                       ccGCATTGGCGTGCGGGCGCTCATTACCGTTGCCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACGGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTGATTTACGTGGTGCGCCATTGTCTGGTGa                                                                             >  2:1395671/1‑150 (MQ=255)
                                                                       ccGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGa                                                                             >  2:1715734/1‑150 (MQ=255)
                                                                        cGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGaa                                                                            >  2:2246079/1‑150 (MQ=255)
                                                                        cGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGaa                                                                            <  1:1466017/150‑1 (MQ=255)
                                                                        cGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGaa                                                                            >  2:1056758/1‑150 (MQ=255)
                                                                        cGCAGAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGaa                                                                            <  1:379652/150‑1 (MQ=255)
                                                                         gCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAAt                                                                           <  1:2725463/150‑1 (MQ=255)
                                                                           aTAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTa                                                                         <  1:1832644/150‑1 (MQ=255)
                                                                              ggCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTg                                                                      <  1:2181566/150‑1 (MQ=255)
                                                                               gCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTgg                                                                     >  2:240450/1‑150 (MQ=255)
                                                                                 atGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAATAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGtt                                                                   <  1:823778/149‑1 (MQ=255)
                                                                                  tGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTc                                                                  >  2:2016654/1‑150 (MQ=255)
                                                                                    cGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAg                                                                <  2:970315/150‑1 (MQ=255)
                                                                                     gggCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGt                                                               <  2:2368918/150‑1 (MQ=255)
                                                                                     gggCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGt                                                               <  2:1620809/150‑1 (MQ=255)
                                                                                       gcgcTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAGTTACTGGTTCAGTgg                                                             <  1:1015142/150‑1 (MQ=255)
                                                                                       gcgcTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTgg                                                             >  1:1831754/1‑150 (MQ=255)
                                                                                       gcgcTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTgg                                                             <  1:674242/150‑1 (MQ=255)
                                                                                        cgcTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAATAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTggg                                                            >  2:615254/1‑150 (MQ=255)
                                                                                         gcTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTgggg                                                           >  1:1239503/1‑150 (MQ=255)
                                                                                           tCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGgag                                                         >  2:2129405/1‑150 (MQ=255)
                                                                                           tCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGgag                                                         >  2:1859648/1‑150 (MQ=255)
                                                                                                  cGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAATAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTCAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGgcgc                                                  <  2:1511709/150‑1 (MQ=255)
                                                                                                     tACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGaa                                               <  2:1197902/150‑1 (MQ=255)
                                                                                                     tACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGaa                                               <  2:290805/150‑1 (MQ=255)
                                                                                                     tACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGaa                                               <  1:2639587/150‑1 (MQ=255)
                                                                                                     tACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGaa                                               <  2:1296662/150‑1 (MQ=255)
                                                                                                        ccAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTc                                            <  1:2042588/150‑1 (MQ=255)
                                                                                                            ggCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGc                                        >  1:1812864/1‑150 (MQ=255)
                                                                                                               ggAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGc                                     >  2:496933/1‑150 (MQ=255)
                                                                                                                 aaaaCCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCAt                                   >  1:1745020/1‑150 (MQ=255)
                                                                                                                   aaCCGGTGCCGTTTGGATCACTGGTACGGGAAAATAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTg                                 <  2:1725215/150‑1 (MQ=255)
                                                                                                                       ggTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCa                             <  2:759087/150‑1 (MQ=255)
                                                                                                                          gCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTa                          <  1:1278240/150‑1 (MQ=255)
                                                                                                                           ccGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTAt                         >  1:1879309/1‑150 (MQ=255)
                                                                                                                           ccGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTAt                         <  1:1021670/150‑1 (MQ=255)
                                                                                                                            cGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATg                        >  2:1272465/1‑150 (MQ=255)
                                                                                                                              tttGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTa                      >  2:936049/1‑150 (MQ=255)
                                                                                                                              tttGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTa                      >  1:1962838/1‑150 (MQ=255)
                                                                                                                                 ggATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTg                   <  2:2443491/150‑1 (MQ=255)
                                                                                                                                   aTCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGcc                 >  2:27945/1‑150 (MQ=255)
                                                                                                                                   aTCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGcc                 >  2:28175/1‑150 (MQ=255)
                                                                                                                                   aTCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGcc                 >  1:2436751/1‑150 (MQ=255)
                                                                                                                                    tCACTGGTACGGGAAAATAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCg                >  2:354604/1‑150 (MQ=255)
                                                                                                                                      aCTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGaa              >  1:441672/1‑150 (MQ=255)
                                                                                                                                        tGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGc            >  1:338216/1‑150 (MQ=255)
                                                                                                                                        tGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGc            >  1:1587160/1‑150 (MQ=255)
                                                                                                                                        tGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGc            >  1:1728169/1‑150 (MQ=255)
                                                                                                                                          gTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCaa          >  1:1378115/1‑150 (MQ=255)
                                                                                                                                           tACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCaaa         >  2:1109352/1‑150 (MQ=255)
                                                                                                                                           tACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCaaa         >  1:739513/1‑150 (MQ=255)
                                                                                                                                            aCGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGTATTGCGCACTATGTATTGCCGAAGCAAAg        <  1:240450/150‑1 (MQ=255)
                                                                                                                                             cGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGc       <  2:1305809/150‑1 (MQ=255)
                                                                                                                                             cGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGc       >  2:1944341/1‑150 (MQ=255)
                                                                                                                                             cGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGc       >  2:1153512/1‑150 (MQ=255)
                                                                                                                                                gAAAATAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTa    <  1:1122941/150‑1 (MQ=255)
                                                                                                                                                gAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTa    >  2:2589439/1‑150 (MQ=255)
                                                                                                                                                  aaaCAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACa  <  1:1280564/150‑1 (MQ=255)
                                                                                                                                                  aaaCAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACa  <  1:684481/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
ATCGATGTTGAAAAAAATATTAGCAGCGTTGTGCCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACA  >  NC_000913/561956‑562251

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