Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,466,210 G→A 6.8% pseudogene (819/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,2100GA 6.8% 206.2 / 5.0 89pseudogene (819/2513 nt)ydbApseudogene, autotransporter homolog, interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base G (34/48);  new base A (3/3);  total (38/51)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGG  >  NC_000913/1465973‑1466445
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                            
cacTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGa                                                                                                                                                                                                                                             <  2:515409/239‑1 (MQ=255)
  ttgtgtgccgcgagggggTGGGCGAAGATGGCGGATCGAAGCAGCAAAATAAACACAGCGGGGAGGAGCTAATGAAGGGGGATAATGCCACGGTAAACAACAATGGCAAAGCCACCGTTGACGGCAAAGATTCCACCGGTACGGAAATTACCGGCAATAGCGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACg                                                                                                                                                                                                                                           <  1:812576/225‑1 (MQ=2)
  gtGTTGCCCTGCAAGGGGTGAGCCCAGAGGGCCCAAATGATGGGCAGTATAATCACAACGGTTAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACTGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCTGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGGGAc                                                                                                                                                                                                                                            <  2:473092/237‑1 (MQ=18)
  cttttggctttGAGGGGGGGAGCGCAGATGGCACAACGAAGGGGGAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCGCCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATTACGGGAAAGGTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACg                                                                                                                                                                                                                                           <  1:482249/230‑1 (MQ=21)
     ttGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGCACGGAAATCAACGGTAATAACGGGAAAGTGATTAAGGACGGCGACCTGGATGTCGCCGGCGGCGGTCACGGGTTTGCTATCACCGGCGGCAGCGCCACgggg                                                                                                                                                                                                                                        >  1:828097/1‑237 (MQ=255)
        cccTTGAAGGGGTGGGCGGAGATGGCGCAACGAAGGGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGAt                                                                                                                                                                                                                                     <  1:58222/239‑1 (MQ=255)
        cccTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCACCGGCGGCGGGCACGGTATTGATATACCCGGGGCAAGCGCAACGggggct                                                                                                                                                                                                                                     >  2:812576/1‑237 (MQ=25)
         ccTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGAt                                                                                                                                                                                                                                     <  2:46969/238‑1 (MQ=255)
          cTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGGGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGGATTGATATCACCGGGGCAAGCGCGCCGGGGGGTac                                                                                                                                                                                                                                   >  2:347486/1‑238 (MQ=25)
          cTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGCAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTACCGGGATTGGTATTCACGGTGACAGCCCGACGGGGGGTaa                                                                                                                                                                                                                                   >  2:83990/1‑239 (MQ=25)
                     gggcgcAGATGGCGCAACGAATTGGCAATTTAATCACAACGGTGAGCTTTGTATTACGGGTGATAAGGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCa                                                                                                                                                                                                                        <  2:751350/237‑1 (MQ=255)
                      aGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGGGCCAGCGCGACGGGGGATAACAAGGGGACCAt                                                                                                                                                                                                                       >  2:974526/1‑239 (MQ=255)
                        cgcAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTACACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCAGGGGAATGGTAGTCCCGGTGGAAGGGCGGCGGGGGACAACAAGGGccccttt                                                                                                                                                                                                                      >  2:568175/1‑235 (MQ=21)
                          cAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGAAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCATCGGCGCCGGTCCCGGCATTGGTATCACCGACGACGGCGCGACGGCGGATAACAAGGGCACCATTAcc                                                                                                                                                                                                                   >  2:839780/1‑239 (MQ=21)
                              aggcaaaaCGAAGTGGTCATATAATCACAACGGTGAGATTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGt                                                                                                                                                                                                                 <  2:1184228/231‑1 (MQ=255)
                                      gaaaGTGGCAATTAAATCAAAACGGTGGGGTTGTTTTTACGGGTGAAAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCg                                                                                                                                                                                                       <  1:535456/237‑1 (MQ=255)
                                      cGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCg                                                                                                                                                                                                       <  2:297063/239‑1 (MQ=255)
                                      cGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGGCATAATGCCACAGTAAACAACAATGACAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGTCGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCg                                                                                                                                                                                                       <  2:356449/239‑1 (MQ=255)
                                      cGAAGTGGCAATATAATCAAAACGGTGAGCTTGTTATTAGGGGTGATAAGTCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATGCCACCGGGACGGAAATCAACGGTAATAACGGGAAAGTGAATCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCg                                                                                                                                                                                                       <  1:862512/239‑1 (MQ=255)
                                      cGAAGTGGCAATATAATCAAAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCg                                                                                                                                                                                                       <  2:335947/239‑1 (MQ=255)
                                      cGAAGGGGAAAGATAATAACAACGTGTATCTTGTTATTACGGGTGATAATGTCACAGTAAACAACAATGGCAAAACCACCGTGGACTGGAAAGGTTCCACCGGTACGGAAATCAACGGGAATAACGGGAAAGGTATTCAGGCCGGCGATCTGGAGGTCAGCTGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCg                                                                                                                                                                                                       <  1:581785/239‑1 (MQ=17)
                                       tagaggtgtaaaagggaaagaCGGGGAGTTTGTATTTACGGGTGATAATGCCAAAGTAAACAACAATGGCAAAACCACAGTTGACGGTAAAGAGTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCGTGACCGTCACCGATCCg                                                                                                                                                                                                       <  1:981645/218‑1 (MQ=17)
                                        aaGTTGCAATATATACTCACCGGTGGTCTTTTTATTACGGGTGATATGTCCACAGATACCACCGAGGGCAAAACCACCGTTGACGGGAAAGATTCCACCGGTCCGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGTTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCg                                                                                                                                                                                                       <  1:568175/237‑1 (MQ=14)
                                        aaGTGGCAATTTAAAACCAACGGGTAGCTTGTTATTACGGGTGAGAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGTTCCAACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCg                                                                                                                                                                                                       <  1:784771/237‑1 (MQ=17)
                                           ttgaaaaattaTCACAACTTTGACATTGTTATTACGGATGAAAATGTCAAAATAAACAACAATGTTAAAACCACCGTTGACGGTAAAGATGCCACCTGTACGGAAATCAACGGTAATAACGGGAAAGTGATGAAGGACGGCGATCTTGATGTCAGCGGCTGCAGTCACGGTATTGATATCGCCGGTGACAGCGCGACGGTGGATAACAAGGACACCATGACCGTCACCGATCCAGAGTc                                                                                                                                                                                                  <  1:136609/229‑1 (MQ=255)
                                            ggCAATATAATCACAACGGTGAGCTTGTTTATACGGTTGATAATGACCAAGTAAACAACAAGGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTc                                                                                                                                                                                                  <  1:966502/238‑1 (MQ=255)
                                            gcaacactcaTCACCGCGGTGGGCTTTTTATTACGGGGGATTACGCCCACGCAAACACAAATGGCAAAACCAACGTTGACGGCAAAGATCCCACCGGTACGGAAATTGTCGGCAATAACGGGAAAGTGATTCAGGACGGCGATCGGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTc                                                                                                                                                                                                  <  1:974526/229‑1 (MQ=11)
                                              cAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCAACGGTGGATAAAAAGTGTACCATGACCCCCAccccgccgcagcc                                                                                                                                                                                                  >  2:988190/1‑225 (MQ=18)
                                                     aTCACAACGGTGAGCTTTTTATTTCGGGTTAAAATGCCACAGTAAACAACAATGGCAAAACCACCTTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATc                                                                                                                                                                                        <  1:114332/239‑1 (MQ=255)
                                                       cacaACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGGAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGAGGGCACGGAATTGAAATTACCCGTGACAGCGCGACGGTGGGACACAAAGGCACCCTGACCGCAACCGAACCCGAGTCCAAggggaaaac                                                                                                                                                                                      >  1:707603/1‑233 (MQ=18)
                                                                 aGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGCCCCCGATCCGGAGGCCATGGGAATCCGGATCGCCGGt                                                                                                                                                                            >  2:1248343/1‑239 (MQ=255)
                                                                     gtttattgTGGGTGATAATGCCAAAGGAAACAAAAATGGAAAAACCCCCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACa                                                                                                                                                                        <  1:733635/237‑1 (MQ=255)
                                                                      gTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGAATCACGATCGACGGGGGCaa                                                                                                                                                                       >  2:1029055/1‑239 (MQ=255)
                                                                      gTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGGATTGATATCACCGGTGACAGGGCGCCGGTGGATAACAAGGGCACCGTGACCGTCACCGATCCGGAGTCCATgggcgtgctcttatacacatctc                                                                                                                                                                        >  2:981645/1‑218 (MQ=17)
                                                                      gTTAATACGGGTTATAATGCCAAAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGAcca                                                                                                                                                                       <  2:828097/239‑3 (MQ=11)
                                                                       ttattaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGAAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTCTTGATATCACCGGTGACAGCGAGACGGGGGCTAACAAGGGCACCATGACCGTCACCGAACCGGAGTACATGGGTATCCAGATCGGCGGTGCCAAg                                                                                                                                                                      >  2:703405/1‑239 (MQ=17)
                                                                           taCGGTTTAAAATGCCACAGTGAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCAGCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCa                                                                                                                                                                  <  2:707603/239‑1 (MQ=255)
                                                                                gTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGAATTGATATCCCCGGTCACAGACCGACGGCCGAGAACAAGGGACCGATGACCGTCACCGATACGGAGTCCTTGGTTAACCAGAACACCGGTGACAAGGCCATCGTc                                                                                                                                                             >  2:604211/1‑239 (MQ=17)
                                                                                 tGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGGATTGATATCACCGGTGACAGCGCGACGGGGGGTAAAAAGGGCACCATGACCGTCACCGATCCGGAGTCCATgg                                                                                                                                                                                             >  2:261192/1‑206 (MQ=255)
                                                                                 agaaaaggCCACCGTTAAAAACAAGCGCCAAACCCCCGTGTGCCGTAAAGATCCCTCCGGTAGGGGAAATAGAGGTAATAGAGGGGAAGTGATTCAAGACCGCGAGCGGGAGTTCAGCGGGGGCGGTCGAGGGATGGAGAACACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCAATGGGTATCCAGATCGACAGTGACAATGCCATCGTc                                                                                                                                                             <  1:703405/231‑1 (MQ=9)
                                                                                         ccACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGGTCGACGGAGACAAGGCCCTAGTCAATAAccac                                                                                                                                                    >  2:1178471/1‑236 (MQ=255)
                                                                                            aactacatagcaaTTGCCAAACACCCGTTGACGGTAAAGATCCCCGCGGGTCGGAAAGCGACGGTAATGACGGGAAAGTGATTCAGGAGGGCGAACTGGATGTCAGCGGCGGCGGTCACGGTATTGAAATCACCGGTGACAGCGCGGCGGTGGATAACAAGGGCACCATGACCGTCTCCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCCATAACGAAGGc                                                                                                                                                 <  1:357468/229‑1 (MQ=17)
                                                                                               tAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCaa                                                                                                                                                           <  2:723240/226‑1 (MQ=255)
                                                                                               tAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGGATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCACTAACGCAGGCgag                                                                                                                                              >  2:211488/1‑239 (MQ=255)
                                                                                                         tGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGAACCGGAGTCCATGGGTATCCAGATGCACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATc                                                                                                                                     >  2:515975/1‑238 (MQ=255)
                                                                                                         tGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGGATTGATATCACCGGTGACAGCGCGACGGGGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGGATCCAGATCCGCCGGGACAAAGCCCCCGCCAAAAAAAAAGGGgaggtcccttata                                                                                                                                    >  1:723240/1‑229 (MQ=14)
                                                                                                                 ccaccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCCATCTGGATGGCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGt                                                                                                                            <  2:640308/239‑1 (MQ=255)
                                                                                                                   tccGTTGGCGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTgg                                                                                                                          <  2:819622/238‑1 (MQ=255)
                                                                                                                       ttGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGAc                                                                                                                                                                               >  1:1248986/1‑182 (MQ=255)
                                                                                                                       ttGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGAc                                                                                                                                                                               <  2:1248986/182‑1 (MQ=255)
                                                                                                                            ggTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGCCCATCGGCAATAACGAAGGCGACAGCACCATCACCAAcacgggcccccgca                                                                                                                  >  2:196000/1‑225 (MQ=11)
                                                                                                                            ggTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCAGCGCCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCCTCACCGAACCGGAGTCCATGGGGATCCAGATCGACGGGGACAAGGCCATCGTCAAAAAAGAAGGCGAGAGCACAATCACCAACGAGggcaccgcca                                                                                                                  >  1:154390/1‑238 (MQ=17)
                                                                                                                              tAAAGTTTCCCACGGTACGGGAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGACGTCGGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcaccggcacc                                                                                                                <  1:347486/238‑1 (MQ=255)
                                                                                                                              tAAAGATTCCGCCTGTACGGAACTCAGCGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCAACGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcaccggcacc                                                                                                                <  2:429196/238‑1 (MQ=255)
                                                                                                                                     tCCCCCGGTAAGGAAATCAACGGGAAAAACGGGAAAGTGATTCAGGACGGGGATCGGGATGGCAGCGGAGGCGGTCACGGTATTGATAACACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTa                                                                                                         <  2:1248714/238‑1 (MQ=17)
                                                                                                                                           gtaggggAAATAACGGGAAAAACGGGAAAGGGTTTCAGGCCGGCGATCTGGATGGCGAGGGCGGGGGGCACGGTATTGATATCACCGGTGACAGCGCGACGGGGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtgat                                                                                                 <  1:717332/234‑1 (MQ=17)
                                                                                                                                                        aCGGTAATTACGGGAAAGTGTTTCAGGACGGGGAGCTGGATGTCAGCGGCGGCGGTCACGGGATTGATATCCCCGGTGACAGCGCGACGGTGGATAACAAGGGCCCCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCAGCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCa                                                                                     <  1:515975/239‑1 (MQ=255)
                                                                                                                                                            taataaCGGGAAAGTTATTCAGGCCGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCCCCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCgaga                                                                                                                                             <  2:1179482/179‑1 (MQ=255)
                                                                                                                                                                                cggtcGGGGATTTGGAAGTCAAGGGTGGGGGTCCAGGTTTTGATTTCACCGGAGGCAGGGCGACGGGGGATAACAAGGGCACCAGGACCGTCACCGATCCGGATTCCATGGGTAGCCAGAGCGACGGTGGCAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCAACAACGGTGGCACCGGGACCCAGATTAATGGTGATGAGGCCACGGAAAACAACAACGGCAAAACCACCGtt                                                            <  1:83990/235‑1 (MQ=18)
                                                                                                                                                                                cggACGGGCATCTGGAGGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGtt                                                            <  2:1127326/238‑1 (MQ=255)
                                                                                                                                                                                    cGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAggg                                                                                                                                                                                                                             <  1:958946/75‑1 (MQ=39)
                                                                                                                                                                                    cGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAggg                                                                                                                                                                                                                             >  2:958946/1‑75 (MQ=39)
                                                                                                                                                                                          tCTGGCTGTCCGCGGCGGCTGTCACGGTATTGATATCACCGGTGACCGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAg                                                   <  2:776731/238‑1 (MQ=255)
                                                                                                                                                                                                 ggcAGCGGCGGCGGTCACGGTATTGATATCCCCGGTGAAAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAc                                           <  2:592637/237‑1 (MQ=255)
                                                                                                                                                                                                 gTCAGCGGCGGCGGTCACGGTATTTATATCCCCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAc                                           <  2:107203/239‑1 (MQ=255)
                                                                                                                                                                                                    aGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTc                                                                                                                                                                                                  >  1:762234/1‑86 (MQ=38)
                                                                                                                                                                                                      cggcggcggTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGGACCATGCCCGTCACCGATCCGGCGTCCATGGGTATCCAGATCGACGGGGACAAAGCCCTCGTCAATAACGACGGCGAGACCAACATCCCCAACGGAGGCACCGGCCACCAGCTTAAAGGGGATGACCCACCGCCAACAAACACCGGCAAACCCACCCTTTACGGCAAGGATCCCACCggaa                                      >  1:987181/1‑237 (MQ=14)
                                                                                                                                                                                                           gcggTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTa                                                                                                                                                                                          <  2:750394/87‑1 (MQ=255)
                                                                                                                                                                                                           gcggTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAAGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTa                                                                                                                                                                                          >  1:750394/1‑87 (MQ=255)
                                                                                                                                                                                                               tCACGGTATTGATATCCCCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATc                              <  2:822353/238‑1 (MQ=255)
                                                                                                                                                                                                                 caggTATTGATATCACCGGTTGCAGCGCGGAGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAAc                           <  2:154390/237‑1 (MQ=255)
                                                                                                                                                                                                                     tATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCCCCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAGCGGt                                                                                                                            <  2:695853/139‑1 (MQ=255)
                                                                                                                                                                                                                       ttGAAATCCCCGGTGACAGGGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGAGCCGGATGCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGtaa                      <  2:987181/238‑1 (MQ=255)
                                                                                                                                                                                                                       ttGAAATCACCGGTGACAGCGAGACGGGGGATAACAAGGGGACCACGACCGGCACCGATCCGGAGGCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGtaat                     <  2:561752/239‑1 (MQ=255)
                                                                                                                                                                                                                            atCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGAACCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACgg                <  2:346059/239‑1 (MQ=255)
                                                                                                                                                                                                                             tCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCACCGGCAAACAACAACGGCAAAACCTCCGTTGACGGCAAGGATTCCACAGGTACGGCAAACaaagccactaccgga               >  2:1232769/1‑226 (MQ=17)
                                                                                                                                                                                                                             tCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACAGTGGCACCGGCACCCAGGTTAATGGTGATGACGCCACGGCAAACACAAACCGCAAAACCACCGTGTACGGCACGGACTCCACCGGGACGCaaacaaaagggaaaaaagga               >  1:346059/1‑222 (MQ=17)
                                                                                                                                                                                                                               aCCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAg                                                                                                                                                                                                    >  1:195164/1‑57 (MQ=255)
                                                                                                                                                                                                                               aCCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATc                                                                                                                                                                                        >  2:1289692/1‑69 (MQ=255)
                                                                                                                                                                                                                               aCCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaacaac                                                                              <  2:562571/174‑1 (MQ=255)
                                                                                                                                                                                                                               aCCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaacaac                                                                              >  1:562571/1‑174 (MQ=255)
                                                                                                                                                                                                                               aCCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCGAAACCCACGTTGACGGCAAGGAATCCACCGGAAAGGAAATAAACGAAAAAAAAGGaaa             >  1:662242/1‑239 (MQ=17)
                                                                                                                                                                                                                               aCCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCACAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGCAATAACGGaaa             >  2:1144488/1‑239 (MQ=255)
                                                                                                                                                                                                                               aCCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACAACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaa             >  2:468080/1‑239 (MQ=255)
                                                                                                                                                                                                                               aCCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCAACGGCAAACAAAAACGGCAAAACCACCGTTGCCGGCAAGGATTCAACCGGTACGGAAATAAACGGTAATACCGGgaa             >  2:1172515/1‑239 (MQ=255)
                                                                                                                                                                                                                               aCCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGGGATGACGCCACGGCAAACAACAACGGCAAAACCACCGGTGCCGTAAAGGAATCCCCCGGTACGGAAACCACCGCAAGAAACGGaaa             >  2:854519/1‑239 (MQ=17)
                                                                                                                                                                                                                               aCCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATAACCAACGGTGGCACCGGCACCCAGATTAATGGTAATGACGCCCAGGAAAACAACAACGGAAAAACCACCGTTAACGGCAAGGATTCCACCGGAACCAAAATCAAGGGGAATAACGGcaa             >  2:1262523/1‑239 (MQ=17)
                                                                                                                                                                                                                               aCCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGGGACAAGGCCATCGTCAACAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGCTTAATGGGGATGACGCCCCGGCAAAAAACAACAGAAAAACCACCGTTGCGGGAAAGAGTTCCCCCGGAAAGAAAAACAACGGGAAAAAAGGaga             >  1:427732/1‑237 (MQ=17)
                                                                                                                                                                                                                                      ccAGCGCGACGGGTGATAACAAGGGCACAATGACCGTCACTGGGCCGGAGTCCAAGGGGAACCAGATCGACGGTGATAAGGGCATCGTCAATAACGAAGGAGAGAGCACCATCACCAACGGTGGCACCGGCATCCAGATTAAGGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTAtt      <  2:139208/238‑1 (MQ=18)
                                                                                                                                                                                                                                          cgcgACGGTGGATAACAAGGGCACTATGACCGTCACCGATCCGGGGTCCATGTGCATCCAGAGCGACGGTGACAAGGCCATTGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGAACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAACCAACCGATGGAGGCCAGGCTGCCACCGGAACGGGAAAAAAAGGGAAAAACGGAAAAGTAAATCAga  >  1:924293/1‑238 (MQ=17)
                                                                                                                                                                                                                                          cgcgACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAg   >  2:529419/1‑238 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                            
CACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGG  >  NC_000913/1465973‑1466445

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