Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,026,541 T→C 7.4% intergenic (+16/+16) helD → / ← mgsA DNA helicase IV/methylglyoxal synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,026,5410TC 7.4% 474.3 / 21.1 162intergenic (+16/+16)helD/mgsADNA helicase IV/methylglyoxal synthase
Reads supporting (aligned to +/- strand):  ref base T (84/65);  new base C (6/6);  total (90/72)
Fisher's exact test for biased strand distribution p-value = 7.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

Only 100 of 162 total aligned reads displayed.

AGAACGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCG  >  NC_000913/1026393‑1026687
                                                                                                                                                    |                                                                                                                                                  
agaaCGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTg                                                                                                                                                   >  1:2013979/1‑150 (MQ=255)
  aaCGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGcc                                                                                                                                                 >  1:1504936/1‑150 (MQ=255)
     ggcggTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCAcc                                                                                                                                              <  2:1967971/150‑1 (MQ=255)
        ggTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCAAACGCCACCtgt                                                                                                                                           <  1:563803/150‑1 (MQ=255)
        ggTTAATGTACGTGGCGCGGACCGGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCAAACGCCACCtgt                                                                                                                                           <  2:2472606/150‑1 (MQ=255)
          ttAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGc                                                                                                                                         >  2:2115632/1‑150 (MQ=255)
           tAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCa                                                                                                                                        <  1:1027660/150‑1 (MQ=255)
           tAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATAATGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCa                                                                                                                                        <  1:2035012/150‑1 (MQ=255)
            aaTGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCAAACGCCACCtgtttct                                                                                                                                       <  1:1992468/150‑5 (MQ=255)
             aTGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAat                                                                                                                                      >  2:1361867/1‑150 (MQ=255)
               gTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAatat                                                                                                                                    >  1:196353/1‑150 (MQ=255)
                 aCGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTa                                                                                                                                  >  2:25884/1‑150 (MQ=255)
                  cGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTAc                                                                                                                                 <  1:84823/150‑1 (MQ=255)
                    tGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACtt                                                                                                                               <  1:159820/150‑1 (MQ=255)
                     ggCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTc                                                                                                                              >  2:2406605/1‑150 (MQ=255)
                      gcgcTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCa                                                                                                                             <  2:389127/150‑1 (MQ=255)
                      gcgcTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCa                                                                                                                             <  1:1588685/150‑1 (MQ=255)
                       cgcTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAg                                                                                                                            >  2:1763031/1‑150 (MQ=255)
                       cgcTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAg                                                                                                                            >  1:453264/1‑150 (MQ=255)
                           gACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACgg                                                                                                                        <  1:499523/150‑1 (MQ=255)
                               cGGGCACGCCATCGGGTATTGGCACTGTATAACAAAGTGAATCCCTCACCCTTTGTGTAAATACTGAAAAATCTGGATGTTCCGTTGGCGAGAAAACCGTAAGAACCAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCcg                                                                                                                    <  2:2353985/150‑1 (MQ=255)
                                 ggCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCcgcg                                                                                                                  <  2:1989462/150‑1 (MQ=255)
                                    aCGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCgaga                                                                                                               >  1:2028458/1‑150 (MQ=255)
                                     cGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGAt                                                                                                              <  2:1536199/150‑1 (MQ=255)
                                      gCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATa                                                                                                             >  1:1564274/1‑150 (MQ=255)
                                       ccATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATaa                                                                                                            >  1:1670523/1‑150 (MQ=255)
                                        cATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAAc                                                                                                           >  2:617662/1‑150 (MQ=255)
                                              gTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGAt                                                                                                     >  1:854215/1‑150 (MQ=255)
                                              gTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGAt                                                                                                     <  2:2013979/150‑1 (MQ=255)
                                              gTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGAt                                                                                                     <  2:1772167/150‑1 (MQ=255)
                                                  gggCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATc                                                                                                 <  2:1446643/150‑1 (MQ=255)
                                                     cACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCggg                                                                                              >  1:452868/1‑150 (MQ=255)
                                                       cTGTTTAACAAAGGGAAGCCCTCTCCCTTTGTGGATATACTGAAAAATCTGGATGTGCCGGTGGCGTGAAAAGCGTACGAAACAGGTGGCGTTTGCCACCGGGGCAATATTACTTCAGACGGTCCGGGAGAGACCGCTGCTAATCggggg                                                                                            >  2:1975208/1‑149 (MQ=255)
                                                       cTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACGTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGa                                                                                            <  1:1403635/150‑1 (MQ=255)
                                                         gTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATc                                                                                          >  2:1380262/1‑150 (MQ=255)
                                                         gTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATc                                                                                          >  1:371742/1‑150 (MQ=255)
                                                          tttAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCa                                                                                         <  1:2235011/150‑1 (MQ=255)
                                                            tAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGa                                                                                       >  1:1824147/1‑150 (MQ=255)
                                                            tAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGa                                                                                       >  1:1428203/1‑150 (MQ=255)
                                                            tAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGa                                                                                       >  2:288263/1‑150 (MQ=255)
                                                               cAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAata                                                                                    >  1:896556/1‑150 (MQ=255)
                                                                aaaGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAatat                                                                                   >  1:1527479/1‑150 (MQ=255)
                                                                    agaATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGAc                                                                               <  2:221632/150‑1 (MQ=255)
                                                                    agaATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGAc                                                                               <  1:593090/150‑1 (MQ=255)
                                                                      aaTCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACcg                                                                             >  2:644733/1‑150 (MQ=255)
                                                                           ctctCCGATTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGcgtcg                                                                        <  2:849959/150‑1 (MQ=255)
                                                                           ctctCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGcgtcg                                                                        >  2:2100800/1‑150 (MQ=255)
                                                                            tctcCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGcgtcgt                                                                       <  1:1906316/150‑1 (MQ=255)
                                                                            tctcCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGcgtcgt                                                                       <  1:1164412/150‑1 (MQ=255)
                                                                            tctcCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCAAACGCCACCtgt                                                                                                                                           >  1:579085/1‑82 (MQ=38)
                                                                            tctcCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCAAACGCCACCtgt                                                                                                                                           <  2:579085/82‑1 (MQ=38)
                                                                             ctcCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGtt                                                                      <  2:1580157/150‑1 (MQ=255)
                                                                              tcCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCGGAATATCGACCGCGTCGTTg                                                                     >  2:101291/1‑150 (MQ=255)
                                                                              tcCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTg                                                                     <  2:453264/150‑1 (MQ=255)
                                                                                ccTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGaa                                                                   >  2:1472956/1‑150 (MQ=255)
                                                                                ccTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGaa                                                                   <  1:1683980/150‑1 (MQ=255)
                                                                                    tgtgGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATgc                                                               >  2:2305594/1‑150 (MQ=255)
                                                                                    tgtgGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATgc                                                               <  1:345862/150‑1 (MQ=255)
                                                                                     gtgGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATgcg                                                              <  1:25884/150‑1 (MQ=255)
                                                                                     gtgGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATgcg                                                              >  1:2007062/1‑150 (MQ=255)
                                                                                      tgGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATgcgg                                                             >  2:2281849/1‑150 (MQ=255)
                                                                                       ggAAATACTGAAAAATCTGGAGGTGCCGGTGCCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATgcggc                                                            <  2:1928470/150‑1 (MQ=255)
                                                                                        gAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATgcggcg                                                           >  1:1288317/1‑150 (MQ=255)
                                                                                          aaTACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGAc                                                         <  1:1314864/150‑1 (MQ=255)
                                                                                            tACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTg                                                       <  2:2028458/150‑1 (MQ=255)
                                                                                             aCTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTgg                                                      <  1:1664611/150‑1 (MQ=255)
                                                                                              cTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTggc                                                     <  1:149965/150‑2 (MQ=255)
                                                                                              cTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGa                                                     >  1:1128809/1‑150 (MQ=255)
                                                                                              cTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGa                                                     >  2:2472253/1‑150 (MQ=255)
                                                                                               tGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATg                                                                <  1:2481267/138‑1 (MQ=255)
                                                                                                gAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGAtt                                                   <  2:371742/150‑1 (MQ=255)
                                                                                                  aaaaTCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTAt                                                 >  2:1897115/1‑150 (MQ=255)
                                                                                                  aaaaTCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTAt                                                 >  1:700767/1‑150 (MQ=255)
                                                                                                    aaTCTGGATGTGCCGGTGGCAAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGa                                               >  2:2029254/1‑150 (MQ=255)
                                                                                                     aTCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGaa                                              >  1:1873243/1‑150 (MQ=255)
                                                                                                     aTCTGGATGTGCCCGTGGCGAGAAAACCGAAAGAAACAGGTGGCGTTTGCCACCTGTGCAATAATAATTCAGAAGGTCCGCGAGATAACGCTCATAATCGGGGATCAGAAAATCGACCGCGTCCTTGAAATGCGGCGAATGGAATATGaa                                              >  2:2393939/1‑150 (MQ=255)
                                                                                                        tGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTc                                           >  1:224797/1‑150 (MQ=255)
                                                                                                          gATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTg                                         <  2:1284206/150‑1 (MQ=255)
                                                                                                              tgCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGt                                     >  2:424010/1‑150 (MQ=255)
                                                                                                               gCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGc                                        >  1:2359068/1‑146 (MQ=255)
                                                                                                                ccGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTg                                   <  2:1504936/150‑1 (MQ=255)
                                                                                                                       cGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGtt                            >  2:118720/1‑150 (MQ=255)
                                                                                                                       cGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGtt                            >  2:1043533/1‑150 (MQ=255)
                                                                                                                       cGAGAAAAACGTAAGAAACAGGTGGCGTTTGCCAACTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATCCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGtt                            >  2:1806910/1‑150 (MQ=255)
                                                                                                                         agaAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTtgg                          <  1:1472956/150‑1 (MQ=255)
                                                                                                                          gaAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTtggt                         >  1:142880/1‑150 (MQ=255)
                                                                                                                           aaaaCCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTtggtg                        <  1:2100800/150‑1 (MQ=255)
                                                                                                                            aaaCCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTtggtgg                       <  2:2007062/150‑1 (MQ=255)
                                                                                                                             aaCCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGc                      <  1:1380262/150‑1 (MQ=255)
                                                                                                                              aCCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCg                     >  2:333169/1‑150 (MQ=255)
                                                                                                                              aCCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCg                     <  2:1861821/150‑1 (MQ=255)
                                                                                                                                 gTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGAcc                  >  1:603326/1‑150 (MQ=255)
                                                                                                                                     gAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGaa              >  2:1133058/1‑150 (MQ=255)
                                                                                                                                      gcaCAGGTGGCAAACGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAAt             >  1:2433819/3‑150 (MQ=255)
                                                                                                                                        aCAGGTGGCAAACGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGt           >  2:272992/1‑150 (MQ=255)
                                                                                                                                         cAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGtt          >  2:1445183/1‑150 (MQ=255)
                                                                                                                                         cAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGtt          >  1:1636831/1‑150 (MQ=255)
                                                                                                                                            gTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCa       >  1:1498674/1‑150 (MQ=255)
                                                                                                                                             tGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCAt      >  2:18839/1‑150 (MQ=255)
                                                                                                                                                 gTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCg  <  2:1858880/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AGAACGGCGGTTAATGTACGTGGCGCTGACCCGGGCACGCCATCGGGTATGGGCACTGTTTAACAAAGAGAATCCCTCTCCCTTTGTGGAAATACTGAAAAATCTGGATGTGCCGGTGGCGAGAAAACCGTAAGAAACAGGTGGCGTTTGCCACCTGTGCAATATTACTTCAGACGGTCCGCGAGATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCG  >  NC_000913/1026393‑1026687

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