Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,170,145 T→C pseudogene (23/475 nt) gatR ← pseudogene, repressor for gat operon; interrupted by IS3; split galactitol utilization operon repressor, fragment 2; split galactitol utilization operon repressor, interrupted

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,170,1450TC99.1% 428.7 / NA 112pseudogene (23/475 nt)gatRpseudogene, repressor for gat operon; interrupted by IS3; split galactitol utilization operon repressor, fragment 2; split galactitol utilization operon repressor, interrupted
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (60/50);  total (61/51)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 112 total aligned reads displayed.

ATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTG  >  NC_000913/2169997‑2170287
                                                                                                                                                    |                                                                                                                                              
aTCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCAc                                                                                                                                              >  2:1020103/1‑151 (MQ=255)
  cGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACt                                                                                                                                             <  4:1104846/150‑1 (MQ=255)
  cGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACt                                                                                                                                             <  4:1662609/150‑1 (MQ=255)
   gAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTc                                                                                                                                           <  3:267591/151‑1 (MQ=255)
   gAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTc                                                                                                                                           >  4:159804/1‑151 (MQ=255)
     gcgcGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCg                                                                                                                                         <  2:830359/151‑1 (MQ=255)
     gcgcGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCg                                                                                                                                         >  4:555996/1‑151 (MQ=255)
       gcgTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCt                                                                                                                                       <  3:475897/151‑1 (MQ=255)
       gcgTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCt                                                                                                                                       >  2:1305392/1‑151 (MQ=255)
        cgTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTg                                                                                                                                      >  4:521474/1‑151 (MQ=255)
            ttATGGCGGACCGTGCCACCACAGGCGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGt                                                                                                                                     >  2:779846/1‑148 (MQ=255)
             tATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCga                                                                                                                                  <  2:396701/150‑1 (MQ=255)
              aTGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCgaga                                                                                                                                <  4:1549345/151‑1 (MQ=255)
                ggCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGaa                                                                                                                               <  2:554332/150‑1 (MQ=255)
                ggCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGaa                                                                                                                               <  1:158421/150‑1 (MQ=255)
                 gCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAAt                                                                                                                              >  4:598413/1‑150 (MQ=255)
                  cGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATa                                                                                                                             >  1:209063/1‑150 (MQ=255)
                   ggACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAAc                                                                                                                           <  2:1310871/151‑1 (MQ=255)
                     aCCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACgg                                                                                                                         >  4:1508828/1‑151 (MQ=255)
                       cGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTc                                                                                                                       <  4:82098/151‑1 (MQ=255)
                        gTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCt                                                                                                                      >  4:351802/1‑151 (MQ=255)
                            caccacAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATc                                                                                                                  >  4:480851/1‑151 (MQ=255)
                             accacAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCt                                                                                                                 >  2:65892/1‑151 (MQ=255)
                             accacAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCt                                                                                                                 <  3:37374/151‑1 (MQ=255)
                                acaGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTAc                                                                                                              <  4:737949/151‑1 (MQ=255)
                                  aGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTAccc                                                                                                            >  3:99039/1‑151 (MQ=255)
                                   gacgacCATAGTAATGTCTTTATTTTCGGAAATGGCAAACGCTGCTGGGAGACTGTTGGTAATAACCGTGATATTTTTGTCAGACATCACTCCCCCCCCGATAAGCAATTTTGTCCCTCCCCTGTTGCCACTACCGGTCTGTTcccaaacc                                                                                                           <  2:1339897/151‑8 (MQ=255)
                                     tgacCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGc                                                                                                         >  2:1526954/2‑151 (MQ=255)
                                      gacCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGc                                                                                                         >  4:1290224/1‑150 (MQ=255)
                                      gacCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCa                                                                                                        <  3:174004/151‑1 (MQ=255)
                                       acCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCaa                                                                                                       <  2:1212038/151‑1 (MQ=255)
                                       acCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCa                                                                                                        >  4:256993/1‑150 (MQ=255)
                                         cAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAAt                                                                                                      >  2:1484089/1‑150 (MQ=255)
                                            agTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTg                                                                                                  <  1:125305/151‑1 (MQ=255)
                                               aaTGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGac                                                                                               >  4:1424663/1‑151 (MQ=255)
                                               aaTGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGac                                                                                               >  1:404652/1‑151 (MQ=255)
                                                aTGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGac                                                                                               >  1:165809/1‑150 (MQ=255)
                                                 tGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGacac                                                                                             >  4:892795/1‑151 (MQ=255)
                                                  gTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACAcg                                                                                            >  3:255962/1‑151 (MQ=255)
                                                    cTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACAcgcg                                                                                          >  3:410313/1‑151 (MQ=255)
                                                      ttATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACAcgcg                                                                                          <  1:83442/149‑1 (MQ=255)
                                                        aTTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCt                                                                                       >  4:1157634/1‑150 (MQ=255)
                                                        aTTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCt                                                                                       <  2:737470/150‑1 (MQ=255)
                                                        aTTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCt                                                                                       <  4:702301/150‑1 (MQ=255)
                                                        aTTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTa                                                                                      <  3:454365/151‑1 (MQ=255)
                                                        aTTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTa                                                                                      <  3:443876/151‑1 (MQ=255)
                                                        aTTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTa                                                                                      <  4:825949/151‑1 (MQ=255)
                                                           ttCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGt                                                                                   >  1:159559/1‑151 (MQ=255)
                                                           ttCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGt                                                                                   >  2:1410088/1‑151 (MQ=255)
                                                            tCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGt                                                                                   <  2:426838/150‑1 (MQ=255)
                                                             cGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTg                                                                                  >  4:189556/1‑150 (MQ=255)
                                                             cGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTg                                                                                  >  2:59588/1‑150 (MQ=255)
                                                                aaaGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGt                                                                               >  4:1280927/1‑150 (MQ=255)
                                                                 aaGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTa                                                                              >  4:1273324/1‑150 (MQ=255)
                                                                 aaGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTa                                                                              >  2:750665/1‑150 (MQ=255)
                                                                  aGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTaaa                                                                            <  3:192329/151‑1 (MQ=255)
                                                                     gCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAActc                                                                         <  2:833252/151‑1 (MQ=255)
                                                                      cAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAActct                                                                        >  1:66496/1‑151 (MQ=255)
                                                                       aaaCGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACtctc                                                                       >  2:183139/1‑151 (MQ=255)
                                                                          cGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCtcagt                                                                    >  4:1310464/1‑151 (MQ=255)
                                                                           gCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCtcagtc                                                                   >  1:69008/1‑151 (MQ=255)
                                                                           gCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCtcagtc                                                                   <  4:859758/151‑1 (MQ=255)
                                                                           gCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGAACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCtcag                                                                     >  3:43159/1‑149 (MQ=255)
                                                                            cTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCtcagtc                                                                   >  4:797348/1‑150 (MQ=255)
                                                                                 gggAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGt                                                             <  4:88785/151‑1 (MQ=255)
                                                                                 gggAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGt                                                             >  4:264361/1‑151 (MQ=255)
                                                                                 gggAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGt                                                             <  1:7576/151‑1 (MQ=255)
                                                                                          ttGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGCGACTCACAt                                                    >  2:695342/1‑151 (MQ=255)
                                                                                           tGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACAt                                                    <  4:1562534/150‑1 (MQ=255)
                                                                                            ggTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATg                                                   <  4:1605390/150‑1 (MQ=255)
                                                                                               aaTCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACa                                                <  2:295237/150‑1 (MQ=255)
                                                                                                 tCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACaaa                                              >  4:660020/1‑150 (MQ=255)
                                                                                                 tCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACaa                                               >  1:282647/1‑149 (MQ=255)
                                                                                                  cACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACaaaa                                             <  2:1682917/150‑1 (MQ=255)
                                                                                                  cACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACaaaa                                             >  4:779576/1‑150 (MQ=255)
                                                                                                    ccGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACa                                          <  2:723561/151‑1 (MQ=255)
                                                                                                    ccGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACa                                          <  3:203467/151‑1 (MQ=255)
                                                                                                     cGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAAc                                           >  2:1563178/1‑149 (MQ=255)
                                                                                                           atTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCa                                    >  4:1182705/1‑150 (MQ=255)
                                                                                                                tGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGt                              >  2:125368/1‑151 (MQ=255)
                                                                                                                tGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGt                              >  2:1444033/1‑151 (MQ=255)
                                                                                                                 ggCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGt                              >  2:674125/1‑150 (MQ=255)
                                                                                                                 ggCGGTCATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTa                             <  1:277117/151‑1 (MQ=255)
                                                                                                                       cATTAATCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTaaaaaa                        <  4:156370/150‑1 (MQ=255)
                                                                                                                           aaTCCTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCg                   >  4:573857/1‑151 (MQ=255)
                                                                                                                           aaTCCTTCGGCGATAAGCATTGCCGCACTTCCGCGGTCGAGAATAACGGTCTGATCTTACCCAGCAGTATTGGACACGCGGATAAGTGAGTAACCTTTCAGTCGGAGGTGACTCACATGACAAAAACAGTATCAACCATTAAAAAccccg                    <  2:1008958/150‑2 (MQ=18)
                                                                                                                              ccTTCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGATACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTa                 <  2:159532/150‑1 (MQ=255)
                                                                                                                                 tCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTaaa                           <  4:852754/137‑1 (MQ=255)
                                                                                                                                 tCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTaaa                           >  2:851639/1‑137 (MQ=255)
                                                                                                                                 tCGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAAc              <  2:102676/150‑1 (MQ=255)
                                                                                                                                  cGGCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAAc              >  4:1638720/1‑149 (MQ=255)
                                                                                                                                   ggCGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGc           >  4:470057/1‑151 (MQ=255)
                                                                                                                                     cGATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCAt         <  4:66015/151‑1 (MQ=255)
                                                                                                                                      gATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCAtt        >  4:1533231/1‑151 (MQ=255)
                                                                                                                                      gATAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCAtt        >  2:455941/1‑151 (MQ=255)
                                                                                                                                        tAAGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCg      >  4:121035/1‑151 (MQ=255)
                                                                                                                                         aaGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCg      <  1:188770/150‑1 (MQ=255)
                                                                                                                                         aaGCATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGc     <  4:1637826/151‑1 (MQ=255)
                                                                                                                                            cATTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTg  >  1:52648/1‑151 (MQ=255)
                                                                                                                                             aTTGTCGCACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTg  >  3:165363/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
ATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTG  >  NC_000913/2169997‑2170287

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