Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,343,645 A→G 6.6% I175I (ATT→ATC yfaL ← putative autotransporter adhesin YfaL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,343,6450AG 6.6% 230.7 / 4.3 77I175I (ATT→ATCyfaLputative autotransporter adhesin YfaL
Reads supporting (aligned to +/- strand):  ref base A (33/38);  new base G (3/2);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCGTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGAAAAAATAACATCAGT  >  NC_000913/2343496‑2343794
                                                                                                                                                     |                                                                                                                                                     
aTCGTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATaa                                                                                                                                                       >  1:1356469/1‑150 (MQ=255)
  cGTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAAt                                                                                                                                                     <  1:963243/150‑1 (MQ=255)
   gTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAGAtt                                                                                                                                                    >  1:6042/1‑150 (MQ=255)
      taACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTgcg                                                                                                                                                 <  2:1046341/150‑1 (MQ=255)
        aCCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGcc                                                                                                                                               <  2:816074/150‑1 (MQ=255)
        aCCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGcc                                                                                                                                               <  1:588355/150‑1 (MQ=255)
        aCCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGAATAAATTGCGcc                                                                                                                                               >  1:820513/1‑150 (MQ=255)
           tGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCAcc                                                                                                                                            >  1:839854/1‑150 (MQ=255)
           tGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCAcc                                                                                                                                            >  1:1395157/1‑150 (MQ=255)
            gAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCt                                                                                                                                           >  2:787443/1‑150 (MQ=255)
                 aTGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGcc                                                                                                                                      <  1:398076/150‑1 (MQ=255)
                  tgcttgctGTCGCTATTATTATCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCca                                                                                                                                     <  2:1395157/150‑1 (MQ=255)
                   gcttgctGTCGCTATTATTATCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCcat                                                                                                                                    >  1:1372987/1‑150 (MQ=255)
                    cttgctGTCGCTATTATTATCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCcatc                                                                                                                                   >  2:280658/1‑150 (MQ=255)
                     ttgctGTCGCTATTATTATCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCcatca                                                                                                                                  <  2:113673/150‑1 (MQ=255)
                      tgctGTCGCTATTATTATCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCcatcat                                                                                                                                 >  2:196164/1‑150 (MQ=255)
                             tcttatattatCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGAt                                                                                                                          <  1:8698/145‑1 (MQ=255)
                              cTATTATTATCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATg                                                                                                                         >  2:328032/1‑150 (MQ=255)
                               tattattatCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGt                                                                                                                        <  2:1318859/150‑1 (MQ=255)
                                      atCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAGATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCg                                                                                                                 <  2:6042/150‑1 (MQ=255)
                                       tCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGa                                                                                                                <  2:1136947/150‑1 (MQ=255)
                                         ggTAACATCGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGaaa                                                                                                              >  1:864405/1‑150 (MQ=255)
                                          gTAACATCGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAAc                                                                                                             <  1:196164/150‑1 (MQ=255)
                                                 cGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGt                                                                                                      <  1:784586/150‑1 (MQ=255)
                                                     tGCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGt                                                                                                  <  2:152009/150‑1 (MQ=255)
                                                      gCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTg                                                                                                 <  2:1170046/150‑1 (MQ=255)
                                                      gCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTg                                                                                                 <  1:553022/150‑1 (MQ=255)
                                                      gCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTg                                                                                                 <  2:1372987/150‑1 (MQ=255)
                                                         ccGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGaca                                                                                              <  2:432311/150‑1 (MQ=255)
                                                         ccGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGaca                                                                                              >  2:1070615/1‑150 (MQ=255)
                                                         ccGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGaca                                                                                              <  2:1356469/150‑1 (MQ=255)
                                                            ccatcgccatcACTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGt                                                                                           <  1:1087892/150‑1 (MQ=255)
                                                                cgccatcACTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAAt                                                                                       >  1:1044372/1‑150 (MQ=255)
                                                                 gccatcACTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTTCTAAACATTGCGTTAGTGACACGTAAAta                                                                                      >  1:613683/1‑149 (MQ=255)
                                                                   catcACTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGa                                                                                    >  1:1186236/1‑150 (MQ=255)
                                                                      cACTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACgg                                                                                 >  2:1081017/1‑150 (MQ=255)
                                                                      cACTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACgg                                                                                 >  2:925204/1‑150 (MQ=255)
                                                                      cACTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACgg                                                                                 >  1:319643/1‑150 (MQ=255)
                                                                       aCTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATCGCGTTAGTGACACGTAAATCGACGGc                                                                                <  1:809569/150‑1 (MQ=255)
                                                                          gTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCAcc                                                                             >  2:1378801/1‑150 (MQ=255)
                                                                            aaCTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAg                                                                           <  2:1186236/150‑1 (MQ=255)
                                                                            aaCTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAg                                                                           <  1:343335/150‑1 (MQ=255)
                                                                             aCTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGt                                                                          <  1:683377/150‑1 (MQ=255)
                                                                             aCTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGt                                                                          >  2:17590/1‑150 (MQ=255)
                                                                             aCTTGTTGATGTATATGCATTGTTGTTATCAAAAATAACATCACTTAAATAAACCTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGATGGCACCAGt                                                                          >  2:27345/1‑150 (MQ=255)
                                                                                tgttgATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAGATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATc                                                                       >  2:389074/1‑150 (MQ=255)
                                                                                  ttgATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGt                                                                     <  2:319643/150‑1 (MQ=255)
                                                                                   tgATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGtt                                                                    <  1:1081017/150‑1 (MQ=255)
                                                                                   tgATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGtt                                                                    <  2:828502/150‑1 (MQ=255)
                                                                                   tgATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGtt                                                                    <  2:823507/150‑1 (MQ=255)
                                                                                   tgACGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAAGGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGtt                                                                    <  1:418322/150‑1 (MQ=255)
                                                                                    gATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTa                                                                   <  1:638979/150‑1 (MQ=255)
                                                                                       gTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCTTTATCGATTTTTTTTCGAAACATGGCGTTAGTTACACGTTAATAGACGGCACCAGTATCGTTAGTa                                                                <  2:798865/150‑1 (MQ=255)
                                                                                        tataTGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTAc                                                               >  1:525529/1‑150 (MQ=255)
                                                                                        tataTGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTAc                                                               >  2:360916/1‑150 (MQ=255)
                                                                                        tataTGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTAc                                                               >  2:196548/1‑150 (MQ=255)
                                                                                         atatGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTAcc                                                              <  1:1246813/150‑1 (MQ=255)
                                                                                            tGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCaga                                                           >  2:79787/1‑150 (MQ=255)
                                                                                               cTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCagaaga                                                        <  1:280717/150‑1 (MQ=255)
                                                                                               cTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCagaaga                                                        >  1:1125907/1‑150 (MQ=255)
                                                                                                tGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAgaagaa                                                       <  2:1044372/150‑1 (MQ=255)
                                                                                                 ggTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAAt                                                      >  1:154953/1‑150 (MQ=255)
                                                                                                         tCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTgcg                                              >  2:483018/1‑150 (MQ=255)
                                                                                                                tAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCAt                                       <  2:1038212/150‑1 (MQ=255)
                                                                                                                  aCATCACTTAAATAAACCTCATTATTAATGGTATAAATTGCCCCACCTTTGCCATCATTAGCGATTTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAg                                     <  2:151443/150‑1 (MQ=255)
                                                                                                                    cgatcTTAAATAAACGTCATTATTAATGGTATATTTTTCGCCACCTTTGCCATCATTATCGATGTTATTTCGAAACATTGCGTTAGTGACACGTAAATCGACGGCAACATTATCGTTATTACTATTTGAATAGATTGCGCCACCATAgcc                                   <  2:1177855/146‑1 (MQ=255)
                                                                                                                      cACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAgccgc                                 >  2:136044/1‑150 (MQ=255)
                                                                                                                      cACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAgccgc                                 >  1:744818/1‑150 (MQ=255)
                                                                                                                       aCTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAgccgcc                                <  2:606718/150‑1 (MQ=255)
                                                                                                                       aCTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAgccgcc                                >  1:606718/1‑150 (MQ=255)
                                                                                                                         ttAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTg                              >  1:578114/1‑150 (MQ=255)
                                                                                                                          taaataaaCGTCATTATTAATGGTATAGATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAgccagatcg                             >  1:1177855/1‑144 (MQ=255)
                                                                                                                                   gTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTAcc                    >  1:1286798/1‑150 (MQ=255)
                                                                                                                                            aaTGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGAAAAAAt           <  1:27345/150‑1 (MQ=255)
                                                                                                                                             aTGGTATAGATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGAAAAAATa          <  1:389074/150‑1 (MQ=255)
                                                                                                                                             aTGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGAAAAAATa          >  2:983284/1‑150 (MQ=255)
                                                                                                                                                 tataAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGAAAAAATAACAt      >  2:836953/1‑150 (MQ=255)
                                                                                                                                                  ataAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGAAAAAATAACATc     <  1:787443/150‑1 (MQ=255)
                                                                                                                                                     aaTTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGAAAAAATAACATCAGt  <  1:925204/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATCGTATAACCTGAAGGATGCTTGCTGTCGCTATTATTATCGGTAACATCGATTGCCCCGCCATCGCCATCACTGTAACTTGTTGATGTATATGCCTGGTTGTTATCAAAAATAACATCACTTAAATAAACGTCATTATTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGAAAAAATAACATCAGT  >  NC_000913/2343496‑2343794

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