Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,466,579 G→C 10.8% pseudogene (1188/2513 nt) ydbA → pseudogene, autotransporter homolog, interrupted by IS2 and IS30

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,466,5790GC10.8% 176.4 / 6.8 70pseudogene (1188/2513 nt)ydbApseudogene, autotransporter homolog, interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base G (35/25);  new base C (6/3);  total (41/29)
Fisher's exact test for biased strand distribution p-value = 7.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGT  >  NC_000913/1466311‑1466841
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                      
ccATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGgttgtt                                                                                                                                                                                                                                                                  <  1:753231/275‑1 (MQ=255)
          cGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCACCGGTATCCAGGTTTACGGCGACCCGGCCCTTGTTAACAACCAAg                                                                                                                                                                                                                                                        >  1:378188/1‑275 (MQ=255)
           ggTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCCAGGATTCCACCGGCACGCAAATCAACGGTAATAACGGTCAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGGCACGGTATTGATATCACCGGCGACCGCGCCACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGTAGCCCATCGGTATCCGGGTTGACGGCGACCAGGCCAGCGTTAACAACGACgg                                                                                                                                                                                                                                                       >  1:343955/1‑275 (MQ=11)
           ggTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAg                                                                                                                                                                                                                                                        >  1:94034/1‑274 (MQ=255)
           ggTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGCCACCGATCCGGAGTTCATCGGTACCCCGGTAGCCGGCGCACAGGCGGTTGTGAACAACGaac                                                                                                                                                                                                                                                        >  2:95918/1‑273 (MQ=17)
               gcaccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCga                                                                                                                                                                                                                                                    >  2:396621/1‑274 (MQ=255)
               gcaccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCga                                                                                                                                                                                                                                                    <  2:94034/274‑1 (MQ=255)
                     gaacgccgcttAATCGTGATGGCGGCACGGCAAACAACAACGGAAAACCCACCGTTGACGGCAAGGGTTCCACCGGTACGGAAATTAACGGTAATAAAGGAAAAGTTATGCAGGACGTCGATCTGGATGTAATCGGCGGAGGGCAGGGTATAGATATCACCGGCGACAGAGAAACGGTGGATAACAAGGGAAACGTGACCGTCAGCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGcc                                                                                                                                                                                                                                             <  1:193510/266‑1 (MQ=255)
                             aTTAATGTTGATGACGCCACGGCAAACAACAACGGAAAAACCACCGTTGACGGCAAGGATTCAACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGCTGTCAGCGGCGGCGGTCACGGTATTGCTATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGAACGTCACCGAACCGGAGTCCATCGGTATACTGGTTGCCGGCGACCCGGCGCGTGGTAACAACGCAGGCGAGAGAGCCATCACCa                                                                                                                                                                                                                                      >  2:875989/1‑274 (MQ=17)
                                       atgaCGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCAACGGCGACAGCGCAACGGTGGATAACAAGGGAACCATGACCGTCACCGATCCGGCGCCCCTCGGTATCCAGGTTGACGGCGACCAGGCCATTGTTAACAACGAAGGCGAGAGCACCATCACCAACGGTgtctcc                                                                                                                                                                                                                           >  2:370586/1‑270 (MQ=14)
                                             ccACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAATGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACg                                                                                                                                                                                                                     <  2:726139/275‑1 (MQ=255)
                                             ccACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACg                                                                                                                                                                                                                     <  1:396620/275‑1 (MQ=255)
                                              cACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAAACACGAAGGCGAGAGCGCTATCACCAACGGCGGCCCCGGCACGc                                                                                                                                                                                                                    >  2:539713/1‑275 (MQ=255)
                                                             cGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAAAAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCa                                                                                                                                                                                                                   >  1:655037/1‑261 (MQ=255)
                                                             cGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCCGTTGTCAACAACGAAGGCGAGGGCGCCCGTCCCAACGGAGGAGCCGGAACCCCGAGTAAccgagacg                                                                                                                                                                                                     >  2:264551/1‑268 (MQ=14)
                                                                 aaaaCCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGAtt                                                                                                                                                                                                               <  2:655038/261‑1 (MQ=255)
                                                                   gaccaccGTTTAGCGTAAGGTTGCAACCGGAACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCAATCTGGAGGTCCGCGGCTGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGGGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCa                                                                                                                                                                                               <  2:758579/274‑1 (MQ=17)
                                                                   aaccgcctgtGACGGTAAAGATTCACCCGGTACGGAAAGAAACGGGAATAACGGAAAAGTTATGCAGGACGGCGATCTGGAGGGCAGCGGCGGCGGTCACGGTATTGATATCACCGGAGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCa                                                                                                                                                                                               <  1:197523/266‑1 (MQ=14)
                                                                   aaCCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGGTACCACCGGGGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCa                                                                                                                                                                                               <  2:600594/275‑1 (MQ=255)
                                                                                               ggTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACAAACGGTGGCACCGGCACGCAGATTAACGGTGATAACGCCACGGCAAACAACAACCGCAAAacccacg                                                                                                                                                                    >  2:553124/1‑270 (MQ=255)
                                                                                                gTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCGATCGTCAATAACGAAGGCGAGAGCACCATCACCAAGGGCGGCACCGGCCCCACGATCAACGGTAACGACGCCACCGCGAAAAAAAACGGCAAAACCACCgg                                                                                                                                                                   >  1:706224/1‑273 (MQ=11)
                                                                                                 tACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCCGTTGTCAACAACGAAGGCGAGAGCGCCATCACCAACGGCGGCACCGGCACGCCGATCAACGGTGACGACGCCACCGCAAACAACAACGGCAacacacacgtt                                                                                                                                                                  >  1:863601/1‑267 (MQ=14)
                                                                                                         acAACGGTAATAGCGGCATAGTTATGCAGACAGGCGACCTGGATGTCAGGGGCGGCGGTCAGTGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACCAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGTCACGCAGATTAACGGTGATTACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAg                                                                                                                                                         <  2:767902/274‑1 (MQ=17)
                                                                                                           aaCGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGCCCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGAAGGCGACCAGGCGGTCGTTAACAAAGAAGACGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGAGGACGCCACGGCAAACAACACCGGCAAAACCACCgcagccagcaacca                                                                                                                                                       >  2:596712/1‑262 (MQ=17)
                                                                                                                               gtttCAGGACGGCGCTCGTGGTGTCAGCGGCGGCGGGTACTGGATTGATATCACCTGTGAATGGGGAACGTTGGATAACAAGGGCACCATTACCGTCACCGTTCCGGAGTCCATCGGTATCCAGAGTGTCGGCTAGATGGCTGTTGATAATAACGAAGGCGCGGACACCATCACTAGCGGGGGCGTCGGCACGGTGAGTTGTGTAGATGAGGCAGCGGTAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTg                                                                                                                                   <  1:867642/273‑1 (MQ=2)
                                                                                                                               gagtCAGGACGGCGGTCTGGATGTCAGCGGCGGCGGTCACGGTATGGATATCGCCGGGGACAGCGCGACGGTTGATAACAAGGGCACCGTGACCGTCACCGATCCGGAGTCCATGGGTATCCAGGTCGACGGCGACCAGGCGATTGTTAAAAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTg                                                                                                                                   <  1:189397/272‑1 (MQ=11)
                                                                                                                                 ttCACGACGGCGATCTGGATGGCAGCGGCGGCGGACACGGTATTGACATCAGCGGCGACAGCGGAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCt                                                                                                                                 <  2:116016/275‑1 (MQ=255)
                                                                                                                                       acgtggtccTGCATGTAAGCGGCGGCGTTCACTGTATTGATACCCCCGGCGACAGCCCCACGGTGGATAACACGGGCACCATGACCGTAACCGATCCGGAGTCCATGGGTATCCAGTTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGGGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCGCCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAAt                                                                                                                           <  2:736404/267‑1 (MQ=17)
                                                                                                                                       aCGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGc                                                                                                                                                                              >  2:740859/1‑224 (MQ=255)
                                                                                                                                        cGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACg                                                                                                                                          >  2:260056/1‑259 (MQ=255)
                                                                                                                                         ggCGATCTGGATGTCAACGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTCGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGCGGCACCGGCACTCAGATACACGCTGTTGACCACACGTCACGCAACAACCGCAAGACCACGGTTGAAGCTAAGGCTTCGCCCGCTCCGGAAAATACTgcacatc                                                                                                                          >  1:797824/1‑268 (MQ=14)
                                                                                                                                           cGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCaaaa                                                                                                                                                                          <  1:740858/224‑1 (MQ=255)
                                                                                                                                           cGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTAAGGAAATGGCTGGCACTAACg                                                                                                                       >  2:439585/1‑275 (MQ=255)
                                                                                                                                            gATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGaaa                                                                                                                                      <  1:260056/259‑1 (MQ=255)
                                                                                                                                             agcTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACg                                                                                                                       <  2:567168/271‑1 (MQ=255)
                                                                                                                                                        aGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGagag                                                                                                                                                                                                                                                 >  1:672269/1‑140 (MQ=255)
                                                                                                                                                            ggggcggTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTGGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGGGATTCAGGa                                                                                                       <  2:152809/272‑1 (MQ=255)
                                                                                                                                                            gcggcggTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGcc                                                                                                                                                                                                                                             <  2:672270/140‑1 (MQ=255)
                                                                                                                                                                 ggtccgggtcagtaatacACCGCAGACAGCGCAACGGTAGATATCAAGGGCACGAAGACCTTCATCGGGCCGTAGTCCATCGGTATCCAGGTTGACGGCGACAAGGCGGTCGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCg                                                                                                  <  2:76653/257‑1 (MQ=17)
                                                                                                                                                                                 tCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCg                                                                                                                                                                                                                                                                        >  2:674776/1‑92 (MQ=34)
                                                                                                                                                                                   aCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcggcgg                                                                               >  1:333789/1‑275 (MQ=255)
                                                                                                                                                                                   aCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCGGGCAATAACGGGCAGGTAATTCAGGACGGCGATCTGGATGTCAgcggcgg                                                                               >  2:830324/1‑275 (MQ=255)
                                                                                                                                                                                   aCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGGCCCCGAACACGCATCTGTCGGTGTCGAGGTTGACGCCAGCCTGGAGCTTGGTAATCACGAAAGCGTGAGCGCCATAACCAAAGGTGGCACCGGAACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGGAGAGCCACAGTTACAGGCAAGAAGTCCAACCGGAGGGAAATGGGAGGCAATAACGGGAAGGAGATTTAGGCAGGCGaacggggttgcagcgtgg                                                                                >  1:171299/1‑257 (MQ=12)
                                                                                                                                                                                     cGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGCGGCACCGGCACGCAGATCAACGGTGACGACGCCACCGCAAATAACAGCGGAAAAACCACCGTTGACGGAAAAGATTCCACGGGTACGAAAATAGCGGGCAATAACGGGACTGTCAATCCGGACGGTAACCTGACTGTCAgaggcggcg                                                                             >  1:420112/1‑275 (MQ=11)
                                                                                                                                                                                           cAGGGCAACAGGGGGTAAAAAGGGCACCATGACCGTATCCGATCCGGAGTCCATGGGTATCCAGTTCGACGGCGACTAGGCCCTCGTGAAAAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAAGCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGTAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCa                                                                         <  1:87995/273‑1 (MQ=11)
                                                                                                                                                                                                    ggTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGata                                                              <  2:333790/275‑1 (MQ=255)
                                                                                                                                                                                                        gATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCa                                                           >  2:660726/1‑274 (MQ=255)
                                                                                                                                                                                                               agagccccgtGACCCTCCCCGCTCCTGATTCCATGGGTATCCAGGTTGGCGGTGACCAGGCGGTTGTCAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGAGTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGa                                                    <  1:842082/265‑1 (MQ=17)
                                                                                                                                                                                                                 agCACCCTGACCGTCACCGATCCGGAGGCCATCGGTAGCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAg                                                 <  1:660725/274‑1 (MQ=255)
                                                                                                                                                                                                                  gCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAgc                                                >  2:10558/1‑275 (MQ=255)
                                                                                                                                                                                                                      cgtGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGGTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGc                                                      <  1:10558/263‑1 (MQ=255)
                                                                                                                                                                                                                      cATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTCTTGATATCCCCGGCGACAGCGCaa                                            >  2:157606/1‑275 (MQ=255)
                                                                                                                                                                                                                        tGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGACGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACg                                          >  2:47774/1‑275 (MQ=255)
                                                                                                                                                                                                                        tGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGAc                                                                                                      >  2:619618/1‑215 (MQ=255)
                                                                                                                                                                                                                            cGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCg                                                                                                  <  1:619617/215‑1 (MQ=255)
                                                                                                                                                                                                                                aCCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGAAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTCATTCAGGACGGCGATCTGGGTGTCAGAGGCGGCGGGAACGGTATTGATATCCACGGCGACAGGGAAACGGCGAATa                                   >  2:759634/1‑274 (MQ=17)
                                                                                                                                                                                                                                aCCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGAAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTCCTGGCAATAACGGGAAGGTGATCCAGACCGGCGATCTGGCTGTCAGCGGCGGCGGTCACGGTCTGACTATCATCGGCGTCATCGCAACGGTGCAtc                                   >  2:443212/1‑273 (MQ=17)
                                                                                                                                                                                                                                    aTCCGGAGCCCATCGGTATCCAGGTTGACGGGGGCCAGGACGTCGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACCCAGATTAACGGGGATGAAGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAg                              <  2:700124/275‑1 (MQ=255)
                                                                                                                                                                                                                                    aTCCAGAGTCCATCGGTATCACGGTTGACGGCGACCCGGCGGTTGTTAACAACGAAGGCGAGAGCGCCACCAACAACGGGGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAg                              <  2:171299/275‑1 (MQ=255)
                                                                                                                                                                                                                                               aTCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCTATCACCAACGGCTGCACCGTCACGCAGATTGACGGTGACGACGCCACGGCTAACAAAAACGGCAAAACTCCCGGTGACGTCAAGGAATCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGGGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCAACGGTCACAGCGAAACCGGTGATAACAAgggacacatgca                   >  1:220990/1‑266 (MQ=17)
                                                                                                                                                                                                                                               aTCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGAAAATTGCTGGCAATAACGGGAAGGTGAATCAGGACGGTGACCTGGATGTCAGCGGCGGCGGTAACGGTATTGACACTGCCGGCGACCGACCCACGTTGGCTAAACAGGGcccccggac                   >  1:416305/1‑270 (MQ=17)
                                                                                                                                                                                                                                                tCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATCAACGTTGATGACGCCACGACCTACAACAACCGCCGAACCACCGTTGACGGAACGGGTTCCACCGGTACGGAAATTGATGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTTATATCACCGGCGAAAGCCAAACGGTGGATCACAAgcgacccctcacc                  >  1:356971/1‑263 (MQ=17)
                                                                                                                                                                                                                                                tCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGCGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACGGGTACGAAAATCGCTGGCAATAACGGCAAGGTCATTCAGGACGGTGATCTGGCTGTCAGCGGCGGTGGGCACGGTGTTGTAAGCACCGGCGACCACGCCACGGTCAACAACAAgggacacattgcg                  >  1:809751/1‑265 (MQ=17)
                                                                                                                                                                                                                                                    tATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGAAGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCa              >  2:592235/1‑275 (MQ=255)
                                                                                                                                                                                                                                                      tCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCACAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAGACTGCTGGCAATAACGGGAAGATGATTCAGGATGGCGATCTGGATGTCCGCGGCGCCGGGCCCGGTATTGATATCACCGGCGATAGCGCAACGGTGGAAAACAAGGGCACCATGACGGTCAc             >  1:49538/1‑274 (MQ=17)
                                                                                                                                                                                                                                                      tCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGCTCTGGATGTCAGCGGCGGCGGGCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCAc             >  2:596853/1‑274 (MQ=255)
                                                                                                                                                                                                                                                            ttGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCTGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACAATGACCGTCACCGCTcc       >  2:276969/1‑274 (MQ=255)
                                                                                                                                                                                                                                                            tgtgaggCGACCCGGCGGTTGTTAACAACGAGGGCGAGGGCGCCATCAACAACGGTGGCACCGGCACGCAGATGAACGGGGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGTTTCAGGACGGCGAGCTGGATGTCAGAGGCGGCGGAGACGGCACGGAGACAACGGGAGACAGGGAAACGGGGGAGAACAAGGGCACCATGACCGTCACCGATCCg      <  1:157606/270‑1 (MQ=17)
                                                                                                                                                                                                                                                              gACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGCCGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCAACCGGGACGGAAATTGCTGGCAATAACGGGAAGGGGATTCCGGACGGCGATCTGGCTGTCAGAGGCGGGGGGCAGGGGATTGAGAACACCGGCGACAGCGCAACGGTGGACACCAAGGGCACCATGTCAGTCCCCGATCCGGa    >  1:840726/1‑275 (MQ=17)
                                                                                                                                                                                                                                                                cGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGCCACGGTTTTGATATCCCCGGCGACCGCGCCACGGTGGACAACAAGGGCACCATGACCGTACCCGATCCGGAGt  >  2:680272/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                cGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTAACGACGCCACGGCCAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACGGGTACGGAAATTGCTGGCAATAACGTGATGGGGAATCAGGAAGGCAATCTGAATGTTACCGGCGGCGGGCACGGTGTTGCTATCACTGGCTACAGCTCAACGACGGATAACACGGCCACCGTGAcagctccacatgtgcagt  >  1:661627/1‑258 (MQ=17)
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                      
CCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGT  >  NC_000913/1466311‑1466841

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: