Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,247,032 T→C 7.5% intergenic (‑263/+31) cirA ← / ← lysP colicin IA outer membrane receptor and translocator; ferric iron‑catecholate transporter/lysine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,247,0320TC 7.5% 336.7 / 13.4 122intergenic (‑263/+31)cirA/lysPcolicin IA outer membrane receptor and translocator; ferric iron‑catecholate transporter/lysine transporter
Reads supporting (aligned to +/- strand):  ref base T (47/62);  new base C (6/3);  total (53/69)
Fisher's exact test for biased strand distribution p-value = 2.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.20e-01

Only 100 of 125 total aligned reads displayed.

ATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGTTGCC  >  NC_000913/2246883‑2247180
                                                                                                                                                      |                                                                                                                                                    
atcCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGt                                                                                                                                                      <  2:778515/150‑1 (MQ=255)
 tcCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGtt                                                                                                                                                     >  2:579069/1‑150 (MQ=255)
 tcCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGtt                                                                                                                                                     <  2:1347329/150‑1 (MQ=255)
  cccAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTg                                                                                                                                                    <  2:1187984/150‑1 (MQ=255)
  cccAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTg                                                                                                                                                    <  2:434464/150‑1 (MQ=255)
  cccAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTg                                                                                                                                                    <  1:1651018/150‑1 (MQ=255)
     aGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTgaga                                                                                                                                                 <  1:579268/150‑1 (MQ=255)
         aGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGc                                                                                                                                             <  1:553022/150‑1 (MQ=255)
          gCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCt                                                                                                                                            >  1:1531421/1‑150 (MQ=255)
          gCCGGAACAGCCTATGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCt                                                                                                                                            <  2:604549/150‑1 (MQ=255)
           ccGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCtt                                                                                                                                           >  2:105226/1‑150 (MQ=255)
              gAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGc                                                                                                                                        <  2:834151/150‑1 (MQ=255)
               aaCAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCtttttc                                                                                                                                       <  1:1699937/150‑5 (MQ=255)
               aaCAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCAACCGAGAGGGCtttttc                                                                                                                                       <  1:692712/150‑5 (MQ=255)
                aCAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCa                                                                                                                                       <  1:1672410/150‑1 (MQ=255)
                aCAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGATTAGCaa                                                                                                                                      >  2:452506/1‑150 (MQ=255)
                 cAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACTATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCaag                                                                                                                                     <  2:398196/150‑1 (MQ=255)
                 cAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCaag                                                                                                                                     <  2:144962/150‑1 (MQ=255)
                   gCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCaagga                                                                                                                                   >  2:800267/1‑150 (MQ=255)
                   gCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCaagga                                                                                                                                   >  2:621446/1‑150 (MQ=255)
                    ccTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCaaggaa                                                                                                                                  >  2:174065/1‑150 (MQ=255)
                       ctGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAggg                                                                                                                               >  2:774467/1‑150 (MQ=255)
                        tGATTTGTCAGCCTTTTTTGATTGCGGAGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGgga                                                                                                                              <  2:412385/150‑1 (MQ=255)
                         gATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGggag                                                                                                                             <  1:1018135/150‑1 (MQ=255)
                          aTTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGggagg                                                                                                                            <  1:621446/150‑1 (MQ=255)
                            ttGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGaa                                                                                                                          <  1:1621501/150‑1 (MQ=255)
                              gTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGaaaa                                                                                                                        >  2:1210940/1‑149 (MQ=255)
                                cAGCCTTTTTTGATAGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACtt                                                                                                                      >  2:1340150/1‑150 (MQ=255)
                                   cattttttGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTAtt                                                                                                                   <  1:326513/148‑1 (MQ=255)
                                    cTTTTTTGATTGCGGCGTAACGAGAATGCAACCGATACCAATTCTCAAGCCAAAGGTTAAATTTTATATCCGCATGATTTGTCAGGGAAAGGGAATTG‑AAAAAGCCCTCTAGTTTGAGAGGGCGTAGCAAGGAAGGGAGTACACTTAatt                                                                                                                  >  2:858573/1‑150 (MQ=255)
                                     ttttttGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTc                                                                                                                 <  2:1085802/150‑1 (MQ=255)
                                     ttttttGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGACAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTc                                                                                                                 >  1:1392565/1‑150 (MQ=255)
                                      tttttGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCt                                                                                                                <  2:451623/150‑1 (MQ=255)
                                        tttGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTa                                                                                                              <  1:30818/150‑1 (MQ=255)
                                          tGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATc                                                                                                            >  2:498425/1‑150 (MQ=255)
                                           gATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCg                                                                                                           <  1:1613900/150‑1 (MQ=255)
                                                     tAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGa                                                                                                 >  2:960468/1‑150 (MQ=255)
                                                     tAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATAG‑AAAAAGCCATCACGGTTCAGAGGGCTAAGCAAGGAATGGAGCAAAATTATTTCTTAACGTTCTGCGGGa                                                                                                 >  2:517399/1‑150 (MQ=255)
                                                     tAACGAAAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAAAGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGa                                                                                                 <  2:1100926/150‑1 (MQ=255)
                                                       aCGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAAc                                                                                               <  1:694040/150‑1 (MQ=255)
                                                        cGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACt                                                                                              <  1:611626/150‑1 (MQ=255)
                                                         gATAATGCAAACGATAACAATTATCAATCCAAATGTTAGATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACtt                                                                                             <  1:504814/150‑1 (MQ=255)
                                                         gATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGATGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACtt                                                                                             >  2:564664/1‑150 (MQ=255)
                                                          aTAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTc                                                                                            >  2:1520007/1‑150 (MQ=255)
                                                               gCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTc                                                                                       <  1:573550/150‑1 (MQ=255)
                                                                cAAGCGATAACAATTATCAATCCAAATGTT‑AATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGc                                                                                     >  2:94103/1‑150 (MQ=255)
                                                                 aaaCGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGc                                                                                     >  1:925709/1‑150 (MQ=255)
                                                                 aaaCGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGc                                                                                     >  1:622141/1‑150 (MQ=255)
                                                                   aCGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTg                                                                                   >  1:579711/1‑150 (MQ=255)
                                                                    cGATAACTATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAATTGTAATTG‑AAAAAGCGCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGt                                                                                  <  2:255595/150‑1 (MQ=255)
                                                                      taaaCAATTATAAATACAAAAGGTAAAATTTACTTCCGCATTCTTTTGCACGGAAATGGAAATT‑AAAAAGCCCTCACGGTTGTGCTGGGCTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAg                                                                                <  1:733789/148‑1 (MQ=255)
                                                                        aaCAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCg                                                                              >  2:1034000/1‑150 (MQ=255)
                                                                         aCAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGATAGGGCTTAGCAAGGAAGGGAGGAAACTTATATCTTATCTTTCTGCGGGAACTTCATTTCGCTGTAGctt                                                                             >  2:1300676/1‑148 (MQ=255)
                                                                          cAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTa                                                                            >  2:905555/1‑150 (MQ=255)
                                                                           aaTTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTAc                                                                           <  2:1425199/150‑1 (MQ=255)
                                                                            aTTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACg                                                                          <  2:1484436/150‑1 (MQ=255)
                                                                              tATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGaa                                                                        >  1:1209403/1‑150 (MQ=255)
                                                                                   aTCCAAATTTTAAAATTTATATCCGGATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAg                                                                   <  1:41237/150‑1 (MQ=255)
                                                                                      cAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGGTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTc                                                                <  1:1554794/150‑1 (MQ=255)
                                                                                         aTGTTAAATTTTATATCCGCATGATTTGTCTGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCtt                                                             <  1:248667/150‑1 (MQ=255)
                                                                                         aTGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCtt                                                             <  2:1218553/150‑1 (MQ=255)
                                                                                         aTGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCtt                                                             <  2:1349230/150‑1 (MQ=255)
                                                                                          tGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCttt                                                            <  1:774467/150‑1 (MQ=255)
                                                                                           gTTAAATTTTATATCCGCATGATTTGACAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTa                                                           <  2:1392565/150‑1 (MQ=255)
                                                                                            ttAAATTTTATATCCTCATGATTTTTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTaa                                                          <  1:930097/150‑1 (MQ=255)
                                                                                            ttAAATTTTATATCCGCATTATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGGTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGTGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTaa                                                          <  1:1222329/150‑1 (MQ=255)
                                                                                                aTTTTATATCCGCATGATTTGACAGGGAAATTGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGTAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGAACGAAGTGAGTTCCTTTAATCac                                                      >  2:389226/1‑149 (MQ=255)
                                                                                                  tttATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCt                                                    <  1:1415858/150‑1 (MQ=255)
                                                                                                   ttATATCCTAATGATTTTTCAGTGAATTTTAATTG‑AAAAAGCACTATCGGGTTAGAGGGCTTAGCAAGGAAGGTAGGAAACTTATTTCTTATCGTTATGCTGGAACTTCATTTAGCTGTAGCGTACTAAGTGAGTTCCTTTAATCAGCtt                                                   <  1:1559511/150‑1 (MQ=255)
                                                                                                     atatCCGGAGGATTTGGCACGGAAATGGAATTG‑AAACAGCCCTCTGGGTTGAGAGGGCTTATCTATGAAGGGTGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGt                                                 <  1:1068496/150‑1 (MQ=255)
                                                                                                       atCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAg                                               >  1:798859/1‑150 (MQ=255)
                                                                                                        tCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGc                                              <  1:473079/150‑1 (MQ=255)
                                                                                                          cGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCg                                            >  1:800726/1‑150 (MQ=255)
                                                                                                             aTGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAAc                                         >  2:514281/1‑150 (MQ=255)
                                                                                                              tGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAAcc                                        >  2:40378/1‑150 (MQ=255)
                                                                                                                    gTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAaataa                                  >  1:1023723/1‑150 (MQ=255)
                                                                                                                    gTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAaataa                                  >  2:1707011/1‑150 (MQ=255)
                                                                                                                    gTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAaataa                                  >  1:620112/1‑150 (MQ=255)
                                                                                                                      cAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATATCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATc                                >  2:651780/1‑150 (MQ=255)
                                                                                                                         ggAAATGGATTTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAgg                             <  2:554033/150‑1 (MQ=255)
                                                                                                                         ggAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAgg                             >  1:1138925/1‑150 (MQ=255)
                                                                                                                          taaaTGTAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGCGAACATAATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGa                            <  1:93798/149‑1 (MQ=255)
                                                                                                                            aaTGGAATTT‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAAc                          <  1:222731/150‑1 (MQ=255)
                                                                                                                            aaTGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAAc                          <  2:525995/150‑1 (MQ=255)
                                                                                                                            aaTGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTGCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAAc                          <  1:1324685/150‑1 (MQ=255)
                                                                                                                               ggAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGc                       <  1:458997/150‑1 (MQ=255)
                                                                                                                                gAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCg                      >  2:51463/1‑150 (MQ=255)
                                                                                                                                 aaTTG‑AACAAGCCCTCTCGCGTGAGAGGGCGTAGCAAGGAAGGGAGGAAACTTATTTCATATCGTTCTGCGGGAACTTCATTTCGCTGGAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCgg                     <  2:255392/150‑1 (MQ=255)
                                                                                                                                  aTTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGATCCTTTAATCAGCTTGTAGCCGAACCAAATAATCATGAACTGCggg                    >  2:302031/1‑150 (MQ=255)
                                                                                                                                      tgctaaGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCGGGATa                 >  1:1342125/5‑150 (MQ=255)
                                                                                                                                       gctaaGCCCTCTCAACCGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCGGGATAc                >  1:914878/4‑150 (MQ=255)
                                                                                                                                       gctaaGCCCTCTCAACCGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGAAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCGGGATAc                >  2:1283582/4‑150 (MQ=255)
                                                                                                                                        ctaaGCCCTCTCAACCGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGCACAGCGGGATAcc               >  2:1587352/3‑150 (MQ=255)
                                                                                                                                        ctaaGCCCTCTCAACCGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCGGGATAcc               >  1:96527/3‑150 (MQ=255)
                                                                                                                                        aaaaGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCGGGATAcc               <  2:1300551/150‑1 (MQ=255)
                                                                                                                                           aGCCCTCTCAACCGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCGGGATACCAat            >  1:1249560/1‑150 (MQ=255)
                                                                                                                                                tctcGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACg       <  1:105226/150‑1 (MQ=255)
                                                                                                                                                 ctcGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGt      >  2:1175634/1‑150 (MQ=255)
                                                                                                                                                    ggTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGTTGc   <  2:622141/150‑1 (MQ=255)
                                                                                                                                                     gTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGTTGcc  >  1:87479/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
ATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGTTGCC  >  NC_000913/2246883‑2247180

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