Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,466,2010TC32.6% 53.4 / 24.0 45pseudogene (810/2518 nt)ydbAputative outer membrane protein, N‑terminal fragment
Reads supporting (aligned to +/- strand):  ref base T (17/12);  new base C (6/8);  total (25/20)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACG  >  NC_000913/1465960‑1466448
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       
atgaAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGCAATAACGGAAAAGCGATTCAGGACGGCGATCTGGATGTCAGCGGCGTCGGACGCGGTCTTGATATCACCGGGGACAGcgcg                                                                                                                                                                                                                                                >  1:147036/1‑251 (MQ=38)
      ggAATAACACTGTTGCCGATGAAGGGGTGAGCACAAATGGAGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCACCGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTg                                                                                                                                                                                                                                          <  2:2742/250‑1 (MQ=21)
       gAATAACACTGTTGCCCTTGAAGGGGTGGTCGCTGATGGCGGAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTgg                                                                                                                                                                                                                                         <  2:173996/251‑1 (MQ=17)
          tAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGAt                                                                                                                                                                                                                                       >  2:108598/1‑250 (MQ=255)
             cacTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGAAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACa                                                                                                                                                                                                                                   <  2:254848/251‑1 (MQ=255)
                      ccTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCAcc                                                                                                                                                                                                                           >  1:335462/1‑250 (MQ=255)
                         tGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATg                                                                                                                                                                                                                        >  1:120279/1‑250 (MQ=255)
                           aaGGGGTGATCACAAATGGCGCGACGAAGTGGAAATATAATCACAACGGCGAGATTGTGATTACGGGTGATAATGCCACAGTAAACAACAACGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAGTGATATCACCGGCGACAGCGCAACGGTTGATAACAAGGGCACCTTGAcc                                                                                                                                                                                                                     <  2:48944/251‑1 (MQ=17)
                                       cAGATGGCGGAACGAAGTGGCAATATAATCGCAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCGCTGGCAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCg                                                                                                                                                                                                         <  1:108598/251‑1 (MQ=17)
                                                      agaggcaaactaaccacgacagTGAGCATGTGATTACGGGTGAGAATGCCACAGTAAAAAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATc                                                                                                                                                                                          <  2:7270/230‑1 (MQ=14)
                                                      aGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCGACGGTGGATAACAAGGGCACCATGACCgacatcggtagggcgtcgtggggtga                                                                                                                                                                                           >  2:39326/1‑225 (MQ=18)
                                                           cAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGGGGCTAACAAGGGCACCACGACAGTCACCGATCCGGAGTCCATGGGTATcc                                                                                                                                                                                         >  1:359038/1‑247 (MQ=255)
                                                            aaTATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATcc                                                                                                                                                                                         >  1:108063/1‑246 (MQ=255)
                                                             atatAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCg                                                                                                                                                                                   <  2:140732/251‑1 (MQ=255)
                                                                  aTCACAACGGTGGGCTTGTTATTACGGGTGATAATGCCACGGTAAACAACAATGGCAAAACCACCGTTGACGGGAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGTGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCGCGGTATGGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACggc                                                                                                                                                                              <  2:147619/251‑2 (MQ=14)
                                                                  aTCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAACGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGGGGATAACAAGGGCACCATGACCGTCACCGACCCGGAGTCCATGGGACTCCAGATTGACGGt                                                                                                                                                                              >  2:159675/1‑251 (MQ=255)
                                                                    cacaACGGTGAGCTTGTTATTACGGGTGATAATGCCCCAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGagatacgg                                                                                                                                                                               >  1:311585/1‑241 (MQ=255)
                                                                    cacaACAGTGAGCTTGTTATTACGGGTGATAACGCCACAGTAAACAACAATGGCAAAACCACCGGTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAg                                                                                                                                                                                       <  2:311585/240‑1 (MQ=14)
                                                                         cggggAGCTTGGTATTACGGGAGATAATGCCACAGTAAACAACTATGGCAAAACCACCATTGACGGTAAAGATACCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATATGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACAGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCTATCCGTAGTCCATGGATATCCAGATCGACGGTGACAAgg                                                                                                                                                                       <  2:72640/247‑1 (MQ=18)
                                                                                               gATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCaa                                                                                                                                                             >  2:207633/1‑239 (MQ=14)
                                                                                               gATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCaa                                                                                                                                                             <  1:207633/239‑1 (MQ=14)
                                                                                                              aacaacAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATCGACGGTGACAAGGCCGTCGTCAATAACGAAGGCGAGGGCCCAATCACCa                                                                                                                                   >  1:322773/1‑250 (MQ=255)
                                                                                                                                              aGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtgatg                                                                                                  <  2:210958/251‑1 (MQ=255)
                                                                                                                                                     aCCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGCGGCACCGGCACCCAGAtagatgggaaggacgcccc                                                                                           >  1:138862/1‑233 (MQ=11)
                                                                                                                                                            cGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCCACGGTGGCACCGGCACCCAGAtcaacggtaactacgccaccgccaa                                                                                     >  2:270263/1‑226 (MQ=11)
                                                                                                                                                                aaTCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAATGATATCACCGGAGCCAGCGCGACGGGGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCGTCGTCAATAACGAAGGCGATAGCAACATCACCAACGGTGGCACCGGCACCCAGATTAACGGGGATGACGCCACGGCAaacaac                                                                                 >  2:362866/1‑250 (MQ=17)
                                                                                                                                                                   cAACAGTAATAGCGGGAAAGTGATTCATGACGTGGATCTGGATGACGGCAGCGGCGGTCACGGTAATGATATCACCGGTGACAGCGCGACGGTGGATAACCAGGGCACCATGAGCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCAGCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTAGCACCGGCACCCCGATTAATGCTGATGACGCCACGGCaaaaaccaca                                                                              <  2:4136/250‑8 (MQ=255)
                                                                                                                                                                     aCGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGc                                                                           <  1:341641/251‑1 (MQ=255)
                                                                                                                                                                            taaCGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGCCCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCAc                                                                                                                                           >  1:52263/1‑180 (MQ=255)
                                                                                                                                                                            taaCGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGGATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCggtagcac                                                                                                                                           <  2:52263/180‑8 (MQ=255)
                                                                                                                                                                                  gAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAAGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTg                                                              >  1:61815/1‑251 (MQ=255)
                                                                                                                                                                                      gTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACgg                                                          >  1:139324/1‑251 (MQ=255)
                                                                                                                                                                                       tGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCCCGGCAAACAACAACGGCAAAACCACCGTTGACGGc                                                         >  1:79947/1‑251 (MQ=255)
                                                                                                                                                                                                 cGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGTACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCa                                               <  2:186695/251‑1 (MQ=255)
                                                                                                                                                                                                   gCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAcc                                             <  1:366869/251‑1 (MQ=255)
                                                                                                                                                                                                                 aGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGGATCCAGATCGACGGTGACAAGCCCATCGTCAAAAACGAAGGCGAGAGCGCCATCACCATCGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCa                                >  1:108679/1‑250 (MQ=255)
                                                                                                                                                                                                                           gTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATTGACGGCGACCAGGCGATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAACGGTGACGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGATGGCAAGGATTCCACCGGTACGGAAATcgctggcaataac                     >  1:10964/1‑239 (MQ=11)
                                                                                                                                                                                                                           gTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCactggcaataac                     >  2:157091/1‑240 (MQ=255)
                                                                                                                                                                                                                           gTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCa                                <  2:315100/240‑1 (MQ=255)
                                                                                                                                                                                                                           gTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCa                                >  1:315100/1‑240 (MQ=255)
                                                                                                                                                                                                                               cGGTATTGATATCACCGGTGACAGGGCGAAGGTGGATAACAAGGACACCATGACCGACACCGATCCGGAGTCCATGGGTAACAAGACCGACGGTGACAAGGGCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCCCCCATATTAATGGTGATGACGCCACAGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaa                 <  2:353966/251‑1 (MQ=18)
                                                                                                                                                                                                                                ggTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaa                 >  1:106434/1‑250 (MQ=255)
                                                                                                                                                                                                                                     tGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGCGGCACCGGCACTCAGATTAACGAGGTTGAAACCACAGCAAACAACACCGGAGAAACCATCGATGACGGCAAAGATCGCCCCGGTACGGAAAACAACAGAAATAACGGAAAAGtt            >  1:231559/1‑250 (MQ=17)
                                                                                                                                                                                                                                            aCCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAGTCCACCGGTACGGAACTCAACGGTAATAACGGAAAAGTTATTCAgg     >  1:221418/1‑250 (MQ=255)
                                                                                                                                                                                                                                              cGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGGCCGGCACCGATCCGGAGGACATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCAACAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACg  <  1:103432/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       
ATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACG  >  NC_000913/1465960‑1466448

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