Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,377,792 A→T 6.0% intergenic (‑70/+19) blc ← / ← ampC outer membrane lipoprotein cell division and growth lipocalin/penicillin‑binding protein; beta‑lactamase, intrinsically weak

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,377,7920AT 6.0% 513.6 / 15.9 146intergenic (‑70/+19)blc/ampCouter membrane lipoprotein cell division and growth lipocalin/penicillin‑binding protein; beta‑lactamase, intrinsically weak
Reads supporting (aligned to +/- strand):  ref base A (64/73);  new base T (3/6);  total (67/79)
Fisher's exact test for biased strand distribution p-value = 5.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 153 total aligned reads displayed.

CACGCGGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTGGAATAAACGCGACATAGCTACCAAATCCGCCG  >  NC_000913/4377644‑4377940
                                                                                                                                                    |                                                                                                                                                    
cacgcGGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGaaaa                                                                                                                                                      <  4:1464849/149‑1 (MQ=255)
   gcggcggcgTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTc                                                                                                                                                   >  4:1132176/1‑149 (MQ=255)
       cggcgTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAc                                                                                                                                               <  2:258511/149‑1 (MQ=255)
         gcgTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAcct                                                                                                                                             <  3:423050/149‑2 (MQ=255)
         gcgTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAcct                                                                                                                                             <  4:504352/149‑2 (MQ=255)
               gagaACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGa                                                                                                                                       >  1:1634496/1‑149 (MQ=255)
                  aaCTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAAtt                                                                                                                                    >  2:1687745/1‑149 (MQ=255)
                   ccTGCGGGCAACGACCAGAAATGCAGTTGTCGCTGCGGTAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACGGGTTTTTTTCTAAGGGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAAttt                                                                                                                                   <  2:1301134/148‑1 (MQ=255)
                    cTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAAtttt                                                                                                                                  >  2:1714486/1‑149 (MQ=255)
                    cTGCAGGCAACGACCAGAAATGCAGCTGTAGCTCCGCCAACGAGAGGGAGCCGGCGCAGAAATGTTTCCTCACTGGGTTTTTTCTAAGTGTAGATGACAGGAAGGAAAAGCGGAGAAACGGTCCGAAAATTCGGACCCCAGGGAAtttt                                                                                                                                  >  2:1574889/1‑149 (MQ=255)
                      gCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTAc                                                                                                                                <  1:857828/149‑1 (MQ=255)
                       cAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACt                                                                                                                               <  3:1328425/149‑1 (MQ=255)
                          gCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTa                                                                                                                            >  3:1025443/1‑149 (MQ=255)
                            aaCGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTaga                                                                                                                          <  1:368272/149‑1 (MQ=255)
                                  cAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTa                                                                                                                    >  3:322983/1‑149 (MQ=255)
                                   aGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTaa                                                                                                                   <  4:631088/149‑1 (MQ=255)
                                   aGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTaa                                                                                                                   <  4:1147793/149‑1 (MQ=255)
                                   aGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTaa                                                                                                                   <  4:187651/149‑1 (MQ=255)
                                    gAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAg                                                                                                                  >  2:1709056/1‑149 (MQ=255)
                                       aTGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAAt                                                                                                               <  3:1389451/149‑1 (MQ=255)
                                       aTGAAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTACTTACTGGTTTATTTCTAAGTGTAGATGACAGCAAGGAAAAGCGAAGAAAAGGTCCGAAGATTCGGACCCGATGGACTTTTACTGTAGAGCGTTAAGAAt                                                                                                               <  4:759066/149‑1 (MQ=255)
                                        tGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATc                                                                                                              <  3:458111/149‑1 (MQ=255)
                                         gCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCt                                                                                                             >  4:1279064/1‑149 (MQ=255)
                                           aGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTc                                                                                                                                                   <  1:1226759/109‑1 (MQ=255)
                                           aGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTc                                                                                                                                                   >  2:1226759/1‑109 (MQ=255)
                                           aGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTaa                                                                                                                   <  3:402923/141‑1 (MQ=255)
                                           aGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGc                                                                                                           >  3:1626835/1‑149 (MQ=255)
                                            gCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGcc                                                                                                          >  4:616075/1‑149 (MQ=255)
                                                tCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCAGAGGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAggc                                                                                                      >  4:797118/1‑149 (MQ=255)
                                                    tGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGgcggcg                                                                                                  >  3:1104411/1‑149 (MQ=255)
                                                        gCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGa                                                                                              <  3:464928/149‑1 (MQ=255)
                                                         cAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGAc                                                                                             >  1:984725/1‑149 (MQ=255)
                                                          aaCGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCgg                                                                                                                                                 <  3:544567/96‑1 (MQ=255)
                                                            cGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTa                                                                                                                    >  2:1643667/1‑123 (MQ=255)
                                                            cGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTa                                                                                                                    <  1:1643667/123‑1 (MQ=255)
                                                             gagaGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTGGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACtctc                                                                                         <  1:796078/149‑1 (MQ=255)
                                                                  ggAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTgg                                                                                    >  3:641989/1‑149 (MQ=255)
                                                                  ggAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTgg                                                                                    >  2:1127059/1‑149 (MQ=255)
                                                                     gCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGAtt                                                                                 >  4:244838/1‑149 (MQ=255)
                                                                     gCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGAtt                                                                                 <  1:275257/149‑1 (MQ=255)
                                                                     gCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGAtt                                                                                 <  1:101490/149‑1 (MQ=255)
                                                                      cAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTg                                                                                <  1:581889/149‑1 (MQ=255)
                                                                      cAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTg                                                                                <  1:1709056/149‑1 (MQ=255)
                                                                       aGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGAccc                                                                                                                                             >  3:652033/1‑87 (MQ=255)
                                                                       aGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGAccc                                                                                                                                             >  4:567675/1‑87 (MQ=255)
                                                                       aGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGAccc                                                                                                                                             <  3:567675/87‑1 (MQ=255)
                                                                       aGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTgg                                                                               >  4:200384/1‑149 (MQ=255)
                                                                         gcgcATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCt                                                                                                             <  4:1552486/117‑1 (MQ=255)
                                                                            cATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAg                                                                          >  4:735092/1‑149 (MQ=255)
                                                                               aaaTGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATa                                                                           >  1:131250/1‑145 (MQ=255)
                                                                               aaaTGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATa                                                                           <  2:131250/145‑1 (MQ=255)
                                                                                aaTGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGtttt                                                                      <  4:641989/149‑1 (MQ=255)
                                                                                aaTGTTTCCTTACTGGCTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGtttt                                                                      <  3:1069460/149‑1 (MQ=255)
                                                                                 aTGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGttttt                                                                     <  2:676213/149‑1 (MQ=255)
                                                                                    tttCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgtt                                                                  >  4:1519520/1‑149 (MQ=255)
                                                                                    tttCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgtt                                                                  <  3:1646685/149‑1 (MQ=255)
                                                                                    tttCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgtt                                                                  <  2:1210643/149‑1 (MQ=255)
                                                                                       ccTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGc                                                               <  3:1550464/149‑1 (MQ=255)
                                                                                        cTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTgg                                                                               >  1:1711812/1‑132 (MQ=255)
                                                                                        cTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGcc                                                              <  3:848006/149‑1 (MQ=255)
                                                                                            cTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTaa                                                                                                                   >  1:1425339/1‑92 (MQ=255)
                                                                                            cTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTaa                                                                                                                   <  2:1425339/92‑1 (MQ=255)
                                                                                             tGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGc                                                                                                           >  4:838547/1‑99 (MQ=255)
                                                                                             tGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCAt                                                         <  2:546793/149‑1 (MQ=255)
                                                                                             tGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCAt                                                         >  3:999976/1‑149 (MQ=255)
                                                                                              ggTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATc                                                        >  1:1009967/1‑149 (MQ=255)
                                                                                               gTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgttt                                                                 <  2:1399050/139‑1 (MQ=255)
                                                                                                tttttttCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTc                                                                                                                                                   >  3:1559745/1‑56 (MQ=255)
                                                                                                tttttttCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGActc                                                                                           >  3:811833/1‑112 (MQ=255)
                                                                                                tttttttCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGttttt                                                                     <  2:1170459/134‑1 (MQ=255)
                                                                                                tttttttCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGACAGGCGGCGTCGACTCTCGCTGGATTGGGATAGttttt                                                                     >  1:1170459/1‑134 (MQ=255)
                                                                                                   ttttCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGAt                                                   <  2:1077941/149‑1 (MQ=255)
                                                                                                      tCTAAGTGTAGATGACAGCAAGGAAAAGGGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGTGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATAca                                                <  3:797118/149‑2 (MQ=255)
                                                                                                                gATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCtttt                                      <  2:1026260/149‑1 (MQ=255)
                                                                                                                      aGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGc                                                                                       >  2:831329/1‑94 (MQ=255)
                                                                                                                      aGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTGGa                                >  4:166111/1‑149 (MQ=255)
                                                                                                                      aGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTGGa                                >  2:863610/1‑149 (MQ=255)
                                                                                                                       gCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCtttt                                      >  1:312164/1‑142 (MQ=255)
                                                                                                                         aaGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCt                                   >  4:497752/1‑143 (MQ=255)
                                                                                                                           ggAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATTGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTGGAATaaa                           >  3:1391366/1‑149 (MQ=255)
                                                                                                                            gAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTaa                                                                                                                   <  3:1173866/60‑1 (MQ=255)
                                                                                                                            gAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTGGAATaaaa                          <  3:490934/149‑2 (MQ=255)
                                                                                                                            gAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTGGAATAAAc                          <  1:127030/149‑1 (MQ=255)
                                                                                                                              aaaGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTGGAATAAAcgc                        <  2:655117/149‑1 (MQ=255)
                                                                                                                                 gCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTGGAATAAACGCGAc                     >  2:1627025/1‑149 (MQ=255)
                                                                                                                                      gaAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTGGAATAAACGCGACATAGc                >  4:672155/1‑149 (MQ=255)
                                                                                                                                      gaAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTGGAATAAACGCGACATAGc                >  2:119778/1‑149 (MQ=255)
                                                                                                                                        tcgggTCCGAATTTTCGGACGCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTATTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTGGAATAAACGCGACATAGCTa              >  4:1393211/4‑149 (MQ=255)
                                                                                                                                        tcgggTCCGAATTTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTGGAATAAACGCGACATAGCTa              >  4:1241005/4‑149 (MQ=255)
                                                                                                                                        aaaGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTGGAATAAACGCGACATAGCTa              <  2:1053999/149‑1 (MQ=255)
                                                                                                                                         aaGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTGGAATAAACGCGACATAGCTAc             >  3:1493008/1‑149 (MQ=255)
                                                                                                                                          gggTCCGAATTTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTGGAATAAACGCGACAt                   >  2:987731/2‑142 (MQ=255)
                                                                                                                                               cGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTGGAATAAACGCGACATAGCTACCAAATc       >  2:953058/1‑149 (MQ=255)
                                                                                                                                               cGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTGGAATAAACGCGACATAGCTACCAAATc       >  3:1367055/1‑149 (MQ=255)
                                                                                                                                                gAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATa                                                  <  3:1663040/105‑1 (MQ=255)
                                                                                                                                                 aaaaTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAggcggc                                                                                                   >  4:229502/1‑55 (MQ=255)
                                                                                                                                                 aaaaTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAggcggc                                                                                                   <  3:229502/55‑1 (MQ=255)
                                                                                                                                                 aaaaTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTGGAATAAACGCGACATAGCTACCAAATccg     >  1:1491065/1‑149 (MQ=255)
                                                                                                                                                  aaaTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTGGAATAAACGCGACATAGCTACCAAATccgc    <  2:659579/149‑1 (MQ=255)
                                                                                                                                                   aaTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTGGAATAAACGCGACATAGCTACCAAATccgcc   >  2:749246/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CACGCGGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTGGAATAAACGCGACATAGCTACCAAATCCGCCG  >  NC_000913/4377644‑4377940

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