Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,466,0900CT53.2% 17.9 / 83.0 62pseudogene (699/2518 nt)ydbAputative outer membrane protein, N‑terminal fragment
Reads supporting (aligned to +/- strand):  ref base C (26/3);  new base T (20/13);  total (46/16)
Fisher's exact test for biased strand distribution p-value = 1.05e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.51e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AAACAGAAAAAACGTTAACTATTCGCGATTCAGTTTTTACTTATACCGAGAATGCTGACGGGACTATATCTCTGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACC  >  NC_000913/1465840‑1466336
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                        
aaaCAGAAAAAACGTTAACTATTCGCGATTCAGTTTTTACTTATACCGAGAATGCTGACGGGACTATATCTCTGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAaccacc                                                                                                                                                                                                                                                          <  1:268354/251‑1 (MQ=255)
   cAGAAAAAACGTTAACTATTCGCGATTCAGTTTTTACTTATACCGAGAATGCTGACGGGACTATATCTCTGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGtt                                                                                                                                                                                                                                                       <  2:149914/251‑1 (MQ=255)
                   tATTCGCGATTCAGTTTTTACTTATACCGAGAATGCTGACGGGACTATATCTCTGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTcc                                                                                                                                                                                                                                        >  1:19068/1‑250 (MQ=255)
                        gcgATTCAGTTTTTACTTATACCGAGAATGCTGACGGGACTATATCTCTGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGt                                                                                                                                                                                                                                  >  1:237687/1‑251 (MQ=255)
                         cgATTCAGTTTTTACTTATACCGAGAATGCTGACGGGACTATATCTCTGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTa                                                                                                                                                                                                                                 <  2:221663/251‑1 (MQ=255)
                          gATTCAGTTTTTACTTATACCGAGAATGCTGACGGGACTATATCTCTGCAAGATAGCACTGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                                                                                                                                                                                                >  2:87089/1‑251 (MQ=255)
                          gATTCAGTTTTTACTTATACCGAGAATGCTGACGGGACTATATCTCTGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGTCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                                                                                                                                                                                                >  2:20919/1‑251 (MQ=255)
                          gATTCAGTTTTTACTTATACCGAGAATGCTGACGGGACTATATCTCTGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTAc                                                                                                                                                                                                                                >  1:312483/1‑251 (MQ=255)
                                  ttttACTTATACCGAGAATGCTGACGGGACTATATCTCTGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTGTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATc                                                                                                                                                                                                                         <  2:456743/250‑1 (MQ=255)
                                               gagaATGCTGACGGGACTATATCTCTGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACgg                                                                                                                                                                                                           >  1:381248/1‑251 (MQ=255)
                                                agaATGCTGACGGGACTATATCTCTGCAAGATAGCACTGGTCGTAACGCAACGATTAATCTTTGGCAGATTGGTGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGCGCTTGTTATTACGGGTGATAATGCCACAGTAAACCACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACggg                                                                                                                                                                                                          >  1:455840/1‑251 (MQ=255)
                                                  aaTGCTGACGAGACTATCTCGCTGCAAGATAGAAATGGTCGTAAGGCAACGATTCATCTTTGGCCGATTGATGAAGCGAATAACACTGCTGCCCTTGACGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACGACGGCAAAACCACTGTTGACGGTAAAGCTTCCACCGTTAAGGAAATCAACGGTCATAACGGGa                                                                                                                                                                                                         >  2:169768/1‑250 (MQ=255)
                                                   aTGCTGACGGGACTATATCTCTGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaa                                                                                                                                                                                                        >  2:265328/1‑250 (MQ=255)
                                                     gCTGACGGGACTATATCTCTGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCCCCGGTACGGAAATCAACGGTAATAACGGGAAAGt                                                                                                                                                                                                     >  2:336046/1‑251 (MQ=255)
                                                           gggACTATATCTCTGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGAAAGGATTCCACCGGTACGGAAATCAACGGCAATAACGGGAAGGGGATTc                                                                                                                                                                                               >  2:341735/1‑250 (MQ=255)
                                                                tataTCTCTGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGCAATACGGGGAAAGTGATTCAgggcg                                                                                                                                                                                         >  2:390685/1‑248 (MQ=255)
                                                                 atatCTCTGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACgg                                                                                                                                                                                        <  1:411565/251‑1 (MQ=255)
                                                                           aaGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGAt                                                                                                                                                                              >  1:113754/1‑251 (MQ=255)
                                                                                 gCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCATATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAgc                                                                                                                                                                        <  2:231886/251‑1 (MQ=255)
                                                                                 gCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAgc                                                                                                                                                                        >  1:91385/1‑251 (MQ=255)
                                                                                  cAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAAtt                                                                                                                                                                                                                         >  2:368947/1‑201 (MQ=255)
                                                                                  cAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAAtt                                                                                                                                                                                                                         <  1:368947/202‑2 (MQ=255)
                                                                                   aaTGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGCAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAgcgg                                                                                                                                                                      >  1:476795/1‑251 (MQ=255)
                                                                                   aaTGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAgcgg                                                                                                                                                                      >  1:161957/1‑251 (MQ=255)
                                                                                                    gATTAATCTTTGCCAGCCTGATGAAACGCCTACCCCTGCTGCCCTTGACCGGGCGAGCGCAGTTGGCGCCCCGAAGTGGCAATATACTCACACCGGCGAGCTTGTTCTTACGGGTGATAATGCCACAGTAAACAACAATGGCCAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCCACGGTAATCCCGGGACAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGACCACGGGATTg                                                                                                                                                     >  2:458016/1‑251 (MQ=16)
                                                                                                                      tgatgaAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAgc                                                                                                                                   <  2:113754/250‑1 (MQ=255)
                                                                                                                      tgatgaAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAgcg                                                                                                                                  >  1:138582/1‑251 (MQ=255)
                                                                                                                           aaGCGACTAACAATGTTGACCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGCAATAAAGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCCACg                                                                                                                             >  2:300316/1‑251 (MQ=25)
                                                                                                                           aaGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCCGATGGCGCAACGAAGTGGCAATATAATCACAAAGGTGAGCTTGTTATTACGGTTGATAATGCCACAGTAAACAACAATGGAAAAACCACCGTTGACGGTAAAGATTCCACAGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATATGGATGTCAGCGGCGCAGGTAACGGTATTGATATCACCGGTGACAGCGCTACg                                                                                                                             >  1:444330/1‑251 (MQ=255)
                                                                                                                            aGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACgg                                                                                                                            >  1:354506/1‑251 (MQ=255)
                                                                                                                               gAATAACACTGTTGCCCTTGAAGGGGAGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGCGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTCAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAACGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGGCACGGTATTGATCTCACCGGGGAACGCGCGACGGTgg                                                                                                                         >  2:325693/1‑251 (MQ=25)
                                                                                                                                aaTAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTgg                                                                                                                         >  1:205264/1‑250 (MQ=255)
                                                                                                                                  tAACAATGTTGCCATTGCAGGGGTGTGCGCAGCTGGCGCTAAGAAGTGGCAATATAGTCACAACGGTGAGCTCATTATTACGGGTGATAATGCCACAGTAAACAACACTGGCACAACCACTGTTGACGGTAAAGATTCCACCGGTACGGCAATCAACGGTATTCACGGGAAAGTGATTCAGGACGGCGCTCTGGATGTCAGCGGCTGCGGGCACGGGCTTGATTAGACCGGGGACAGGGCGACGGTgggt                                                                                                                       >  2:476993/1‑248 (MQ=255)
                                                                                                                                            gCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCAc                                                                                                            >  1:70207/1‑251 (MQ=255)
                                                                                                                                                  gAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGAc                                                                                                      >  1:92446/1‑251 (MQ=255)
                                                                                                                                                    aGGGGTGAGCGCAGCTGGCGCAACGCAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGCTGTCAGCGGCAGAGGGGACGGTATTCATATCACGGGTGAGAGCGCGACGGGGGGGACCAAGGGCACCATGaac                                                                                                     >  2:450847/1‑248 (MQ=25)
                                                                                                                                                         ggAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCAcc                                                                                               <  2:19068/250‑1 (MQ=255)
                                                                                                                                                          gAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGAGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGGGGATAACAAGGGCAACATGACCGTCCCCg                                                                                              >  2:168486/1‑251 (MQ=255)
                                                                                                                                                             cgcAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATACCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACTGTCACCGATc                                                                                           >  1:460191/1‑251 (MQ=255)
                                                                                                                                                               aaaaTGGAGCAACGAAGTGGCAAAACAACCGCAACGGCGAGACTGTTATTACGGGTGATAATGCCACAGCAAACAACAATGGCAAAACCACTGTTGACGGTAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCg                                                                                         <  2:381248/248‑1 (MQ=11)
                                                                                                                                                                    ggCGCAACGAAGTGGCAATATAATCACAACGTTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGGGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGt                                                                                     <  1:306475/250‑1 (MQ=255)
                                                                                                                                                                         aaCGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGGTAACAAGGGCACCATGACCGTCACCCATCCGGAGTCCacgg                                                                               >  1:279126/1‑248 (MQ=255)
                                                                                                                                                                              aGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGAAAAACCACCGTTGCCGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGCCAGCGCAACGGTGGATAACAAGGGCACCATGAACGTCACCGATCCGGAGtaccgcggtatc                                                                          >  2:466441/1‑240 (MQ=18)
                                                                                                                                                                                    aaTATAATCACAAAGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTAGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGCTGTCAGCGGCGGCGGGCACGGGATTGATATCACCGGCGACAGCGCGACGGGGGATCACAAGGGCACCCTGACCGTCACCGATCCGGAATCCATGGGGATCCAggtg                                                                    >  2:227221/1‑248 (MQ=18)
                                                                                                                                                                                       ataATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCGTGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGAc                                                                 <  1:446346/251‑1 (MQ=255)
                                                                                                                                                                                           tCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATAGACGGTg                                                             >  1:262198/1‑251 (MQ=255)
                                                                                                                                                                                              caACGGCGAGCTTGTTATTACGGGTGATAATGCCACAGTCACCAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCCCCGGTACGGAAATCAACGGTAATAACGGGACCGTGATTATCGACGGCGATCTGGCTGTAACCGGCGGCGCACAAGGAATTCATATCACCGGTTACGCCGCGACGGTGGCTAACAAGGGCACCATGCCCGACACCGATCCGGCCACCCTGGGCATCCAGATCGACGGCGAc                                                           >  1:69526/1‑250 (MQ=255)
                                                                                                                                                                                                aCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAg                                                        <  2:92446/251‑1 (MQ=255)
                                                                                                                                                                                                  ggTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGACGACACGCACGGAATTGATATCCCCGGATACAGCGCGACGGGGCTTAACAATGGCACAATGACCGTCACCGATCCGGAGACCATGGGTATCCAGATCGACGGTGACAAGGc                                                      >  1:160715/1‑251 (MQ=18)
                                                                                                                                                                                                   gTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGcc                                                     <  1:168486/251‑1 (MQ=255)
                                                                                                                                                                                                      aGCTTGTTATTACGGGTGATAATGCCACAGTAAACCACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGCAATAACGGGAAAGTGCGTTCGGACGGCCGTCTGGATGTCAGCGGCGGCGGACCCGGGATTGAGACGAACCGTTACAGAGCGACGGTGGATAACAAGGGCACTATGACAGCCACCACTACGGAGTACATGGGAATCCAGGTACACCGTGACAAGGccctc                                                  >  2:29878/1‑248 (MQ=255)
                                                                                                                                                                                                                 aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTAACGGTAATAACGGGAAAGTTCTTCAGAGCGGCGCTCTGAATGTCAGCGACGGCTGTCACGGGATTGCTGTCTCCGGTGTCAGCGCGACGGTGGGTAACAAGGGCACCATGACCGTCATCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATatcg                                        >  1:337874/1‑247 (MQ=255)
                                                                                                                                                                                                                 aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGATGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACg                                        >  1:288321/1‑250 (MQ=255)
                                                                                                                                                                                                                  cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCGCCGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGa                                       >  1:213716/1‑250 (MQ=255)
                                                                                                                                                                                                                     gTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGGTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGAGCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATTTGTTTCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGc                                   <  1:470537/251‑1 (MQ=255)
                                                                                                                                                                                                                             gCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGAGGCGGCGGTCACGGAATTGATCTCACCGGAGAAAGCCGGATGCCGGATAACAAGGGCACCAAGACCGTAACCGACACGGAGTCCATGGGAATGCAGATCGCCGGTGACAAGGACGTCGTCCATAACTAAGGCGAGAGCAc                           >  1:274573/1‑251 (MQ=17)
                                                                                                                                                                                                                                   gTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCAc                     <  1:64724/251‑1 (MQ=255)
                                                                                                                                                                                                                                    tAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGACACGGTATTGATATCACCGGGGACAGCGCGACGGTGGATACCAAGGGCACCCTTACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGGGACAAGGCCATCGACAATAACGAAGGCGAGAGCACCATCAcc                    >  1:385055/1‑251 (MQ=17)
                                                                                                                                                                                                                                             aTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGGAATCGCCGGGGGCCGCGGGGGGGAGGAGAGCAGGGGCACGTTGTCCGCGACCGATCCGGAGTACATGGGGATCCAGATCGACGGGGACAAGGCACTCGTCAATAACGAAGGCGAGGGAACCATCACCAACggcgg            >  2:317565/1‑247 (MQ=12)
                                                                                                                                                                                                                                             aTGGCAAAACAACCGTTGACGGAAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGGCGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggc           >  1:370781/1‑251 (MQ=255)
                                                                                                                                                                                                                                               ggCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGTTGGAGAACAAGGGc                                                                                                              <  1:383538/150‑1 (MQ=255)
                                                                                                                                                                                                                                               ggCAAAACCACTGTTGAAGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGc                                                                                                              >  2:383538/1‑150 (MQ=255)
                                                                                                                                                                                                                                                   aaaCCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGCATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGAGCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCACCGGCACTAACGAAGGCGAGAGCACCACCACCACCGGTggcacccgc     >  2:260066/1‑248 (MQ=18)
                                                                                                                                                                                                                                                      ccaccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCCGGATGTCCGCGGCGGAGGGCACGGGATTGAATTTCGCGGTTCCCCGGAGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCCTCGTCAATAACGCAGGCGAGAGCACCCTCACCAACGGTggcaccggcacc  >  2:196056/1‑251 (MQ=18)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                        
AAACAGAAAAAACGTTAACTATTCGCGATTCAGTTTTTACTTATACCGAGAATGCTGACGGGACTATATCTCTGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACC  >  NC_000913/1465840‑1466336

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