Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 752,647 T→C 7.5% D50G (GAT→GGT)  ybgD ← putative fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913752,6470TC 7.5% 212.4 / 5.9 80D50G (GAT→GGT) ybgDputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base T (23/50);  new base C (3/3);  total (27/53)
Fisher's exact test for biased strand distribution p-value = 3.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.47e-01

TACGCCAACTTTGGATACCGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAACCTCCCCCTTAAAGTGAATTTCGCCAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGTCCTTTAAAC  >  NC_000913/752499‑752793
                                                                                                                                                    |                                                                                                                                                  
tACGCCAACTTTGGATACCGGCATTCCGCTACCGTTGTCAGAAGCACGCAGATCACAGTTGATCAAGCCAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTAt                                                                                                                                                    <  2:760982/149‑1 (MQ=255)
  cGCCAACTTTGGATACCGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATGACAGTTGATCAAGCGAATGTCGACGGGCACGTTATTGCAATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCa                                                                                                                                                  <  1:907809/149‑1 (MQ=255)
  cGCCAACTTTGGATACCGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCa                                                                                                                                                  <  1:405849/149‑1 (MQ=255)
   gCCAACTTTGGATACCGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCaa                                                                                                                                                 >  2:318429/1‑149 (MQ=255)
   gCCAACTTTGGATACCGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCaa                                                                                                                                                 >  2:1105186/1‑149 (MQ=255)
   gCCAACTTTGGATACCGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGCTTATATGGGTTGTCGTGACTTGTCCAAGAGCTATGTGTTTATCaa                                                                                                                                                 >  1:168737/1‑149 (MQ=255)
    ccAACTTTGGATACCGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCGCCCGGTTTATATGGCTTGTCGTGACTTGTCCAAGATCTATGTTTTATTCAat                                                                                                                                                >  2:41927/1‑149 (MQ=255)
          ttGGATACCGGCATTCCGCTACCGTTGTCAGGAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATAtcttc                                                                                                                                          <  1:180785/149‑1 (MQ=255)
           tGGATACCGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATAtcttct                                                                                                                                         <  1:535724/149‑1 (MQ=255)
            ggATACCGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTg                                                                                                                                        <  1:773616/149‑1 (MQ=255)
            ggATACCGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTg                                                                                                                                        <  1:67559/149‑1 (MQ=255)
            ggATACCGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTg                                                                                                                                        <  2:486581/149‑1 (MQ=255)
            ggATACCGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTg                                                                                                                                        <  1:130205/149‑1 (MQ=255)
             gATACCGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTgg                                                                                                                                       >  2:150922/1‑149 (MQ=255)
                  cGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGaa                                                                                                                                  >  1:206647/1‑149 (MQ=255)
                  cGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGaa                                                                                                                                  >  1:422002/1‑149 (MQ=255)
                  cGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGaa                                                                                                                                  >  1:1181618/1‑149 (MQ=255)
                   ggCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAAt                                                                                                                                 <  1:1290157/149‑1 (MQ=255)
                   ggCATTCCGCTACCGGTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGCTCTATGTGTTTATCAATATCTTCTGGATGAAt                                                                                                                                 <  2:1251766/149‑1 (MQ=255)
                    gCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTAGCAATATCTTCTGGATGAAtt                                                                                                                                >  2:919731/1‑149 (MQ=255)
                       ttCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTca                                                                                                                             >  1:1109420/1‑149 (MQ=255)
                        tCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGACGAATGTcac                                                                                                                            <  2:1122364/149‑1 (MQ=255)
                         ccGCTACCGTTTTCAGAAGCAGGCCGATCACAGTTGATCAAGCTAATGTCGACGGCCACTTTATTGCTATGACGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTcaca                                                                                                                           <  2:823137/149‑1 (MQ=255)
                          cGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTCACaa                                                                                                                          >  1:709514/1‑149 (MQ=255)
                           gCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTCACAAg                                                                                                                         <  2:773437/149‑1 (MQ=255)
                           gCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTCACAAg                                                                                                                         <  1:687732/149‑1 (MQ=255)
                           gCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTCACAAg                                                                                                                         <  1:9418/149‑1 (MQ=255)
                            cTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTCACAAgg                                                                                                                        <  2:1218527/149‑1 (MQ=255)
                            cTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTCACAAgg                                                                                                                        <  1:771358/149‑1 (MQ=255)
                                  tttgcAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTCACAAGGTGCTTc                                                                                                                  <  2:177853/145‑1 (MQ=255)
                                    gTCAGAAGCAGGCAGATTACAGTTGATCAAGCGAATTTCGACTGCCACTTGATTGTTATGATGCTCCCGTTTTATATGTGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCaa                                                                                                                <  2:577154/149‑1 (MQ=255)
                                    gTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCaa                                                                                                                <  2:402732/149‑1 (MQ=255)
                                    gTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCaa                                                                                                                <  1:318429/149‑1 (MQ=255)
                                     tCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTCACAGGTTGCTTCatt                                                                                                               >  1:925319/1‑147 (MQ=255)
                                     tCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAAt                                                                                                               <  1:1212117/149‑1 (MQ=255)
                                      cAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTAATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTCTCAATATCTTCTGGATGAATTTCACCAGGTGCTTCAATa                                                                                                              <  2:793111/149‑1 (MQ=255)
                                        gAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAAc                                                                                                            <  1:473271/149‑1 (MQ=255)
                                         aaGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAAcc                                                                                                           >  2:152558/1‑149 (MQ=255)
                                                   aTCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTACCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAACCTCCCCCTTaa                                                                                                 >  2:219341/1‑149 (MQ=255)
                                                   aTCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGGTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTACCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAACCTCCCCCTTaa                                                                                                 >  2:257553/1‑149 (MQ=255)
                                                   aTCACAGTTGATCAAGCGAATGTCGACGGCCACTTGATTGCTATGATGCTCCCGGGTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTACCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAACCTCCCCCTTaa                                                                                                 >  2:955145/1‑149 (MQ=255)
                                                    tCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAACCTCCCCCTTaaa                                                                                                >  2:377241/1‑149 (MQ=255)
                                                    gcacaGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAACCTCCCCCTTaaa                                                                                                <  1:706/148‑1 (MQ=255)
                                                     cacaGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAACCTCCCCCTTAAAg                                                                                               <  2:1226525/149‑1 (MQ=255)
                                                     cacaGTTAATCAAGCGAATGTCGACGGCAACTCTATTGCTACGATGCACCCGGTTTATATGGCTCGTCGTGACTTGCCCAAGATCTATGCTTTTATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAACCCCCCCCTTAAAg                                                                                               <  1:879036/149‑1 (MQ=255)
                                                     aataGTTGATCACGCTTATCGCGACGGCCACTTGATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAACCTCCCCCTTAAAg                                                                                               <  2:1220028/146‑1 (MQ=255)
                                                      acaGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAACCTCCCCCTTAAAGt                                                                                              <  2:752162/149‑1 (MQ=255)
                                                        cgctgatgaATCTAATGTCCACAGCCACTTTATTGCTATGATGCCCCCGGTTTATTTGGGTTGATGTGGCTTTTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATACCCCGCCCCTTATAGtca                                                                                            <  2:968361/146‑3 (MQ=255)
                                                        aGTTGATCAAGCTAATGTCGACGGCCACTTTATTGCTATGCTGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAACCTCCCCCTTAAAGTGa                                                                                            <  2:935405/149‑1 (MQ=255)
                                                             aTCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAACCTCCCCCTTAAAGTGAATTTc                                                                                       <  2:78156/149‑1 (MQ=255)
                                                              tCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAACCTCCCCCTTAAAGTGAATTTCg                                                                                      >  2:346144/1‑149 (MQ=255)
                                                                aaGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGCCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAACCTCCCCCTTAAAGTGAATTTAGcc                                                                                    >  1:178604/1‑149 (MQ=255)
                                                                 aGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAACCTCCCCCTTAAAGTGAATTTCGCCa                                                                                   <  1:65421/149‑1 (MQ=255)
                                                                   cGAATGTCGACGGCCAATTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAACCTCCCCCTTAGACTGAATTTCGCCAGa                                                                                 >  1:813462/1‑149 (MQ=255)
                                                                    gAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAACCTCCCCCTTAAAGTGAATTTCGCCAGaa                                                                                <  2:883221/149‑1 (MQ=255)
                                                                    gAATGTCGACGGCCACTTTATTGCTAGGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAACCTCCCCCTTAAAGTGAATTTCGCCAGaa                                                                                <  2:626269/149‑1 (MQ=255)
                                                                     aaTGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAACCTCCCCCTTAAAGTGAATTTCGCCAGAAc                                                                               <  2:918714/149‑1 (MQ=255)
                                                                         tCGGCGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCCATATCCGCTGGATGAATTTCACAAGGTGCTTCAATAACCTCCCCCTTAAAGTGAATTTCGCCAGAACCTtc                                                                           <  1:596951/149‑1 (MQ=255)
                                                                               gCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAACCTCCCCCTTAAAGTGAATTTCGCCAGAACCTtccacacc                                                                     <  2:665430/149‑7 (MQ=255)
                                                                                  aCTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAACCTCCCCCTTAAAGTGAATTTCGCCAGAACCTTCACCagcagc                                                                  >  2:317320/1‑149 (MQ=255)
                                                                                     tgtttGCTATGATGCCTCCGGGTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAACCTCCCCCTTAAAGTGAATTTCGCCAGAACCTTCATCAGCAGCATa                                                               <  2:3970/146‑1 (MQ=255)
                                                                                           cTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAACCTCCCCCTTAAAGTGAATTTCGCCAGAACCTTCATCAGCAGCATAAACAgg                                                         <  1:1105186/149‑1 (MQ=255)
                                                                                                              tatGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAACCTCCCCCTTAAAGTGAATTTCGCCAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGAt                                      <  1:345699/149‑1 (MQ=255)
                                                                                                              tatGGGTTGTCGTGACTTGTCCAAGATCTATGTTGTTATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAACCTCCCCCTTAAAGTGAATTTCGCCAGAACCTTCATCAGCAGCATAAACACTTCCAGTGAAGAGAAGAGAt                                      >  1:330529/1‑149 (MQ=255)
                                                                                                                  ggTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAACCTCCCCCTTAAAGTGAATTTCGCCAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACgg                                  >  1:709162/1‑149 (MQ=255)
                                                                                                                  ggTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAACCTCCCCCTTAAAGTGAATTTCGCCAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACgg                                  >  1:432423/1‑149 (MQ=255)
                                                                                                                  ggTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAACCTCCCCCTTAAAGTGAATTTCGCCAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACgg                                  >  2:1040553/1‑149 (MQ=255)
                                                                                                                    ttGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAACCTCCCCCTTAAAGTGAATTTCGCCAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGcc                                >  1:897635/1‑149 (MQ=255)
                                                                                                                               tGTCCAAGATCTATGTTTTTACCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAACCTCCCCCTTAAAGTGAATTTCGCCAGAACCTTCATCAGCAGCATAAATAGATCCAGTGAACAGCAGAGTTACGGCCAGTGCGGCCaa                     <  1:955145/149‑1 (MQ=255)
                                                                                                                               tGTCCAAGATCTATGTTTTTACCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAACCTCCCCCTTAAAGTGAATTTCGCCAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCaa                     <  1:257553/149‑1 (MQ=255)
                                                                                                                               tGTCCAAGATCTATGTTTTTACCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAACCTCCCCCTTAAAGTGAATTTCGCCAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCaa                     <  1:219341/149‑1 (MQ=255)
                                                                                                                               gttCCAAGCTCTAGGTTTTTTTCAATATCTTCTGGATGAATTTCTCACGGTGCTTCAATAACCTCCCCCTTAAAGGGAATTTCGCCAGGACCTTCATCAGCAGCATAACCAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCaa                     <  2:1228165/147‑1 (MQ=255)
                                                                                                                                       aTCTATGTTTTTATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAACCTCCCCCTTAAAGTGAATTTCGCCAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGtttttt             >  1:104461/1‑149 (MQ=255)
                                                                                                                                         cTATGTTTTTATCAATACCTTCTGGATGAATTTCACAAGGTGCTTCAATAACCTCCCCCTTAAAGTGAATTTCGCCAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGt           <  2:851306/149‑1 (MQ=255)
                                                                                                                                          tATGTTTTTATCAATATCTTCTGGATGAATTTCATAAGGTGCTTCAATAACCTCCCCCTTAAAGTGAATTTCGCCAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGTc          <  1:200735/149‑1 (MQ=255)
                                                                                                                                          tATGTTTTTATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAACCTCCCCCTTAAAGTGAATTTCGCCAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGTc          <  2:721424/149‑1 (MQ=255)
                                                                                                                                            tGTTTTTATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAACCTCCCCCTTAAAGTGAATTTCGCCAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGTCCt        <  1:152558/149‑1 (MQ=255)
                                                                                                                                               ttttATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAACCTCCCCCTTAAAGTGAATTTCGCCAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGTCCTTTa     <  2:709514/149‑1 (MQ=255)
                                                                                                                                                  tATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATACCCTCCCCCTTAAAGTGAATTTCGCCAGAACCTTCAGCAGCAGCATAAACAGGTGCAGTTTTCAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGTCCTTTAAAc  <  2:1181618/149‑1 (MQ=255)
                                                                                                                                                  tATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAACCTCCCCCTTAAAGTGAATTTCGCCAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGTCCTTTAAAc  <  2:422002/149‑1 (MQ=255)
                                                                                                                                                  tATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAACCTCCCCCTTAAAGTGAATTTCGCCAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGTCCTTTAAAc  <  2:1109420/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TACGCCAACTTTGGATACCGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGGTTGTCGTGACTTGTCCAAGATCTATGTTTTTATCAATATCTTCTGGATGAATTTCACAAGGTGCTTCAATAACCTCCCCCTTAAAGTGAATTTCGCCAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGTCCTTTAAAC  >  NC_000913/752499‑752793

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