Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,740,5270AT36.1% 95.2 / 105.7 127intergenic (+36/+35)hutF/PP_5037formiminoglutamate deiminase/lipocalin family lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (24/45);  new base T (17/22);  total (41/86)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.73e-02
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 136 total aligned reads displayed.

CAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTC  >  NC_002947/5740378‑5740676
                                                                                                                                                     |                                                                                                                                                     
cAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCTTGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCaa                                                                                                                                                       <  2:3582518/150‑1 (MQ=255)
cAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCaa                                                                                                                                                       <  1:3910673/150‑1 (MQ=255)
cAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCaa                                                                                                                                                       <  1:4178509/150‑1 (MQ=255)
cAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCaa                                                                                                                                                       <  1:4516919/150‑1 (MQ=255)
  ggTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGc                                                                                                                                                     <  1:2594207/150‑1 (MQ=255)
  ggTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGc                                                                                                                                                     <  1:2594197/150‑1 (MQ=255)
         gATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGaa                                                                                                                                              <  2:2064640/150‑1 (MQ=255)
         gATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGaa                                                                                                                                              <  2:2578778/150‑1 (MQ=255)
         gATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGaa                                                                                                                                              <  2:2325061/150‑1 (MQ=255)
         gATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGGGCTTGCACGCGaa                                                                                                                                              <  2:4960964/150‑1 (MQ=255)
         gATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGaa                                                                                                                                              <  2:1089048/150‑1 (MQ=39)
          aTGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAGAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAg                                                                                                                                             <  1:4929582/150‑1 (MQ=39)
          aTGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAg                                                                                                                                             <  2:703884/150‑1 (MQ=255)
          aTGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGGGCTTGCACGCGAAg                                                                                                                                             <  1:943267/150‑1 (MQ=255)
          aTGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAg                                                                                                                                             <  2:374744/150‑1 (MQ=39)
          aTGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAg                                                                                                                                             <  1:2311997/150‑1 (MQ=39)
          aTGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAg                                                                                                                                             <  1:2312066/150‑1 (MQ=39)
          aTGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAg                                                                                                                                             <  2:374747/150‑1 (MQ=39)
          aTGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAg                                                                                                                                             <  1:3597193/150‑1 (MQ=39)
             tgctGGTGACTGGGCATTGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGGGCTTGCACGCGAAgggc                                                                                                                                          <  1:104620/147‑2 (MQ=38)
              gATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAATCGGGCGGGAGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGtt                                                                                                                                         <  1:1136209/150‑1 (MQ=255)
              gATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAgggct                                                                                                                                         <  2:1590678/150‑3 (MQ=38)
                  gTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgcaa                                                                                                                                     <  1:84229/150‑1 (MQ=255)
                   tgcaTGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGACAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgcaac                                                                                                                                    <  1:403718/147‑1 (MQ=255)
                   tGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgcaac                                                                                                                                    <  2:4232353/150‑1 (MQ=255)
                   tGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgcaac                                                                                                                                    <  1:1476511/150‑1 (MQ=255)
                        gggCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGGTGCATCCCATCCCCTGCAAGCCCCCCAAAGGTTtcaaaccaac                                                                                                                               >  2:4808296/1‑150 (MQ=255)
                                     tgCGCCTGGGGCGGCTTGCCCTGTAATATTAACGGGGGGGGGGTTTTGCCAGGGTGTTGCGCCAGTTTTTGGGCTGAGGTGACGGGGTTTGCTTTGCAGCCATTCGGGTGAAAGCCCGCGATGGTTTGAAACGCAACCCCAAAAAACTCa                                                                                                                  <  2:4052339/150‑1 (MQ=255)
                                     tgCGCCGGGGGCGGCTTGCCCTGTACGATGACCGCGGGGGGGCGTTGGCCATGTTGTTGCGCCTGCTTTTGGGCTGCGGGTACGGGGGGTGGGTTGCAGCCCTGGGCTGGAACGCCGGCGAAGTGTTCCAACGCAACCCCAAAAAACTCa                                                                                                                  <  2:4051977/150‑1 (MQ=255)
                                        gcCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAg                                                                                                                                             >  2:4537789/1‑120 (MQ=255)
                                        gcCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAg                                                                                                                                             <  1:4537789/120‑1 (MQ=255)
                                         ccAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAgggc                                                                                                                                          >  2:1570719/1‑121 (MQ=255)
                                         ccAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAgggc                                                                                                                                          <  1:1570719/122‑2 (MQ=255)
                                            ggggCGGCATGCCCAGGAAGATGAACGCGGGCGGGCGTTTGCCACGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGACGGGTTGCAACGCAACCCCAAAAAACTCAGTGCAtt                                                                                                           <  2:1665344/150‑1 (MQ=255)
                                            ggggCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAggg                                                                                                                                           <  2:5277411/118‑1 (MQ=255)
                                            ggggCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAggg                                                                                                                                           >  1:5277411/1‑118 (MQ=255)
                                                 ggctggCCCCGTAAGTTGAACGGGGGGGGGCGTTTGCCATGGTGTTGCGCCAGCTTTGGGGGTGCGCGTACGGGGGTTGCGTGGCCGCCGTGCGGTGGCACGCCCGCGATGGTTTGCACCGCAACCCCAAAAAACTCAGTGCATTTTTAc                                                                                                      <  2:1919505/145‑1 (MQ=255)
                                                      gCCCTGGAAGATTAACGGGGGGGGGCGTGTGCCATGGTGTTGCGCCAGTTTTGGGGCTGAGCGCACGGGGGTGGCGTTGCAGCCCTGCGCGTGCCCGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTAttt                                                                                                 <  2:4552862/150‑1 (MQ=255)
                                                        cctggAAGATGAACGCGGGGGGGGGTTGGCCAGGGTGTTGCGCCAGCTTTTGGGCTGCGCGTACGGGGGTTGGGTGGCAGCCCTCCGCGTGCACGCCGGCGAGGGTTTGCACCGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGc                                                                                               <  2:1829852/147‑1 (MQ=255)
                                                             aaGATGAATGCGGGCGGGCTTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGCGCGTACGGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATc                                                                                          <  2:793017/150‑1 (MQ=255)
                                                                 ttaaaGCGGGCGGGGGTTTGCCATGTTGTTGCGCCTGCTTTTGGGCTGAGCTTACGGGGGTGGCTTTGCAGCCCTTGGCGTGCAAGCCCGCGATGGTTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGAc                                                                                      <  2:2442957/148‑1 (MQ=255)
                                                                           ggggcgTGGGCCACGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACGGGGTTTGCGTTGCAGCCCTTCGCGTCCAAGCCCGCTATGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGa                                                                            <  2:1105156/149‑1 (MQ=255)
                                                                           cgggcgTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGa                                                                            <  2:4677887/150‑1 (MQ=255)
                                                                            ggggttttGCCTGGTTGTGGCCCCCGTTTTTGGGGTGAGCGTCCGGGGGTTGCTTTCCACCCCTTCGGGCCAAAGCCCGCTAGGGTTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGAt                                                                           <  2:1096405/145‑1 (MQ=255)
                                                                                 tttGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGc                                                                      <  2:4807755/150‑1 (MQ=255)
                                                                                   tGCCATGGTGTGGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCGTTCGGGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCgg                                                                    <  2:5118579/150‑1 (MQ=255)
                                                                                      ctgggTGTTGGGCCTGTTTTTGGGGTGCGGGTCCGGGGGTTGGGTTGCAGCCCTGGGCGTGCACGCCCGGGATGGTTTGCACCGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTg                                                                 <  1:2095895/147‑1 (MQ=255)
                                                                                      cAAGGTGTTGCGCCAGCTTTTGGGCTTAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCATGCCCGCGATGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTg                                                                 <  2:2480375/150‑1 (MQ=255)
                                                                                      cAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAgg                                                                                                                                            >  2:1063600/1‑75 (MQ=255)
                                                                                      cAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAgg                                                                                                                                            <  1:1063309/75‑1 (MQ=255)
                                                                                      cAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAgg                                                                                                                                            >  2:1063309/1‑75 (MQ=255)
                                                                                          gtgtGGCGCCCTCTTTTGGGCTGTGGGTACGGGGGTTGCGGTGCAGCCCTTCGGGTGCACGCCGGCGATGGTTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGgtgt                                                             <  2:654833/150‑1 (MQ=255)
                                                                                          gggtTGCGCCTGTTTTGGGGCTGCGGGTACGGGGTTGGCTTTGCAGCCCTGCGGGTGCACGCCCGCGACGGTTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGgtgt                                                             <  1:3854667/148‑1 (MQ=255)
                                                                                            gtTGCGCCAGCTTTTGGGCTTAGCGTACCGGGGTTGGGTTGCGGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCa                                                           <  2:2954790/150‑1 (MQ=255)
                                                                                              tGCGCCTGCTTTTGGGCTGAGCGTACCGGGTTTGCGTTGCAGCCCTTCGCTTGCACGCCGGCGACGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATa                                                         <  2:2718299/150‑1 (MQ=255)
                                                                                                                 cgCGTACGGGGGTGGCGTTGCAGCCCTTCGCGTGCAAGCCCGCCATGGGTTGCAACGCAACCCCAAAGAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTg                                      <  2:2830376/149‑1 (MQ=255)
                                                                                                                 cgCGTACCGGGTTTGGGTTGCAGCCCTTCGCTTGAAAGCCCGCTACGGGTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTg                                      <  2:927030/149‑1 (MQ=255)
                                                                                                                 aGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCACGCCCGCGACGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTg                                      <  2:2830393/150‑1 (MQ=255)
                                                                                                                   cGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGATGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGa                                    <  2:3065919/150‑1 (MQ=255)
                                                                                                                        cGGGGTTGCGTTGCCGCCCTTCGCGGGCAAGCCCGCGACGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagca                               <  2:3342496/150‑1 (MQ=255)
                                                                                                                        cGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGACGGTTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagca                               <  1:4808296/150‑1 (MQ=255)
                                                                                                                            gttgcgtggcAGCCCTGCGCGTGCACGCCCGCGACGGGTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGtt                           <  2:958694/150‑1 (MQ=255)
                                                                                                                             ttgcgttgcAGCCCTTCGCGTGCAAGCCCGCGATGGGTTGCAACGCAACCCCAAACAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtc                          <  2:1169617/150‑1 (MQ=255)
                                                                                                                             ttgcgttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtc                          >  1:942378/1‑150 (MQ=255)
                                                                                                                               gctttgcAGCCTTTGGCTTGCAAGCCCGCGATGGGTTCAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtctc                        <  1:2548027/147‑1 (MQ=255)
                                                                                                                                 gttgcAACCCTTCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCAc                      >  2:5283118/1‑150 (MQ=38)
                                                                                                                                  ttgcAACCCTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGa                                                                                       <  1:107849/84‑1 (MQ=18)
                                                                                                                                      aacccTTCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTg                                                                                                <  1:3222497/69‑1 (MQ=18)
                                                                                                                                      aacccTTCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTg                                                                                                >  2:3222497/3‑71 (MQ=18)
                                                                                                                                      aacccTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGc                                                                   <  1:4514359/98‑1 (MQ=25)
                                                                                                                                      aacccTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTcc                 >  2:1239704/3‑150 (MQ=38)
                                                                                                                                      aacccTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTcc                 >  1:2451160/3‑150 (MQ=38)
                                                                                                                                      aacccTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTcc                 >  2:2673549/3‑150 (MQ=38)
                                                                                                                                      aGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTcc                 >  1:4137052/1‑150 (MQ=255)
                                                                                                                                      aGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTcc                 >  1:2543055/1‑150 (MQ=255)
                                                                                                                                       gCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTccc                >  1:5164001/1‑150 (MQ=255)
                                                                                                                                       gCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTccc                >  2:2942604/1‑150 (MQ=255)
                                                                                                                                        cccTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCa               >  2:1056838/1‑150 (MQ=255)
                                                                                                                                        cccTTCGCGTCCAAGCCCGCGAAGGGTTGCACCGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCa               <  1:5157660/150‑1 (MQ=255)
                                                                                                                                        cccTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCa               >  1:5185352/1‑150 (MQ=39)
                                                                                                                                         ccTTGGGGTGCACGCCGGCTTTGGTTTCCACCGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAg              <  1:3647861/150‑1 (MQ=255)
                                                                                                                                         ccTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtctc                        >  2:588782/1‑140 (MQ=255)
                                                                                                                                         ccTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtctc                        <  1:588782/140‑1 (MQ=255)
                                                                                                                                          cTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAgg             >  2:5085530/1‑150 (MQ=255)
                                                                                                                                          cTTCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAgg             >  1:1188616/1‑150 (MQ=255)
                                                                                                                                          cTTCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAgg             >  2:2723542/1‑150 (MQ=255)
                                                                                                                                          cTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAgg             >  1:376702/1‑150 (MQ=39)
                                                                                                                                          cTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAgg             >  2:1562623/1‑150 (MQ=39)
                                                                                                                                           ttCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGc            >  1:3661352/1‑150 (MQ=255)
                                                                                                                                           ttCGCGGGCTTGCACGCGGAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGc            >  1:3930021/1‑150 (MQ=255)
                                                                                                                                           ttCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGa                                    >  1:2643408/1‑126 (MQ=38)
                                                                                                                                           ttCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGc            >  1:1308059/1‑150 (MQ=39)
                                                                                                                                            tCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGcc           >  1:2196592/1‑150 (MQ=255)
                                                                                                                                              gcgTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCaa                                                                        >  1:1057958/1‑87 (MQ=255)
                                                                                                                                              gcgGGCACGCCCGCTACGGGTTCCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCgg         <  1:2972897/150‑1 (MQ=255)
                                                                                                                                                   cAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAgtg    >  1:1781013/1‑150 (MQ=255)
                                                                                                                                                    aaGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGa                                    <  1:3399558/117‑1 (MQ=255)
                                                                                                                                                    aaGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGa                                    >  2:3399558/1‑117 (MQ=255)
                                                                                                                                                    aaGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTAAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATAACCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAgtgt   >  1:3118101/1‑150 (MQ=255)
                                                                                                                                                     aGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTc  >  2:2386774/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTC  >  NC_002947/5740378‑5740676

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