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.

AAACAGAAAAAACGTTAACTATTCGCGATTCAGTTTTTACTTATACCGAGAATGCTGACGGGACTATATCTCTGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTG‑ATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGA‑TATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACC  >  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)
                                                           gggACTATATCTCTGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGAAAGGATTCCACCGGTACGGAAATCAACGGCAATAACGGGAAGGGG‑ATTc                                                                                                                                                                                               >  2:341735/1‑250 (MQ=255)
                                                                tataTCTCTGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGCAATACGGGGAAAGTG‑ATTCAgggcg                                                                                                                                                                                         >  2:390685/1‑248 (MQ=255)
                                                                 atatCTCTGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTG‑ATTCAGGACgg                                                                                                                                                                                        <  1:411565/251‑1 (MQ=255)
                                                                           aaGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTG‑ATTCAGGACGGCGATCTGGAt                                                                                                                                                                              >  1:113754/1‑251 (MQ=255)
                                                                                 gCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCATATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTG‑ATTCAGGACGGCGATCTGGATGTCAgc                                                                                                                                                                        <  2:231886/251‑1 (MQ=255)
                                                                                 gCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTG‑ATTCAGGACGGCGATCTGGATGTCAgc                                                                                                                                                                        >  1:91385/1‑251 (MQ=255)
                                                                                  cAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAAtt                                                                                                                                                                                                                         >  2:368947/1‑201 (MQ=255)
                                                                                  cAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAAtt                                                                                                                                                                                                                         <  1:368947/202‑2 (MQ=255)
                                                                                   aaTGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGCAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTG‑ATTCAGGACGGCGATCTGGATGTCAgcgg                                                                                                                                                                      >  1:476795/1‑251 (MQ=255)
                                                                                   aaTGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTG‑ATTCAGGACGGCGATCTGGATGTCAgcgg                                                                                                                                                                      >  1:161957/1‑251 (MQ=255)
                                                                                                    gATTAATCTTTGCCAGCCTGATGAAACGCCTACCCCTGCTGCCCTTGACCGGGCGAGCGCAGTTGGCGCCCCGAAGTGGCAATATACTCACACCGGCGAGCTTGTTCTTACGGGTGATAATGCCACAGTAAACAACAATGGCCAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCCACGGTAATCCCGGGACAGTG‑ATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGACCACGGGATTg                                                                                                                                                     >  2:458016/1‑251 (MQ=16)
                                                                                                                      tgatgaAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTG‑ATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGA‑TATCACCGGTGACAgc                                                                                                                                   <  2:113754/250‑1 (MQ=255)
                                                                                                                      tgatgaAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTG‑ATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGA‑TATCACCGGTGACAgcg                                                                                                                                  >  1:138582/1‑251 (MQ=255)
                                                                                                                           aaGCGACTAACAATGTTGACCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGCAATAAAGGGAAAGTG‑ATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGA‑TATCACCGGCGACAGCGCCACg                                                                                                                             >  2:300316/1‑251 (MQ=25)
                                                                                                                           aaGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCCGATGGCGCAACGAAGTGGCAATATAATCACAAAGGTGAGCTTGTTATTACGGTTGATAATGCCACAGTAAACAACAATGGAAAAACCACCGTTGACGGTAAAGATTCCACAGGTACGGAAATCAACGGTAATAACGGGAAAGTG‑ATTCAGGACGGCGATATGGATGTCAGCGGCGCAGGTAACGGTATTGA‑TATCACCGGTGACAGCGCTACg                                                                                                                             >  1:444330/1‑251 (MQ=255)
                                                                                                                            aGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTG‑ATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGA‑TATCACCGGTGACAGCGCGACgg                                                                                                                            >  1:354506/1‑251 (MQ=255)
                                                                                                                               gAATAACACTGTTGCCCTTGAAGGGGAGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGCGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTCAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAACGTG‑ATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGGCACGGTATTGA‑TCTCACCGGGGAACGCGCGACGGTgg                                                                                                                         >  2:325693/1‑251 (MQ=25)
                                                                                                                                aaTAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTG‑ATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGA‑TATCACCGGTGACAGCGCGACGGTgg                                                                                                                         >  1:205264/1‑250 (MQ=255)
                                                                                                                                  tAACAATGTTGCCATTGCAGGGGTGTGCGCAGCTGGCGCTAAGAAGTGGCAATATAGTCACAACGGTGAGCTCATTATTACGGGTGATAATGCCACAGTAAACAACACTGGCACAACCACTGTTGACGGTAAAGATTCCACCGGTACGGCAATCAACGGTATTCACGGGAAAGTG‑ATTCAGGACGGCGCTCTGGATGTCAGCGGCTGCGGGCACGGGCTTGA‑TTAGACCGGGGACAGGGCGACGGTgggt                                                                                                                       >  2:476993/1‑248 (MQ=255)
                                                                                                                                            gCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTG‑ATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGA‑TATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCAc                                                                                                            >  1:70207/1‑251 (MQ=255)
                                                                                                                                                  gAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTG‑ATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGA‑TATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGAc                                                                                                      >  1:92446/1‑251 (MQ=255)
                                                                                                                                                    aGGGGTGAGCGCAGCTGGCGCAACGCAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTG‑ATTCAGGACGGCGATCTGGCTGTCAGCGGCAGAGGGGACGGTATTCA‑TATCACGGGTGAGAGCGCGACGGGGGGGACCAAGGGCACCATGaac                                                                                                     >  2:450847/1‑248 (MQ=25)
                                                                                                                                                         ggAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTG‑ATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGA‑TATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCAcc                                                                                               <  2:19068/250‑1 (MQ=255)
                                                                                                                                                          gAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTG‑ATTCAGGACGGAGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGA‑TATCACCGGTGACAGCGCGACGGGGGATAACAAGGGCAACATGACCGTCCCCg                                                                                              >  2:168486/1‑251 (MQ=255)
                                                                                                                                                             cgcAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATACCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTG‑ATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGA‑TATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACTGTCACCGATc                                                                                           >  1:460191/1‑251 (MQ=255)
                                                                                                                                                               aaaaTGGAGCAACGAAGTGGCAAAACAACCGCAACGGCGAGACTGTTATTACGGGTGATAATGCCACAGCAAACAACAATGGCAAAACCACTGTTGACGGTAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTT‑ATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGA‑TATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCg                                                                                         <  2:381248/248‑1 (MQ=11)
                                                                                                                                                                    ggCGCAACGAAGTGGCAATATAATCACAACGTTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTG‑ATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGA‑TATCACCGGTGACAGCGCGACGGGGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGt                                                                                     <  1:306475/250‑1 (MQ=255)
                                                                                                                                                                         aaCGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTG‑ATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGA‑TATCACCGGTGACAGCGCGACGGTGGGTAACAAGGGCACCATGACCGTCACCCATCCGGAGTCCacgg                                                                               >  1:279126/1‑248 (MQ=255)
                                                                                                                                                                              aGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGAAAAACCACCGTTGCCGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTG‑ATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGA‑TATCACCGGCGCCAGCGCAACGGTGGATAACAAGGGCACCATGAACGTCACCGATCCGGAGtaccgcggtatc                                                                          >  2:466441/1‑240 (MQ=18)
                                                                                                                                                                                    aaTATAATCACAAAGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTAGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTG‑ATTCAGGACGGCGATCTGGCTGTCAGCGGCGGCGGGCACGGGATTGA‑TATCACCGGCGACAGCGCGACGGGGGATCACAAGGGCACCCTGACCGTCACCGATCCGGAATCCATGGGGATCCAggtg                                                                    >  2:227221/1‑248 (MQ=18)
                                                                                                                                                                                       ataATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTG‑ATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGA‑TATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCGTGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGAc                                                                 <  1:446346/251‑1 (MQ=255)
                                                                                                                                                                                           tCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTG‑ATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGA‑TATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATAGACGGTg                                                             >  1:262198/1‑251 (MQ=255)
                                                                                                                                                                                              caACGGCGAGCTTGTTATTACGGGTGATAATGCCACAGTCACCAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCCCCGGTACGGAAATCAACGGTAATAACGGGACCGTG‑ATTATCGACGGCGATCTGGCTGTAACCGGCGGCGCACAAGGAATTCA‑TATCACCGGTTACGCCGCGACGGTGGCTAACAAGGGCACCATGCCCGACACCGATCCGGCCACCCTGGGCATCCAGATCGACGGCGAc                                                           >  1:69526/1‑250 (MQ=255)
                                                                                                                                                                                                aCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTG‑ATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGA‑TATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAg                                                        <  2:92446/251‑1 (MQ=255)
                                                                                                                                                                                                  ggTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTG‑ATTCAGGACGGCGATCTGGATGTCAGCGACGACACGCACGGAATTGA‑TATCCCCGGATACAGCGCGACGGGGCTTAACAATGGCACAATGACCGTCACCGATCCGGAGACCATGGGTATCCAGATCGACGGTGACAAGGc                                                      >  1:160715/1‑251 (MQ=18)
                                                                                                                                                                                                   gTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTG‑ATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGA‑TATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGcc                                                     <  1:168486/251‑1 (MQ=255)
                                                                                                                                                                                                      aGCTTGTTATTACGGGTGATAATGCCACAGTAAACCACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGCAATAACGGGAAAGTGCGTTC‑GGACGGCCGTCTGGATGTCAGCGGCGGCGGACCCGGGATTGA‑GACGAACCGTTACAGAGCGACGGTGGATAACAAGGGCACTATGACAGCCACCACTACGGAGTACATGGGAATCCAGGTACACCGTGACAAGGccctc                                                  >  2:29878/1‑248 (MQ=255)
                                                                                                                                                                                                                 aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATTAACGGTAATAACGGGAAAGTT‑CTTCAGAGCGGCGCTCTGAATGTCAGCGACGGCTGTCACGGGATTGC‑TGTCTCCGGTGTCAGCGCGACGGTGGGTAACAAGGGCACCATGACCGTCATCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATatcg                                        >  1:337874/1‑247 (MQ=255)
                                                                                                                                                                                                                 aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGATGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTG‑ATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGA‑TATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACg                                        >  1:288321/1‑250 (MQ=255)
                                                                                                                                                                                                                  cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTG‑ATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGA‑TATCACCGGCGACAGCGCGCCGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGa                                       >  1:213716/1‑250 (MQ=255)
                                                                                                                                                                                                                     gTGATAATGCCACAGTAAACAACAATGGCAAAACCACTGGTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTG‑ATTCAGGACGGCGAGCTGGATGTCAGCGGCGGCGGTCACGGTATTGA‑TATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATTTGTTTCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGc                                   <  1:470537/251‑1 (MQ=255)
                                                                                                                                                                                                                             gCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTG‑ATTCAGGACGGCGATCTGGATGTCAGAGGCGGCGGTCACGGAATTGA‑TCTCACCGGAGAAAGCCGGATGCCGGATAACAAGGGCACCAAGACCGTAACCGACACGGAGTCCATGGGAATGCAGATCGCCGGTGACAAGGACGTCGTCCATAACTAAGGCGAGAGCAc                           >  1:274573/1‑251 (MQ=17)
                                                                                                                                                                                                                                   gTAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTG‑ATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGA‑TATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCAc                     <  1:64724/251‑1 (MQ=255)
                                                                                                                                                                                                                                    tAAACAACAATGGCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTG‑ATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGACACGGTATTGA‑TATCACCGGGGACAGCGCGACGGTGGATACCAAGGGCACCCTTACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGGGACAAGGCCATCGACAATAACGAAGGCGAGAGCACCATCAcc                    >  1:385055/1‑251 (MQ=17)
                                                                                                                                                                                                                                             aTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTG‑ATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGG‑AATCGCCGGGGGCCGCGGGGGGGAGGAGAGCAGGGGCACGTTGTCCGCGACCGATCCGGAGTACATGGGGATCCAGATCGACGGGGACAAGGCACTCGTCAATAACGAAGGCGAGGGAACCATCACCAACggcgg            >  2:317565/1‑247 (MQ=12)
                                                                                                                                                                                                                                             aTGGCAAAACAACCGTTGACGGAAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTG‑ATTCAGGGCGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGA‑TATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggc           >  1:370781/1‑251 (MQ=255)
                                                                                                                                                                                                                                               ggCAAAACCACTGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTG‑ATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGA‑TATCACCGGTGACAGCGCGACGTTGGAGAACAAGGGc                                                                                                              <  1:383538/150‑1 (MQ=255)
                                                                                                                                                                                                                                               ggCAAAACCACTGTTGAAGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTG‑ATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGA‑TATCACCGGTGACAGCGCGACGGTGGATAACAAGGGc                                                                                                              >  2:383538/1‑150 (MQ=255)
                                                                                                                                                                                                                                                   aaaCCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTG‑ATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGCATTGA‑TATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGAGCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCACCGGCACTAACGAAGGCGAGAGCACCACCACCACCGGTggcacccgc     >  2:260066/1‑248 (MQ=18)
                                                                                                                                                                                                                                                      ccaccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTG‑ATTCAGGACGGCGATCCGGATGTCCGCGGCGGAGGGCACGGGATTGAATTTC‑GCGGTTCCCCGGAGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCCTCGTCAATAACGCAGGCGAGAGCACCCTCACCAACGGTggcaccggcacc  >  2:196056/1‑251 (MQ=18)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                        
AAACAGAAAAAACGTTAACTATTCGCGATTCAGTTTTTACTTATACCGAGAATGCTGACGGGACTATATCTCTGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTG‑ATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGA‑TATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACC  >  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: ‑