Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,892,270 A→G P19P (CCA→CCG queD → 6‑pyruvoyl tetrahydrobiopterin synthase (PTPS)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,892,2700AG95.4% 380.4 / ‑4.3 109P19P (CCA→CCGqueD6‑pyruvoyl tetrahydrobiopterin synthase (PTPS)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (43/61);  minor base C (5/0);  total (48/61)
Fisher's exact test for biased strand distribution p-value = 1.47e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 110 total aligned reads displayed.

CAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCG  >  NC_000913/2892122‑2892416
                                                                                                                                                     |                                                                                                                                                  
cAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACcgca                                                                                                                                                  <  3:16098/151‑1 (MQ=255)
        ttGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCTGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCg                                                                                                                                           <  2:1047349/150‑1 (MQ=255)
         tGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCgg                                                                                                                                          <  1:838040/150‑1 (MQ=255)
           cgcAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTTTCGTCCTGTAGAGAAATTATGATGTCCACCACGGTTTTTTAAGATTTTACCTTCGAAACCGCCCACCCCCTACCGCACGTccccg                                                                                                                                          >  1:1699/1‑146 (MQ=255)
           cgcAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAg                                                                                                                                       <  4:378550/151‑1 (MQ=255)
             caacaaAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAgg                                                                                                                                      <  2:1268288/150‑1 (MQ=255)
              aacaaAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTTTTGTCCTGTAGAGAAATTATGATGTCCACCACCTTATTTTAAGATTTCACCTTTTAAACCGCCCACCCCTTACCGCCAGTCCCCGGAgggg                                                                                                                                    >  2:101907/1‑150 (MQ=255)
               acaaAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTTAAGATTTCACCTTCGAAGCCGCTCACCCCTTACCCCACGTCCCGGAAggg                                                                                                                                     >  1:135294/1‑149 (MQ=255)
                caaAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGGTTTATCGTTCTTTAAAGAAATTATGATGTCCACCCCGGTTTTTAAAGAATTCCCCTTCCCAACCGCTCACCGCTTACCCCCCCGCCCGGacagg                                                                                                                                     >  3:185141/1‑144 (MQ=255)
                 aaaGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATTGTCCTGTAGAGAAAATATTATGTCCACCACGGTATTTTAAGATTTTACCTTCCAAGCCGCTCACCGCCTTCCGCCCGTCCCGGAAGGGCa                                                                                                                                   >  2:1100348/1‑149 (MQ=255)
                 aaaGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCCGTAGAGAAATTATGATGTCCACCACGTTATTTAAAAATTTTACCTTTGAAGCCGCCCACCCCTTACCGCACGTCCCGGGAGGGGAt                                                                                                                                  >  1:186974/1‑150 (MQ=255)
                      cGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATTCTTCTGTAAAGAAAATATGATGTCCCCCAAGTTATTTTAAGAATTCACCTTTAAACCCGCCCACCCCCTACCCCCCCGCCCGGGAGGGCATaaa                                                                                                                               >  1:11343/1‑148 (MQ=255)
                      cGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAAtgt                                                                                                                            <  1:825651/151‑1 (MQ=255)
                       gTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATTATGTCCACCACGTTATTTAAAGATTTTACCTTCGAAGCCGCTCACCCCTTACCCCCCGTCCCGGAAGGGCATAAAtgtg                                                                                                                           >  2:82158/1‑151 (MQ=255)
                           tAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTc                                                                                                                        <  1:467924/150‑1 (MQ=255)
                                aaTGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTg                                                                                                                   <  1:587858/150‑1 (MQ=255)
                                 aTGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTTTTGTCCTGTAGAGAAAATATGATGTTCCCCCCGCTATTTTAAGATTTTACCTTCGAAGCCGCCCACCGCTTTCCGCACGTCCCCGGAGGGCATAAATGTGGGCGCCt                                                                                                                    >  1:618090/1‑148 (MQ=255)
                                 aTGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCa                                                                                                                 <  2:929219/151‑1 (MQ=255)
                                   gCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCAc                                                                                                                <  1:581713/150‑1 (MQ=255)
                                   gCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCAc                                                                                                                <  1:500480/150‑1 (MQ=255)
                                   gCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACg                                                                                                               <  2:159952/151‑1 (MQ=255)
                                    cGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTTTCGTCCTGTAGAGGAAATTTGATTTCCACCCCGGTTTTTTAAAGTTTTCCCTTCGGAGCCGCTCACCCCTTTCCGCACGTCCCGGGAGGGGATAAAAGAGGTGGCCTGCa                                                                                                                 >  2:878612/1‑148 (MQ=255)
                                    cGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTTTCGTCCTGTAGAAAAAATATTATGTCCACCCCGTTATTTTAGGATTTTACCTTCTAAGCCGCCCACCCCTTACCCCCCCTCCCCGAAAGGCAAAAATGTGGGTGCCCGCACgg                                                                                                              >  2:121267/1‑151 (MQ=255)
                                     gAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACggg                                                                                                             <  2:726728/151‑1 (MQ=255)
                                           aCAGCCTTTCCGGTACTCTACGGCGGTTTTTTCGTCCTGTAGAGAAATTATGATTTCCACCCCGGTTTTTAAAGATTTTACCTTCGAAGCCGCCCACCCCTTACCGCACGTCCCGGGAGGGCATAAATGTGGGCGCCTGCACGGGCATTcc                                                                                                       >  1:423913/1‑151 (MQ=255)
                                             aGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCtt                                                                                                     <  2:1073616/151‑1 (MQ=255)
                                             aGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCtt                                                                                                     <  1:296070/151‑1 (MQ=255)
                                                   tCCGGTACTCTACGGCGGATTTTTCGTCCTGTAGAGAAATTTTTATGTCCCCCACGTTATTTAAAGATTTTCCCTTCGAAGCCCCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACCGGCATTCCTTTAt                                                                                                  >  1:470767/1‑148 (MQ=255)
                                                    ccGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGt                                                                                               <  2:1133361/150‑1 (MQ=255)
                                                      ggTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCg                                                                                            <  2:81382/151‑1 (MQ=255)
                                                          ctctACGGCGGTTTTTTCGTTCTGTAGAGAAATTATGATTTTCACCACGTTATTTAAAGATTTTACCTTTGAAAACGCTCACCGCTTACCGCACGTCCCGGGAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTTTGGTGCGACTg                                                                                        >  3:351195/1‑151 (MQ=255)
                                                              aCGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTgg                                                                                       <  2:509669/148‑1 (MQ=255)
                                                                ggcggTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAAtt                                                                                  <  1:242504/151‑1 (MQ=255)
                                                                   ggTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTAcc                                                                               >  2:406789/1‑151 (MQ=255)
                                                                   ggTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTAcc                                                                               >  2:365406/1‑151 (MQ=255)
                                                                    gTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCg                                                                              <  2:348060/151‑1 (MQ=255)
                                                                      tttATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCgg                                                                             <  2:1188901/150‑1 (MQ=255)
                                                                       ttATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCgggg                                                                           >  4:25272/1‑151 (MQ=255)
                                                                       ttATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCgggg                                                                           >  3:119540/1‑151 (MQ=255)
                                                                         aTCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGa                                                                          >  1:458776/1‑150 (MQ=255)
                                                                         aTCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGa                                                                          <  1:1180974/150‑1 (MQ=255)
                                                                         aTCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGa                                                                          >  2:829912/1‑150 (MQ=255)
                                                                            gTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTc                                                                      <  4:338998/151‑1 (MQ=255)
                                                                             tCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCg                                                                     >  4:34415/1‑151 (MQ=255)
                                                                              ccTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCg                                                                     >  2:860076/1‑150 (MQ=255)
                                                                               cTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGAt                                                                   >  2:774113/1‑151 (MQ=255)
                                                                                tGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATc                                                                  >  3:250244/1‑151 (MQ=255)
                                                                                 gTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGAt                                                                   <  1:1162605/149‑1 (MQ=255)
                                                                                 gTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGAt                                                                   <  3:309701/149‑1 (MQ=255)
                                                                                 gTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATcc                                                                 <  4:345039/151‑1 (MQ=255)
                                                                                   agagAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCg                                                                <  2:57687/150‑1 (MQ=255)
                                                                                   agagAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGc                                                               <  2:832256/151‑1 (MQ=255)
                                                                                    gagaAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCa                                                              <  2:1171880/151‑1 (MQ=255)
                                                                                     agaAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCAt                                                             <  1:438997/151‑1 (MQ=255)
                                                                                      gaAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCAt                                                             >  4:159627/1‑150 (MQ=255)
                                                                                       aaaTTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATAc                                                           >  1:236241/1‑151 (MQ=255)
                                                                                       aaaTTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATAc                                                           >  3:46993/1‑151 (MQ=255)
                                                                                        aaTTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACg                                                          <  1:909136/151‑1 (MQ=255)
                                                                                         aTTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACgg                                                         <  3:15764/151‑1 (MQ=255)
                                                                                           tATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACggg                                                        <  4:374342/150‑1 (MQ=255)
                                                                                            atgatgTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCt                                                      >  1:1168318/1‑151 (MQ=255)
                                                                                               atgTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGa                                                   <  2:137141/151‑1 (MQ=255)
                                                                                                 gTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGAt                                                  <  2:828528/150‑1 (MQ=255)
                                                                                                  tCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGAtt                                                 >  1:779472/1‑150 (MQ=255)
                                                                                                   ccaccaCGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTa                                                <  2:22886/150‑1 (MQ=255)
                                                                                                    caccacGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATc                                              >  1:1143231/1‑151 (MQ=255)
                                                                                                      ccacGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCg                                             >  1:1049579/1‑150 (MQ=255)
                                                                                                      ccacGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGa                                            <  2:253629/151‑1 (MQ=255)
                                                                                                      ccacGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGa                                            <  2:233026/151‑1 (MQ=255)
                                                                                                       cacGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGTTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGAt                                           <  1:949840/151‑1 (MQ=255)
                                                                                                       cacGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGAt                                           >  1:1270408/1‑151 (MQ=255)
                                                                                                       cacGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGAt                                           <  1:1233116/151‑1 (MQ=255)
                                                                                                        acGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGAtt                                          >  1:603524/1‑151 (MQ=255)
                                                                                                         cGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGAttt                                         <  1:981528/151‑1 (MQ=255)
                                                                                                           gTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTc                                        <  2:42834/151‑1 (MQ=255)
                                                                                                           gTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTc                                        <  4:382244/151‑1 (MQ=255)
                                                                                                            ttATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCg                                       >  4:150919/1‑151 (MQ=255)
                                                                                                                ttAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTg                                    <  3:140197/150‑1 (MQ=255)
                                                                                                                  aaaGATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGaa                                  <  2:554259/150‑1 (MQ=255)
                                                                                                                     gATTTCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTaa                              <  4:408445/151‑1 (MQ=255)
                                                                                                                        ttCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTaaaa                            <  1:179071/150‑1 (MQ=255)
                                                                                                                        ttCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTaaa                             >  3:453527/1‑149 (MQ=255)
                                                                                                                         tCACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAg                           >  1:363462/1‑150 (MQ=255)
                                                                                                                          cACCTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAgcg                         <  2:571464/151‑1 (MQ=255)
                                                                                                                            ccTTCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGCCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAgcggc                       >  1:130595/1‑151 (MQ=255)
                                                                                                                              ttCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGt                     >  3:355747/1‑151 (MQ=255)
                                                                                                                               tCGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGt                     <  4:380717/150‑1 (MQ=255)
                                                                                                                                cGAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGtt                    <  2:851361/150‑1 (MQ=255)
                                                                                                                                 gAAGCCGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGttt                   >  1:463675/1‑150 (MQ=255)
                                                                                                                                     ccGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAaacc              <  2:1123206/151‑1 (MQ=255)
                                                                                                                                      cGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAaacct             <  2:534145/151‑2 (MQ=255)
                                                                                                                                      cGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAaacc              >  3:170963/1‑150 (MQ=255)
                                                                                                                                      cGCTCACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAaacc              >  3:110145/1‑150 (MQ=255)
                                                                                                                                          cACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAaaccaacc          <  1:884161/150‑1 (MQ=255)
                                                                                                                                          cACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAaaccaacc          <  4:268313/150‑1 (MQ=255)
                                                                                                                                          cACCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCt         >  1:367570/1‑151 (MQ=255)
                                                                                                                                           aCCGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTa        <  2:1046373/151‑1 (MQ=255)
                                                                                                                                            ccGCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTAc       >  2:961319/1‑151 (MQ=255)
                                                                                                                                              gCTTACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACg      >  1:357438/1‑150 (MQ=255)
                                                                                                                                                 tACCGCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAgcg  >  2:573752/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CAGCGATTTTGCGCAACAAAGTCGTTTTAGATAATGCGAAAAAACAGCCTTTCCGGTACTCTACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCG  >  NC_000913/2892122‑2892416

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