Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,247,031 G→A 7.4% intergenic (‑262/+32) 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,0310GA 7.4% 372.6 / 10.7 123intergenic (‑262/+32)cirA/lysPcolicin IA outer membrane receptor and translocator; ferric iron‑catecholate transporter/lysine transporter
Reads supporting (aligned to +/- strand):  ref base G (46/66);  new base A (6/3);  total (53/70)
Fisher's exact test for biased strand distribution p-value = 1.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 126 total aligned reads displayed.

CATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGTTGCC  >  NC_000913/2246882‑2247180
                                                                                                                                                      |                                                                                                                                                     
catcCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCgg                                                                                                                                                       <  1:685786/150‑1 (MQ=255)
 atcCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGt                                                                                                                                                      <  2:778515/150‑1 (MQ=255)
  tcCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGtt                                                                                                                                                     <  2:1347329/150‑1 (MQ=255)
  tcCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGtt                                                                                                                                                     >  2:579069/1‑150 (MQ=255)
   cccAGAAAGCCGGAACAGCCTCTTATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTg                                                                                                                                                    <  2:482186/150‑1 (MQ=255)
   cccAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTg                                                                                                                                                    <  2:434464/150‑1 (MQ=255)
   cccAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTg                                                                                                                                                    <  1:1651018/150‑1 (MQ=255)
    ccAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTga                                                                                                                                                   >  2:714573/1‑150 (MQ=255)
      aGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTgaga                                                                                                                                                 <  1:579268/150‑1 (MQ=255)
         aaGCCGGAACTGCCGCTGATTTGTCAGCCTTTTTAGATATCTGCGTAACGATAATGCAAACGATAACAATTATCAAACCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGCAATGGAATTG‑AAAACGCCCTCTCGGTTGAGAggg                                                                                                                                              <  2:256982/150‑1 (MQ=255)
          aGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGc                                                                                                                                             <  1:553022/150‑1 (MQ=255)
           gCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCt                                                                                                                                            >  1:1531421/1‑150 (MQ=255)
           gCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCAACCGAGAGGGCt                                                                                                                                            <  2:1438972/150‑1 (MQ=255)
            ccGGAACAGCCTCTGATTTGTCATCCTTTTTTGATTGCGGCGTAAAGATAATTCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCAACCGAGAGGGCtt                                                                                                                                           <  2:161315/150‑1 (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)
                          gATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGggag                                                                                                                             <  1:1018135/150‑1 (MQ=255)
                           aTTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGggagg                                                                                                                            <  1:621446/150‑1 (MQ=255)
                             ttGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAATCGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGaa                                                                                                                          <  2:1344758/150‑1 (MQ=255)
                             ttGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGaa                                                                                                                          <  1:1621501/150‑1 (MQ=255)
                               gTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGaaaa                                                                                                                        >  2:1210940/1‑149 (MQ=255)
                               gTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAAc                                                                                                                        >  2:522232/1‑150 (MQ=255)
                                 cAGCCTTTTTTGATAGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACtt                                                                                                                      >  2:1340150/1‑150 (MQ=255)
                                    cattttttGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTAtt                                                                                                                   <  1:326513/148‑1 (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)
                                            gATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCg                                                                                                           <  1:1613900/150‑1 (MQ=255)
                                                     gTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCggg                                                                                                  <  2:1259453/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)
                                                          gATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGATGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACtt                                                                                             >  2:564664/1‑150 (MQ=255)
                                                           aTAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTc                                                                                            >  2:1520007/1‑150 (MQ=255)
                                                            tAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCa                                                                                           <  2:408577/150‑1 (MQ=255)
                                                                gCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTc                                                                                       <  1:573550/150‑1 (MQ=255)
                                                                gCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTc                                                                                       >  1:602289/1‑150 (MQ=255)
                                                                 cAAGCGATAACAATTATCAATCCAAATGTT‑AATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGc                                                                                     >  2:94103/1‑150 (MQ=255)
                                                                 cAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCg                                                                                      <  1:786218/150‑1 (MQ=255)
                                                                  aaaCGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGc                                                                                     >  1:925709/1‑150 (MQ=255)
                                                                    aCGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTg                                                                                   <  1:508191/150‑1 (MQ=255)
                                                                     cGATAACTATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAATTGTAATTG‑AAAAAGCGCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGt                                                                                  <  2:255595/150‑1 (MQ=255)
                                                                       aTAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAg                                                                                <  2:1109950/150‑1 (MQ=255)
                                                                         aaCAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCg                                                                              >  2:1034000/1‑150 (MQ=255)
                                                                         aaCAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCATATCGGTTGAGAGGGCTTACCAAGGAAGGGAGGAAACTTTTTTCTTATCGTTCTCCGGGAACTTAAGTTCGCTGTAGCg                                                                              >  2:850664/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)
                                                                               tATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGaa                                                                        >  1:1209403/1‑150 (MQ=255)
                                                                                    aTCCAAATTTTAAAATTTATATCCGGATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAg                                                                   <  1:41237/150‑1 (MQ=255)
                                                                                    aTCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGat                                                                   >  1:986610/1‑149 (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)
                                                                                           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)
                                                                                                   tttATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCt                                                    <  1:924832/150‑1 (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                                  >  1:620112/1‑150 (MQ=255)
                                                                                                                     gTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAaataa                                  >  2:1707011/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)
                                                                                                                             aaTGGAATTT‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAAc                          <  1:222731/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)
                                                                                                                                            aGCCCTCTCAACCGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCGGGATACCAat            >  1:1249560/1‑150 (MQ=255)
                                                                                                                                            aGCCCTCTCAACCGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCGGGATACCAat            >  1:627464/1‑150 (MQ=255)
                                                                                                                                                ctctCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATAc        <  2:602289/150‑1 (MQ=255)
                                                                                                                                                 tctcGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACg       <  1:105226/150‑1 (MQ=255)
                                                                                                                                                  ctcGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGt      >  2:1175634/1‑150 (MQ=255)
                                                                                                                                                     ggTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGTTGc   <  2:622141/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTG‑AAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGTTGCC  >  NC_000913/2246882‑2247180

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