Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,549,509 G→A A234A (GCG→GCA fimH → minor component of type 1 fimbriae

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,549,5090GA98.7% 275.3 / ‑4.1 75A234A (GCG→GCAfimHminor component of type 1 fimbriae
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (38/36);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01

CGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCGTCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAATAACACGGTATCGTTAGGAGCAGTAGGGACTTCGGCGGTGAGTCTGGGATTAACGGCAAATTATGCACGTACCGGAGGG  >  NC_000913/4549366‑4549653
                                                                                                                                               |                                                                                                                                                
cGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGtt                                                                                                                                             <  1:386656/149‑1 (MQ=255)
cGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGtt                                                                                                                                             <  1:674007/149‑1 (MQ=255)
     tGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCt                                                                                                                                      >  1:346743/1‑151 (MQ=255)
      gTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCt                                                                                                                                      <  1:806238/150‑1 (MQ=255)
         aCCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGcac                                                                                                                                  >  2:1042687/1‑151 (MQ=255)
           cGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAg                                                                                                                                <  2:144032/151‑1 (MQ=255)
             ttACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCACAGATGCGGGCAACTCGATTTTCAACAATACCGCATCGTTTTCACCTGCACAggg                                                                                                                              >  2:685971/1‑151 (MQ=255)
              tACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGc                                                                                                                             >  2:814818/1‑151 (MQ=255)
               aCTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGc                                                                                                                             <  2:1124275/150‑1 (MQ=255)
                ctctGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATTCCGCATCGTTTTCACCTGCACAggg                                                                                                                              <  1:304123/148‑1 (MQ=255)
                ctctGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGc                                                                                                                             <  1:81169/149‑1 (MQ=255)
                  ctGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGt                                                                                                                           <  1:899711/149‑1 (MQ=255)
                    gCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTc                                                                                                                          <  1:886815/148‑1 (MQ=255)
                    gCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCg                                                                                                                         <  1:757801/149‑1 (MQ=255)
                    gCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGc                                                                                                                       <  1:1220847/151‑1 (MQ=255)
                       ggACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGGGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGAGGCGGGCAACTCGATTGTCACCAATACCGCATCGTTTTCACCTGCACAGGGGGTCGGCGTa                                                                                                                    <  1:1224000/151‑1 (MQ=255)
                         aCTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTAc                                                                                                                   >  1:1014620/1‑150 (MQ=255)
                         aCTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACa                                                                                                                  <  2:84121/151‑1 (MQ=255)
                            aCCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAg                                                                                                                 <  2:259901/149‑1 (MQ=255)
                             cccTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCACATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTTCAGt                                                                                                                <  1:829065/149‑1 (MQ=255)
                               cTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGAc                                                                                                            >  1:366018/1‑151 (MQ=255)
                                 ggTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGAcgc                                                                                                          >  2:75610/1‑151 (MQ=255)
                                  gTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGAcgcg                                                                                                         >  1:644129/1‑151 (MQ=255)
                                  gTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGAcgcg                                                                                                         >  2:501902/1‑151 (MQ=255)
                                         gCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACgg                                                                                                   <  2:571819/150‑1 (MQ=255)
                                         gCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGt                                                                                                  >  2:691746/1‑151 (MQ=255)
                                          ccAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGTCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCTCCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACgg                                                                                                   <  1:713286/149‑1 (MQ=255)
                                            aaTTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACg                                                                                               <  1:954882/151‑1 (MQ=255)
                                             aTTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACg                                                                                               >  2:788320/1‑150 (MQ=255)
                                                ccTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGatt                                                                                            >  1:736562/1‑150 (MQ=255)
                                                 ctctTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGattat                                                                                          >  1:1106060/1‑151 (MQ=255)
                                                     tACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTcc                                                                                       <  1:208511/150‑1 (MQ=255)
                                                     tACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTcc                                                                                       >  2:960045/1‑150 (MQ=255)
                                                        cGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGTCGCGCAACGGTACGATTATTCCAGCg                                                                                   <  1:557505/151‑1 (MQ=255)
                                                        cGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCg                                                                                   <  1:860141/151‑1 (MQ=255)
                                                           ttATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAAt                                                                                >  1:68269/1‑151 (MQ=255)
                                                           ttATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAAt                                                                                >  2:470668/1‑151 (MQ=255)
                                                            tATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAAt                                                                                <  1:1080603/150‑1 (MQ=255)
                                                             aTTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAATa                                                                               <  2:797211/150‑1 (MQ=255)
                                                                gtgCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAATAacac                                                                           >  2:564071/1‑151 (MQ=255)
                                                                    gAAACGCCACGACCTGGGGCTTTACCTCTTCGCGACAACCGCAGATGCGGGCAACTCGATTATCACCCATACCGCATCGTTTTCACCTGCACAGGTCGTCCGCGTACAGTTGACGCGCACCGGAACGATTATTCGAGCGAATACCACGGt                                                                        >  1:1143127/1‑150 (MQ=255)
                                                                     aaaaGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAATAACACGGTAt                                                                      <  1:501887/151‑1 (MQ=255)
                                                                       aaGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCGTCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAATAACACGGTATCg                                                                    >  2:1047208/1‑151 (MQ=255)
                                                                          ccAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAATAACACGGTATCg                                                                    >  1:1098452/1‑148 (MQ=255)
                                                                           cAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAATAACACGGTATCGtt                                                                  >  2:350499/1‑149 (MQ=255)
                                                                           cAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAATAACACGGTATCGTTAg                                                                <  2:744868/151‑1 (MQ=255)
                                                                            aaaaCCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAATAACACGGTATCGTTAgg                                                               >  1:184759/1‑151 (MQ=255)
                                                                                    gggTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAATAACACGGTATCGTTAGGAGCAGTAg                                                       <  1:470850/151‑1 (MQ=255)
                                                                                     ggTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAATAACACGGTATCGTTAGGAGCAGTAgg                                                      >  2:1196402/1‑151 (MQ=255)
                                                                                        aTTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAATAACACGGTATCGTTAGGAGCAGTAGGGAc                                                   <  1:789330/151‑1 (MQ=255)
                                                                                           aCCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAATAACACGGTATCGTTAGGAGCAGTAGGGACTTc                                                >  1:1090658/1‑151 (MQ=255)
                                                                                           aCCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAATAACACGGTATCGTTAGGAGCAGTAGGGACTTc                                                >  2:597731/1‑151 (MQ=255)
                                                                                           aCCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAATAACACGGTATCGTTAGGAGCAGTAGGGACTTc                                                <  2:258704/151‑1 (MQ=255)
                                                                                               ctcCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAATAACACGGTATCGTTAGGAGCAGTAGGGACTTcggcg                                            >  1:130226/1‑151 (MQ=255)
                                                                                                tcCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAATAACACGGTATCGTTAGGAGCAGTAGGGACTTcggcgg                                           >  1:275800/1‑151 (MQ=255)
                                                                                                 ccGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAATAACACGGTATCGTTAGGAGCAGTAGGGACTTcggc                                             <  2:296076/148‑1 (MQ=255)
                                                                                                 ccGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAATAACACGGTATCGTTAGGAGCAGTAGGGACTTCGGCGGt                                          <  1:451740/151‑1 (MQ=255)
                                                                                                  cGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAATAACACGGTATCGTTAGGAGCAGTAGGGACTTCGGCGGTg                                         >  1:389393/1‑151 (MQ=255)
                                                                                                  cGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAATAACACGGTATCGTTAGGAGCAGTAGGGACTTCGGCGGTg                                         >  2:479145/1‑151 (MQ=255)
                                                                                                  cGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAATAACACGGTATCGTTAGGAGCAGTAGGGACTTCGGCGGTg                                         <  2:1241239/151‑1 (MQ=255)
                                                                                                      acaACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAATAACACGGTATCGTTAGGAGCAGTAGGGACTTCGGCGGTGAg                                       >  1:644694/1‑149 (MQ=255)
                                                                                                        aaCCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAATAACACGGTATCGTTAGGAGCAGTAGGGACTTCGGCGGGGAGTCTg                                   >  2:249293/1‑151 (MQ=255)
                                                                                                              aGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAATAACACGGTATCGTTAGGAGCAGTAGGGACTTCGGCGGTGAGTCTGGGATTa                             <  2:307281/151‑1 (MQ=255)
                                                                                                              aGATGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAATAACACGGTATCGTTAGGAGCAGTAGGGACTTCGGCGGTGAGTCTGGGATTa                             <  1:685946/151‑1 (MQ=255)
                                                                                                                aTGCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAATAACACGGTATCGTTAGGAGCAGTAGGGACTTCGGCGGTGAGTCTGGGATTAAc                           <  2:223562/151‑1 (MQ=255)
                                                                                                                  gCGGGCAACTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAATAACACGGTATCGTTAGGAGCAGTAGGGACTTCGGCGGTGAGTCTGGGATTAACg                          <  2:932262/150‑1 (MQ=255)
                                                                                                                         aCTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAATAACACGGTATCGTTAGGAGCAGTAGGGACTTCGGCGGTGAGTCTGGGATTAACGGCAAAtt                   >  2:347617/1‑150 (MQ=255)
                                                                                                                          cTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAATAACACGGTATCGTTAGGAGCAGTAGGGACTTCGGCGGTGAGTCTGGGATTAACGGCAAATTAt                 >  1:1100355/1‑151 (MQ=255)
                                                                                                                          cTCGATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAATAACACGGTATCGTTAGGAGCAGTAGGGACTTCGGCGGTGAGTCTGGGATTAACGGCAAATTAt                 >  2:337971/1‑151 (MQ=255)
                                                                                                                             gATTTTCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAATAACACGGTCTCGTTAGGAGCAGTAGGGACTTCGGCGGTGAGTCTGGGATTAACGGCAAATTATg                >  2:1043804/1‑149 (MQ=255)
                                                                                                                                 ttCACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAATAACACGGTATCGTTAGGAGCAGTAGGGACTTCGGCGGTGAGTCTGGGATTAACGGCAAATTATGCACGt           >  1:1186429/1‑150 (MQ=255)
                                                                                                                                   cACCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAATAACACGGTATCGTTAGGAGCAGTAGGGACTTCGGCGGTGAGTCTGGGATTAACGGCAAATTATGCACGTAcc        >  2:487870/1‑151 (MQ=255)
                                                                                                                                    aCCAATACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAATAACACGGTATCGTTAGGAGCAGTAGGGACTTCGGCGGTGAGTCTGGGATTAACGGCAAATTATGCACGTACCg       >  1:7641/1‑151 (MQ=255)
                                                                                                                                       aaTACCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAATAACACGGTATCGTTAGGAGCAGTAGGGACTTCGGCGGTGAGTCTGGGATTAACGGCAAATTATGCACGTACCgg      >  2:106349/1‑149 (MQ=255)
                                                                                                                                          aCCGCATCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAATAACACGGTATCGTTAGGAGCAGTAGGGACTTCGGCGGTGAGTCTGGGATTAACGGCAAATTATGCACGTACCGGAggg  <  1:283708/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
CGTGATGTCACCGTTACTCTGCCGGACTACCCTGGTTCAGTGCCAATTCCTCTTACCGTTTATTGTGCGAAAAGCCAAAACCTGGGGTATTACCTCTCCGGCACAACCGCAGATGCGGGCAACTCGATTTTCACCAATACCGCGTCGTTTTCACCTGCACAGGGCGTCGGCGTACAGTTGACGCGCAACGGTACGATTATTCCAGCGAATAACACGGTATCGTTAGGAGCAGTAGGGACTTCGGCGGTGAGTCTGGGATTAACGGCAAATTATGCACGTACCGGAGGG  >  NC_000913/4549366‑4549653

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