Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,247,033 T→C 7.4% intergenic (‑264/+30) 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,0330TC 7.4% 350.6 / 14.2 122intergenic (‑264/+30)cirA/lysPcolicin IA outer membrane receptor and translocator; ferric iron‑catecholate transporter/lysine transporter
Reads supporting (aligned to +/- strand):  ref base T (47/66);  new base C (6/3);  total (53/69)
Fisher's exact test for biased strand distribution p-value = 1.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.48e-01

Only 100 of 124 total aligned reads displayed.

TCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGTTGCC  >  NC_000913/2246884‑2247180
                                                                                                                                                      |                                                                                                                                                   
tcCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGtt                                                                                                                                                     <  2:1347329/150‑1 (MQ=255)
 cccAGAAAGCCGGAACAGCCTCTTATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTg                                                                                                                                                    <  2:482186/150‑1 (MQ=255)
 cccAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTg                                                                                                                                                    <  1:1651018/150‑1 (MQ=255)
 cccAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTg                                                                                                                                                    <  2:434464/150‑1 (MQ=255)
 cccAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTg                                                                                                                                                    <  2:1187984/150‑1 (MQ=255)
  ccAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTga                                                                                                                                                   >  2:714573/1‑150 (MQ=255)
    aGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTgaga                                                                                                                                                 <  1:579268/150‑1 (MQ=255)
     taaaGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGagag                                                                                                                                                <  1:1710298/149‑1 (MQ=255)
     gAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGagag                                                                                                                                                <  1:1419360/150‑1 (MQ=255)
       aaGCCGGAACTGCCGCTGATTTGTCAGCCTTTTTAGATATCTGCGTAACGATAATGCAAACGATAACAATTATCAAACCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGCAATGGAATTGAAAACGCCCTCTCGGTTGAGAggg                                                                                                                                              <  2:256982/150‑1 (MQ=255)
        aGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGc                                                                                                                                             <  1:553022/150‑1 (MQ=255)
         gCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCAACCGAGAGGGCt                                                                                                                                            <  2:1438972/150‑1 (MQ=255)
         gCCGGAACAGCCTATGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCt                                                                                                                                            <  2:604549/150‑1 (MQ=255)
          ccGGAACAGCCTCTGATTTGTCATCCTTTTTTGATTGCGGCGTAAAGATAATTCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCAACCGAGAGGGCtt                                                                                                                                           <  2:161315/150‑1 (MQ=255)
             gAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGc                                                                                                                                        <  2:834151/150‑1 (MQ=255)
              aaCAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCtttttc                                                                                                                                       <  1:1699937/150‑5 (MQ=255)
              aaCAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCAACCGAGAGGGCtttttc                                                                                                                                       <  1:692712/150‑5 (MQ=255)
               aCAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCa                                                                                                                                       <  1:1672410/150‑1 (MQ=255)
                cAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCaag                                                                                                                                     <  2:144962/150‑1 (MQ=255)
                cAGCCTCTGATTTGTAAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCaag                                                                                                                                     <  1:111176/150‑1 (MQ=255)
                  gCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCaagga                                                                                                                                   >  2:621446/1‑150 (MQ=255)
                  gCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCaagga                                                                                                                                   >  2:800267/1‑150 (MQ=255)
                   ccTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCaaggaa                                                                                                                                  >  2:174065/1‑150 (MQ=255)
                      ctGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAggg                                                                                                                               >  2:774467/1‑150 (MQ=255)
                       tGATTTGTCAGCCTTTTTTGATTGCGGAGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGgga                                                                                                                              <  2:412385/150‑1 (MQ=255)
                        gATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGggag                                                                                                                             <  1:1018135/150‑1 (MQ=255)
                         aTTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGggagg                                                                                                                            <  1:621446/150‑1 (MQ=255)
                           ttGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAATCGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGaa                                                                                                                          <  2:1344758/150‑1 (MQ=255)
                             gTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAAc                                                                                                                        >  2:522232/1‑150 (MQ=255)
                               cAGCCTTTTTTGATAGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACtt                                                                                                                      >  2:1340150/1‑150 (MQ=255)
                                  cattttttGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTAtt                                                                                                                   <  1:326513/148‑1 (MQ=255)
                                   cTTTTTTGATTGCGGCGTAACGAGAATGCAACCGATACCAATTCTCAAGCCAAAGGTTAAATTTTATATCCGCATGATTTGTCAGGGAAAGGGAATTGAAAAAGCCCTCTAGTTTGAGAGGGCGTAGCAAGGAAGGGAGTACACTTAatt                                                                                                                  >  2:858573/1‑150 (MQ=255)
                                    ttttttGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTc                                                                                                                 <  2:1085802/150‑1 (MQ=255)
                                    ttttttGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGACAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTc                                                                                                                 >  1:1392565/1‑150 (MQ=255)
                                     tttttGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCt                                                                                                                <  2:451623/150‑1 (MQ=255)
                                       tttGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTa                                                                                                              <  1:30818/150‑1 (MQ=255)
                                         tGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATc                                                                                                            >  2:498425/1‑150 (MQ=255)
                                         tGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGCGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTAta                                                                                                            >  1:1264695/1‑149 (MQ=255)
                                          gATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCg                                                                                                           <  1:1613900/150‑1 (MQ=255)
                                                ggagTAACGATAATGCAAACGATAACAATTATCAAACCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAATCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGc                                                                                                     <  1:345561/147‑1 (MQ=255)
                                                   gTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCggg                                                                                                  <  2:1259453/150‑1 (MQ=255)
                                                    tAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGa                                                                                                 >  2:960468/1‑150 (MQ=255)
                                                    tAACGAAAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAAAGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGa                                                                                                 <  2:1100926/150‑1 (MQ=255)
                                                       cGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACt                                                                                              <  1:611626/150‑1 (MQ=255)
                                                        gATAATGCAAACGATAACAATTATCAATCCAAATGTTAGATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACtt                                                                                             <  1:504814/150‑1 (MQ=255)
                                                        gATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGATGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACtt                                                                                             >  2:564664/1‑150 (MQ=255)
                                                          tAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCa                                                                                           <  2:408577/150‑1 (MQ=255)
                                                              gCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTc                                                                                       >  1:602289/1‑150 (MQ=255)
                                                               cAAGCGATAACAATTATCAATCCAAATGTTAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGc                                                                                     >  2:94103/1‑150 (MQ=255)
                                                               cAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCg                                                                                      <  1:786218/150‑1 (MQ=255)
                                                                aaaCGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGc                                                                                     >  1:622141/1‑150 (MQ=255)
                                                                aaaCGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGc                                                                                     >  1:925709/1‑150 (MQ=255)
                                                                  aCGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTg                                                                                   <  1:508191/150‑1 (MQ=255)
                                                                  aCGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTg                                                                                   >  1:579711/1‑150 (MQ=255)
                                                                   cGATAACTATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAATTGTAATTGAAAAAGCGCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGt                                                                                  <  2:255595/150‑1 (MQ=255)
                                                                     taaaCAATTATAAATACAAAAGGTAAAATTTACTTCCGCATTCTTTTGCACGGAAATGGAAATTAAAAAGCCCTCACGGTTGTGCTGGGCTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAg                                                                                <  1:733789/148‑1 (MQ=255)
                                                                     aTAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAg                                                                                <  2:1109950/150‑1 (MQ=255)
                                                                       aaCAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCg                                                                              >  2:1034000/1‑150 (MQ=255)
                                                                         cAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTa                                                                            >  2:905555/1‑150 (MQ=255)
                                                                          aaTTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTAc                                                                           <  2:1425199/150‑1 (MQ=255)
                                                                           aTTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACg                                                                          <  2:1484436/150‑1 (MQ=255)
                                                                             tATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGaa                                                                        >  1:1209403/1‑150 (MQ=255)
                                                                                  aTCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGat                                                                   >  1:986610/1‑149 (MQ=255)
                                                                                     cAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGGTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTc                                                                <  1:1554794/150‑1 (MQ=255)
                                                                                        aTGTTAAATTTTATATCCGCATGATTTGTCTGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCtt                                                             <  1:248667/150‑1 (MQ=255)
                                                                                        aTGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCtt                                                             <  2:1349230/150‑1 (MQ=255)
                                                                                        aTGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCtt                                                             <  2:1218553/150‑1 (MQ=255)
                                                                                         tGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCttt                                                            <  1:774467/150‑1 (MQ=255)
                                                                                           ttAAATTTTATATCCGCATTATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGGTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGTGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTaa                                                          <  1:1222329/150‑1 (MQ=255)
                                                                                               aTTTTATATCCGCATGATTTGACAGGGAAATTGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGTAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGAACGAAGTGAGTTCCTTTAATCac                                                      >  2:389226/1‑149 (MQ=255)
                                                                                                 tttATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCt                                                    <  1:1415858/150‑1 (MQ=255)
                                                                                                 tttATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCt                                                    <  1:924832/150‑1 (MQ=255)
                                                                                                  ttATATCCTAATGATTTTTCAGTGAATTTTAATTGAAAAAGCACTATCGGGTTAGAGGGCTTAGCAAGGAAGGTAGGAAACTTATTTCTTATCGTTATGCTGGAACTTCATTTAGCTGTAGCGTACTAAGTGAGTTCCTTTAATCAGCtt                                                   <  1:1559511/150‑1 (MQ=255)
                                                                                                    atatCCGGAGGATTTGGCACGGAAATGGAATTGAAACAGCCCTCTGGGTTGAGAGGGCTTATCTATGAAGGGTGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGt                                                 <  1:1068496/150‑1 (MQ=255)
                                                                                                      atCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAg                                               >  1:798859/1‑150 (MQ=255)
                                                                                                       tCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGc                                              <  1:473079/150‑1 (MQ=255)
                                                                                                            aTGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAAc                                         >  2:514281/1‑150 (MQ=255)
                                                                                                             tGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAAcc                                        >  2:40378/1‑150 (MQ=255)
                                                                                                                   gTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAaataa                                  >  1:1023723/1‑150 (MQ=255)
                                                                                                                   gTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAaataa                                  >  2:1707011/1‑150 (MQ=255)
                                                                                                                       gggAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAg                              >  1:1616809/1‑150 (MQ=255)
                                                                                                                        ggAAATGGATTTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAgg                             <  2:554033/150‑1 (MQ=255)
                                                                                                                         taaaTGTAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGCGAACATAATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGa                            <  1:93798/149‑1 (MQ=255)
                                                                                                                           aaTGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAAc                          <  2:525995/150‑1 (MQ=255)
                                                                                                                              ggAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGc                       <  1:458997/150‑1 (MQ=255)
                                                                                                                               gAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCg                      >  2:51463/1‑150 (MQ=255)
                                                                                                                                 aTTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGATCCTTTAATCAGCTTGTAGCCGAACCAAATAATCATGAACTGCggg                    >  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)
                                                                                                                                          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)
                                                                                                                                                    gTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGTTGcc  >  1:87479/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
TCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGTTGCC  >  NC_000913/2246884‑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: