Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,245,130 A→T I547N (ATC→AAC)  cirA ← colicin IA outer membrane receptor and translocator; ferric iron‑catecholate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,245,1300AT98.0% 175.8 / ‑6.0 51I547N (ATC→AAC) cirAcolicin IA outer membrane receptor and translocator; ferric iron‑catecholate transporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (22/28);  minor base G (1/0);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

TCGTCACGACTGAGATCCTTGTCGCCAAGGTTCAGCACGCCTGCGCGCAGTTTGACGTCTTTAGTCACCTGCCAGGCCGCGCCGGTATTCCAGATGGTATAACCGCCCGGTGTTTTAGCCGTCGCGCTGTCGGCGCGTTTCTGCCCGGTATAGTGCCCAGAAACATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGTTACGGAAAACGGTAACGCTGGATTCAACCCCTTCCAGCCAGCCTTCTTCACCCATGTAGTAAAGCCCCAGCTCCCAACTTTCGCTGGTTTCT  >  NC_000913/2244836‑2245423
                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                                     
tcgtcACGACTGAGATCCTTGTCGCCAAGGTTCAGCACGCCTGCGCGCAGTTTGACGTCTTTAGTCACCTGCCAGGCCGCGCCGGTATTCCAGATGGTATAACCGCCCGGTGTTTTAGCCGTCGCGCTGTCAGCGCGTTTCTGCCCGGTATAGTGCCCAGAAACATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGa                                                                                                                                                                                                                                                                                                   <  1:447556/299‑1 (MQ=255)
acgtgtgctcttccgatcTTGTCGCCAAGGTTCAGCACGCCTGCGCGCAGTTTGACGTCTTTAGTCACCTGCCAGGCCGCGCCGGTATTCCAGATGGTATAACCGCCCGGTGTTTTAGCCGTCGCGCTGTCGGCGCGTTTCTGCCCGGTATAGTGCCCAGAAACATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGACa                                                                                                                                                                                                                                                                                                 <  2:111439/284‑1 (MQ=255)
           gagaTCCTTGTCGCCAAGGTTCAGCACGCCTGCGCGCAGTTTGACGTCTTTAGTCACCTGCCAGGCCGCGCCGGTATTCCAGATGGTATAACCGCCCGGTGTTTTAGCCGTCGCGCTGTCGGCGCGTTTCTGCCCGGTATAGTGCCCAGAAACATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCTTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACCTCACGACCATCGTTGTCGGTGTCGTTTTTCGCCAGTTTCCAtt                                                                                                                                                                                                                                                                                        >  2:74686/1‑299 (MQ=255)
                cctggtgccgacggTTCAGCACTCCTGCGCGCAGTTTGAGTTCTATAGGCACCTGCCAGGCCGCGCCGGTATTCCAGATGGTATAACCGCCCGGTGTTTTAGCCGTCGCGCTGTCGGCGCGTTTATGCCCGGTATAGTGCCCAGAAACATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGtt                                                                                                                                                                                                                                                                                 <  1:220201/289‑1 (MQ=255)
                 cTTGTCGCCAAGGTTCAGCACGCCTGCGCGCAGTTTGACGTCTTTAGTCACCTGCCAGGCCGCGCCGGTATTCCAGATGGTATAACCGCCCGGTGTTTTAGCCGTCGCGCTGTCGGCGCGTTTCTGCCCGGTATAGTGCCCAGAAACATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTg                                                                                                                                                                                                                                                                                <  2:278703/301‑1 (MQ=255)
                   tGTCGCCAAGGTTCAGCACGCCTGCGCGCAGTTTGACGTCTTTAGTCACCTGCCAGGCCGCGCCGGTATTCCAGATGGTATAACCGCCCGGTGTTTTAGCCGTCGCGCTGTCGGCGCGTTTCTGCCCGGTATAGTGCCCAGAAACATAGACTGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTGGCAGTATGGAACGGCGCATCGTATAGCGGATTGTTTTTGCCGTTGCTGACATCACGACCCGCGTTGTCGGACTAGTTATTCGAACGTTTCGTTTATTCGTTGaa                                                                                                                                                                                                                                                                              >  1:202968/1‑301 (MQ=255)
                   tGTCGCCAAGGTTCAGCACGCCTGCGCGCAGTTTGACGTCTTTAGTCACCTGCCAGGCCGCGCCGGTATTCCAGATGGTATAACCGCCCGGTGTTTTAGCCGTCGCGCTGTCGGCGCGTTTCTGCCCGGTATAGTGCCCAGAAACATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGACaagatcggatgagcgtcgtg                                                                                                                                                                                                                                                                              >  1:111299/1‑282 (MQ=255)
                       gCCAAGGTTCAGCACGCCTGCGCGCAGTTTGACGTCTTTAGTCACCTGCCAGGCCGCGCCGGTATTCCAGATGGTATAACCGCCCGGTGTTTTAGCCGTCGCGCTGTCGGCGCGTTTCTGCCCGGTATAGTGCCCAGAAACATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTTAACGGa                                                                                                                                                                                                                                                                          >  1:190162/1‑301 (MQ=255)
                        ccAAGGTTCAGCACGCCTGCGCGCAGTTTGACGTCTTTAGTCACCTGCCAGGCCGCGCCGGTATTCCAGATGGTATAACCGCCCGGTGTTTTAGCCGTCGCGCTGTCGGCGCGTTTCTGCACGGTATAGTGCCCAGAAACATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCAGCGTCGTAGGTGTGGTTGGTCGACAGTTTCCATTAATCGTTGAACggca                                                                                                                                                                                                                                                                         >  2:360688/1‑299 (MQ=255)
                             gTTCAGCACGCCTGCGCGCAGTTTGACGTCTTTAGTCACCTGCCAGGCCGCGCCGGTATTCCAGATGGTATAACCGCCCGGTGTTTTAGCCGTCGCGCTGTCGGCGCGTTTCTGCCCGGTATAGTGCCCAGAAACATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTGAACGGAAttt                                                                                                                                                                                                                                                                      <  2:401912/299‑1 (MQ=255)
                              ttCAGCACGCCTGCGCGCAGTTTGACGTCTTTAGTCACCTGCCAGGCCGCGCCGGTATTCCAGATGGTATAACCGCCCGGTGTTTTAGCCGTCGCGCTGTCGGCGCGTTTCTGCCCGGTATAGTGCCCAGAAACATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTGAACGGAATTTTc                                                                                                                                                                                                                                                                    <  2:190389/300‑1 (MQ=255)
                              tgcAGCGCGCCTGCGCGCAGTTTGACTTCGTTCGTCACCTGCCAGGCCGCGCCGGTATTCCAGATGGTATAACCGCCCGGTGTTTTAGCCGTCGCGCTGTCGGCGCGTTTCTGCCCGGTATAGTGCCCAGAAACATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTGAACGGAATTTTc                                                                                                                                                                                                                                                                    <  1:338622/298‑1 (MQ=255)
                                       ccTGCGCGCAGTTTGACGTCTTTAGTCACCTGCCAGGCCGCGCCGGTATTCCAGATGGTATAACCGCCCGGTGTTTTAGCCGTCGCGCTGTCGGCGCGTTTCTGCCCGGTATAGTGCCCAGAAACATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGttt                                                                                                                                                                                                                                                          >  2:48541/1‑301 (MQ=255)
                                               cAGTTTGACGTCTTTAGTCACCTGCCAGGCCGCGCCGGTATTCCAGATGGTATAACCGCCCGGTGTTTTAGCCGTCGCGCTGTCGGCGCGTTTCTGCCCGGTATAGTGCCCAGAAACATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGcc                                                                                                                                                                                                                                                   <  2:518058/300‑1 (MQ=255)
                                                                    tccgatcgCCGCGCCGGTATTCCAGATGGGATAACCGCCCGGTGGTTTAGCCGTCGCGCTGTCGGCGCGTTTCTGCCCGGTATAGTGCCCAGAAACATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACg                                                                                                                                                                                                                             <  1:383993/294‑1 (MQ=255)
                                                                            ccGCGCCGGTATTCCAGATGGTATAACCGCCCGGTGTTTTAGCCGTCGCGCTGTCGGCGCGTTTCTGCCCGGTATAGTGCCCAGAAACATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACgagatcgga                                                                                                                                                                                                                     >  2:384439/1‑293 (MQ=255)
                                                                            ccGCGCCGGTATTCCAGATGGTATAACCGCCCGGTGTTTTAGCCGTCGCGCTGTCGGCGCGTTTCTGCCCGGTATAGTGCCCAGAAACATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTgtagta                                                                                                                                                                                                                     >  2:358398/1‑301 (MQ=255)
                                                                                      aTTCCAGATGGTATAACCGCCCGGTGTTTTAGCCGTCGCGCTGTCGGCGCGTTTCTGCCCGGTATAGTGCCCAGAAACATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAcc                                                                                                                                                                                                           <  2:586033/301‑1 (MQ=255)
                                                                                                         cccGGTGTTTTCGCTGTCGCGCTGTCGGCGCGTTTCCGCCCGGTATAGTGCCCAGAAACATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCATACCGTTAGCAGTATGGAGCGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAg                                                                                                                                                                                        <  2:203195/301‑1 (MQ=255)
                                                                                                          ccGGTGTTTTAGACGTCGGGCGGGCGGCGCGTTTCGTCCCGGTATAGGGCCCAGAAACATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAgc                                                                                                                                                                                       <  1:583059/301‑1 (MQ=255)
                                                                                                             gtgtTTTAGCCGTCGCGCTGTCGGCGCGTTTCTGCCCGGTATAGTGCCCAGAAACATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAgcgc                                                                                                                                                                                     <  1:537676/300‑1 (MQ=255)
                                                                                                               gtTTTAGCCGTCGCGCTGTCGGCGCGTTTCTGCCCGGTATAGTGCCCAGAAACATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCg                                                                                                                                                                                  >  2:565030/1‑301 (MQ=255)
                                                                                                                       cGTCGCGCTGTCGGCGCGTTTCTGCCCGGTATAGTGCCCAGAAACATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCCCGACCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCaaaa                                                                                                                                                                          >  2:420054/1‑301 (MQ=255)
                                                                                                                       cGTCGCGCTGTCGGCGCGTTTCTGCCCGGTATAGTGCCCAGAAACATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCaaaa                                                                                                                                                                          <  2:68615/301‑1 (MQ=255)
                                                                                                                                gTCGGCGCGTTTCTGCCCGGTATAGTGCCCAGAAACATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAg                                                                                                                                                                 >  2:265794/1‑301 (MQ=255)
                                                                                                                                 tCGGCGCGTTTCTTCCCGGTATAGTGCCCAGAAACATCGAAGGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAg                                                                                                                                                                 <  1:556335/300‑1 (MQ=255)
                                                                                                                                       cgTTTCTGCCCGGTATAGTGCCCAGAAACATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCCCAAACCCAACAAAGTTTTTGt                                                                                                                                                          >  1:227000/1‑301 (MQ=255)
                                                                                                                                                       aGTGCCCAGAAACATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGtt                                                                                                                                          >  1:341248/1‑301 (MQ=255)
                                                                                                                                                                aaaCATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGACTACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAACCCACCAAAGTTTTTGTAGCCCGGTGCAGCGTTGACGTCAGa                                                                                                                                 >  1:399975/1‑301 (MQ=255)
                                                                                                                                                                        aaTGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACgg                                                                                                                          <  2:136236/300‑1 (MQ=255)
                                                                                                                                                                                    tCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATa                                                                                                              >  2:306965/1‑300 (MQ=255)
                                                                                                                                                                                          gccagcgGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCtt                                                                                                       >  1:135738/1‑301 (MQ=255)
                                                                                                                                                                                                              cGTACCGTTAGCAGTATTGAACGGCAGATCGGATAGCCGTTTGTTTTCGCCTTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGTTACGGaaa                                                                                      <  1:265481/298‑1 (MQ=255)
                                                                                                                                                                                                                 cccGTTAGCAGTATGGAACGGCATGGTGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGTTACGGAAAACGGTa                                                                                <  1:74600/300‑1 (MQ=255)
                                                                                                                                                                                                                      tAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGCCCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTGTACGGAAGTTTCAGTTCGGTTTCCCCGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCGTTTCGCACCCGTCTCAAAACCAACAAAGTTTTCGTAGCCCGTTGCAGCGTTGCCGTCCGACTTACGGCTGCTGTGGCTACGTTAGTTCACCTCGTTCCGCAAATCGGTTACGc                                                                            >  1:401051/1‑300 (MQ=255)
                                                                                                                                                                                                                             aTGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGTTACGGAAAACGGTAACGCTGGATTCa                                                                    <  2:341625/301‑1 (MQ=255)
                                                                                                                                                                                                                                gagcGGCAGATCTGATAGCGGTTTTTTTTCGCCTTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGTTACGGAAAACGGTAACGCTGGATTCAAc                                                                  <  1:306620/297‑1 (MQ=255)
                                                                                                                                                                                                                                gAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGTTACGGAAAACGGTAACGCTGGATTCAAc                                                                  >  1:336339/1‑300 (MQ=255)
                                                                                                                                                                                                                                 aaCGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGTTACGGAAAACGGTAACGCTGGATTCAAccc                                                                <  2:135901/301‑1 (MQ=255)
                                                                                                                                                                                                                                     gCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGTTACGGAAAACGGTAACGCTGGATTCAACCCCtt                                                             >  1:222850/1‑300 (MQ=255)
                                                                                                                                                                                                                                                     tttgtttTCGCCGTTGCTGACATCCCGACCGTCGTTGTAGGTGTAGTTGTTCGACAGTGTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGGGTTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGTTACGGAAAACGGTAACGCTGGATTCAACCCCTTCCAGCCAGCCTTCTTCa                                            <  1:48480/301‑1 (MQ=255)
                                                                                                                                                                                                                                                         ttttCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGTTACGGAAAACGGTAACGCTGGATTCAACCCCTTCCAGCCAGCCTTCTTCACCCa                                        <  2:435526/301‑1 (MQ=255)
                                                                                                                                                                                                                                                          tttCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGTTACGGAAAACGGTAACGCTGGATTCAACCCCTTCCAGCCAGCCTTCTTCACCCAt                                       >  1:97696/1‑301 (MQ=255)
                                                                                                                                                                                                                                                          tttCGCCGTTGCTGACATCACGACCATCGTTGTAGGAGTAGTTGTTCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCTTGAATATGAGCTTTGTTAACGTTGTAGTAGCTAATTACCGCTATGCGCCGTCCGTTAGCACACGCCTCAAAACCCACAAAGTTTTGATAGCTCGGTGCAGAGTATACGTCAGAAGTTCGGCGGATGCTGATCCGACCGTTCACAGCGTTACGGAAACAGGCAACTCAGGATTCTACCCCtactcgacggcactcctccaccat                                       >  1:156116/1‑278 (MQ=255)
                                                                                                                                                                                                                                                           ttCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGTTACGGAAAACGGTAACGCTGGATTCAACCCCTTCCAGCCAGCCTTCTTCACCCATg                                      <  2:227282/301‑1 (MQ=255)
                                                                                                                                                                                                                                                           ttCGCCGTGGCTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGTTACGGAAAACGGTAACGCTGGATTCAACCCCTTCCAGCCAGCCTTCTTCACCCATg                                      <  1:419559/301‑1 (MQ=255)
                                                                                                                                                                                                                                                              gCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGTTACGGAAAACGGTAACGCTGGATTCAACCCCTTCCAGCCAGCCTTCTTCACCCATgtag                                   >  2:380240/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                     cTGACATCACGACCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGTTACGGAAAACGGTAACGCTGGCTTCAACCCCTTCCAGCCAGCCTTCTTCACACATGTAGTCAAGc                             >  1:42215/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                              cGACCATCGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGTTACGGAAAACGGTAACGCTGGATTCAACCCCTTCCAGCCAGCCTTCTTCACCCATGTAGTAAAGCCCCAGCTccc                   <  2:73500/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                     cGTTGTAGGTGTAGTTGTTCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGTTACGGAAAACGGTAACGCTGGATTCAACCCCTTCCAGCCAGCCTTCTTCACCCATGTAGTAAAGCCCCAGCTCCCAACTTTc            <  2:336716/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                tAGTTGTTCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGTTACGGAAAACGGTAACGCTGGATTCAACCCCTTCCAGCCAGCCTTCTTCACCCATGTAGTAAAGCCCCAGCTCCCAACTTTCGctggtttct  <  1:357997/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                                     
TCGTCACGACTGAGATCCTTGTCGCCAAGGTTCAGCACGCCTGCGCGCAGTTTGACGTCTTTAGTCACCTGCCAGGCCGCGCCGGTATTCCAGATGGTATAACCGCCCGGTGTTTTAGCCGTCGCGCTGTCGGCGCGTTTCTGCCCGGTATAGTGCCCAGAAACATAGAATGACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGTACGGCTGATGCTGATACGATCTTTCACATCGTTACGGAAAACGGTAACGCTGGATTCAACCCCTTCCAGCCAGCCTTCTTCACCCATGTAGTAAAGCCCCAGCTCCCAACTTTCGCTGGTTTCT  >  NC_000913/2244836‑2245423

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑