Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,472,791 A→G 17.6% pseudogene (2275/3497 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,472,7910AG17.6% 179.3 / 5.6 86pseudogene (2275/3497 nt)ydbApseudogene, autotransporter homolog, interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base A (24/44);  new base G (16/2);  total (40/46)
Fisher's exact test for biased strand distribution p-value = 4.72e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGGGCGACAGCATCAATGTTGAAGGTATGAACGGTAACAACGGTAATAGCAGCGAAGTGCATTATGGCGACTATACGTTGCCGGATGTGCCGAAGCCCAATACGGTTAGTGTAACGTCGGGAAGTGATGAGGCTGGTGGCAGCATGAACAACCTCAACGGCTATGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGATGCGG  >  NC_000913/1472555‑1473024
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                               
gCAGGGCGACAGCATCAATGTTGAAGGTATGAACGGTAACAACGGTAATAGCAGCGAAGTGCATTATGGCGACTATACGTTGCCGGATGTGCCGAAGCCCAATACGGTTAGTGTAACGTCGGGAAGTGATGAGGCTGGTGGCAGCATGAACAACCTCAACGGCTATGTCGTCGGAACCAACGTTAACGGCAGCGCCGGGAAGCGCAAGGTTGACAATGCCGGCATGAACGGCGGGggga                                                                                                                                                                                                                                               >  1:319305/1‑239 (MQ=255)
     gCGACAGCATCAATGTTGAAGGTATGAACGGTAACAACGGTAATAGCAGCGAAGTGCATTATGGCGACTATACGTTGCCGGATGTGCCGAAGCCCAATACGGTTAGTGTAACGTCGGGAAGTGATGAGGCTGGTGGCAGCATGAACAACCTCAACGGCTATGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGGAGCTGAAGGTTAACAATGCCAGCCTGAACGGCGTGGGGATTaa                                                                                                                                                                                                                                           >  1:179603/1‑238 (MQ=255)
             aTCATTGTTGAAGGGTTGAAGGGGAACAAGGGGAAAGGCAGGGAAGTGGATTATGGGGACTATACGTTGCCGGAGGGGCGGAAGCCCAATACGGTTAGGGTAACGGCGGGAAGTGATGAGGCTGGTGGCAGCATGAACAACCTCAACGGGTATGTCGTCGGTACCAACGGTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCtt                                                                                                                                                                                                                                 <  2:342934/239‑1 (MQ=255)
              tCAATGTTGAAGGTATGAACGGTAACAACGGTAATAGCAGCGAAGTGCATTATGGCGACTATACGTTGCCGGATGTGCCGAAGCCCAATACGGTTAGTGTAACGTCGGGAAGTGATGAGGCTGGTGGCAGCATGAACAACCTCAACGGCTATGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGGAGCTGAAGGTTAACAATGCCAGCATGAACCGCGGGGGGATTAACACGGGCttt                                                                                                                                                                                                                                >  1:114865/1‑239 (MQ=255)
                       aaGGTATGAACGGTAACAACGGTAATAGCAGCGAAGTGCATTATGGCGACTATACGTTGCCGGATGTGCCGAAGCCCAATACGGTTAGTGTAACGTCGGGAAGTGATGAGGCTGGTGGCAGCATGAACAACCTCAACGGCTATGTCGTCGGGACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCAAGCAGAAACGGGGGGGGGATTAAAACgggacttaacccgggg                                                                                                                                                                                                                        >  2:148285/1‑226 (MQ=255)
                          gTATGAACGTGAAGTAGGTTAATGGCAGCGAATGGGATTAGGGCGATGATACGTTTCCGGATGTGCCGAGGCCAAATACGGTTAGGTTAACGTCGGGAAGTGATGAGGCGGGTGGCAGCATGAACAACCTCAACGGCTATGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTa                                                                                                                                                                                                                      <  2:406561/237‑1 (MQ=255)
                            aTGAACGGTAACAACGGTAATAGCAGCGAAGTGCATTATGGCGACTATACGTTGCCGGATGTGCCGAAGCCCAATACGGTTAGTGTAACGTCGGGAAGTGATGAGGCTGGTGGCAGCATGAACAACCTCAACGGCTATGTCGTCGGGTCCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAAAAATGCCAGCATGAACGGCGTGGGGGTTAACACGGGCTTTTCACCTGGCACGGc                                                                                                                                                                                                                  >  1:591674/1‑239 (MQ=255)
                                          cGGTAATAGCAGCGAAGTGCATTATGGCGACTATACGTTGCCGGATGTGCCGAAGCCCAATACGGTTAGTGTAACGTCGGGAAGTGATGAGGCTGGTGGCAGAATGAACAACCTCAACGGCTATGTCGTCGGGACCAACGTTAACGGCAGCGCCGGGAAACTGAGGGTTAACAAGGCCAGCCGGAAGAGGGTGGGGGTTAACAAGGGCTTTAACGCGGGGAGGGccggccacaccggga                                                                                                                                                                                                    >  2:108856/1‑225 (MQ=255)
                                                tAGCAGCGAAGTGCATTATGGCGACTATACGTTGCCGGATGTGCCGAAGCCCAATACGGTTAGTGTAACGTCGGGAAGTGATGAGGCTGGTGGCAGCATGAACAACCTCAACGGCTATGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGGAATTAAAACGGGCTTTACCCCTGGCACGGCAGaaccccctggggggtttt                                                                                                                                                                                             >  1:848961/1‑222 (MQ=255)
                                                   cagcGAAGTGCATTATGGCGACTATACGTTGCCGGATGTGCCGAAGCCCAATACGGTTAGTGTAACGTCGGGAAGTGATGAGGCTGGTGGCAGCATGAACAACCTCAACGGCTATGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCCTGAACGGCGTGGAGATTAAAACGGGCTTTTCCGCTGGGACGGGAGAcaccaacttggggttttgta                                                                                                                                                                                          >  1:539360/1‑224 (MQ=255)
                                                    cgcGAAGTGCATTATGGCGACTATACGTTGCCGGATGTGCCGAAGCCCAATCCGGTTAGTGTAACGTCGGGAAGTGATGAGGCTGGTGGCAGCATGAACAACCTCAACGGCTATGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGAt                                                                                                                                                                                          <  1:709301/236‑1 (MQ=255)
                                                     gcGAAGTGCATTATGGCGACTATACGTTGCCGGATGTGCCGAAGCCCAATACGGTTAGTGTAACGTCGGGAAGTGATGAGGCTGGTGGCAGCATGAACAACCTCAACGGCTATGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCCCTGGTACGGCGGAACCCCCTTGGGGTTTTGTTAAc                                                                                                                                                                                       >  1:159921/1‑239 (MQ=255)
                                                      cGAAGTGCATTATGGCGACTATACGTTGCCGGATGTGCCGAAGCCCAATACGGTTAGTGTAACGTCGGGAAGTGATGAGGCTGGTGGCAGCATGAACAACCTCAACGGCTATGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGGGGGGATTAACACGGGCTTTACCGCGGGGACGGCAGACACCACTGGGGGTTTTGATAACg                                                                                                                                                                                      >  1:770065/1‑239 (MQ=255)
                                                         tgTGCTTTTTGGGGTCTATACGTTGCCGGGTGGGCCGAAGCCAAATACGGTTAGTGTAACGTCGGGAAGTGATGAGGCTGGTGGCAGCATGAACAACCTTAACGGCTATGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAg                                                                                                                                                                                  <  2:233612/238‑1 (MQ=255)
                                                         aGTGCATTATGGCGACTATACGTTGCCGGATGTGCCGAAGCCCAATACGGTTAGTGTAACGTCGGGAAGTGATGAGGCTGGTGGCAGCATGAACAACCTCAACGGCGATGTCGTCGGTACCAACGTTAACCGCAGCGCCGGGAAGCTGAAAGTTACACATCCCCACCTGAACGGGGGGGGGATTACACCCGGCTTTTCCCCTGGTCAGGCGAGAcaccccgcgaaagttttaaaaaagg                                                                                                                                                                                     >  2:486217/1‑218 (MQ=255)
                                                             tttttgggggtttttaCGTTGCCGGATTTGCCGAAGCCCAATACGGTTAGGGTAACGGCGGGAAGTGATGAGGCTGGTGGCAGCATGAACAACCTCAACGGCTATGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGa                                                                                                                                                                              <  2:782443/225‑1 (MQ=255)
                                                             gaTTTTGGCGACTTTGCGTTTCCGGGTGTGGCGAAGCCCAATGCGGTTAGTGTAACGGCGGGAAGTGATGGGGCTGGTGGCAGCATGAACAACCTCAACGGCTATGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGa                                                                                                                                                                              <  2:179603/238‑1 (MQ=255)
                                                                tATGGCGACTATACGTTGCCGGATGTGCCGAAGCCCAATACGGTTAGTGTAACGTCGGGAAGTGATGAGGCTGGTGGCAGCATGAACAACCTCAACGGCTATGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGc                                                                                                                                                                                                                           <  2:236371/194‑1 (MQ=255)
                                                                tATGGCGACTATACGTTGCCGGATGTGCCGAAGCCCAATACGGTTAGTGTAACGTCGGGAAGTGATGAGGCTGGTGGCAGCATGAACAACCTCAACGGCTATGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGc                                                                                                                                                                                                                           >  1:236371/1‑194 (MQ=255)
                                                                  tGGCGACTATACGTTGCCGGATGTGCCGAAGCCCAATACGGTTAGTGTAACGTCGGGAAGTGATGAGGCTGGTGGCAGCATGAACAACCTCAACGGCTATGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTa                                                                                                                                                                                                                      <  2:226420/197‑1 (MQ=255)
                                                                  tGGCGACTATACGTTGCCGGATGTGCCGAAGCCCAATACGGTTAGTGTAACGTCGGGAAGTGATGAGGCTGGTGGCAGCATGAACAACCTCAACGGCTATGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTa                                                                                                                                                                                                                      >  1:226420/1‑197 (MQ=255)
                                                                            actgtGCCGGATGTGCCTAAGCCCAATGCGGTTAGTGTAAGGTCGGGAAGTGAGGAGGCTGGTGGCAGGATGAACAACCTCAACGGCTATGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGAc                                                                                                                                                              <  1:401473/235‑1 (MQ=255)
                                                                                 gCCGGATGTGCCGAAGCCCAATACGGTTAGTGTAACGTCGGGAAGTGATGAGGCTGGTGGCAGCATGAACAACCTCAACGGCTATGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGAAGCTCAAGGTTCACAATTCCAGCATGAACGGCGTGGGAGTTAACACGGGCTTTACCGCTGGGACGGCCGAAACCACCGTGAGTTTTGATAAAGTAGGGGGAGGAAGAGACCCGAccacc                                                                                                                                                          >  2:219573/1‑235 (MQ=255)
                                                                                  ccGGATGTGCCGAAGCCCAATACGGTTAGTGTAACGTCGGGAAGTGATGAGGCTGGTGGCAGCATGAACAACCTCAACGGCTATGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGGGATTAACACGGGCTTTACCGCTGGTACGGCCGCCACCCCTGGGGGTTTTGTTAACGTAGTGGAAGGTAGCAAACTTACCGACGc                                                                                                                                                        >  1:76109/1‑239 (MQ=255)
                                                                                        gtgCCGAAGCCCAATACGGTTAGTGTAACGTCGGGAAGTGATGAGGCTGGTGGCAGCATGAACAACCTCAACGGCTATGTCGTCGGTACCAACGTTAACGGCAGCCCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGCTTAACACGGGGGTTTCCGCTGGTACGGCCGAGCCCACTGGGGGTTTTGTTAACGTAGAGGGAGGGAGAAACACGACCGACGCGGaac                                                                                                                                                   >  2:834785/1‑236 (MQ=255)
                                                                                           ccGAAGCCCAATACGGTTAGGGTAACGTCGGGAAGTGATGAGGCTGGTGGCAGCATGAACAACCTCAACGGCTATGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCa                                                                                                                                                <  1:131174/238‑1 (MQ=255)
                                                                                                    aaTTCGGTTAGTGTAACGGCGGGAAGTGATGAGGCTGGTGGCAGCATGAACAACCTCAACGGCTATGTCGTCGGTTCCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAAc                                                                                                                                      <  1:266004/239‑1 (MQ=255)
                                                                                                      tACGGTTAGTTTAACGGCGGGAAGTGATGGGGCTGGGGGCAGGATGAACAACCTCAACGGCTATGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGt                                                                                                                                    <  1:306268/239‑1 (MQ=255)
                                                                                                      tACGGTTAGTGTAACGTCGGGAAGTGATGAGGCTGGTGGCAGCATGAACAACCTCAACGGCTATGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGAAGCTGGAGGTTAACAATGCCAGCATGAACGGCGGGGGGATTAAAACGGGCTTTACCGCTGGGACGGCAGACACCACCGTGAGTTTTGGTAAAGGAGTGGAAGGGGGCAACCCGACCGACGCGAAAGCCCTCCCCCCAAca                                                                                                                                     >  2:552593/1‑237 (MQ=255)
                                                                                                       aCGGTTAGTGTAACGTCGGGAAGTGATGAGGCTGGTGGCAGCATGAACAACCTCAACGGCTATGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGa                                                                                                                                                                                           <  1:289789/185‑1 (MQ=255)
                                                                                                        cGGTTAGTGTAACGTCGGGAAGTGATGAGGCTGGTGGCAGCATGAACAACCTCAACGGCTATGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTcc                                                                                                                                  <  1:483469/239‑1 (MQ=255)
                                                                                                           ttAGTGTAACGTCGGGAAGTGATGAGGCTGGTGGCAGCATGAACAACCTCAACGGCTATGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGGTAACGTAATGGAAGGTGGCAACCTGACCGACCCTGACGCCATCACCTCAACGGCCGTg                                                                                                                              >  1:584959/1‑239 (MQ=255)
                                                                                                           ttAGTGTAACGTCGGGAAGTGATGAGGCTGGTGGCAGCATGAACAACCTCAACGGCTATGTCGTCGGTACCAACGTTAACGGCAGCGACGGGAAGCTGAAGGTTAACAATGCCAGCCTGAAAGGCCGGGGGGTGAACACGGGCTTTACCGCTGGGACGGCAAACACCACCGTGAGTTTTGTTAACATAGTGGAAAGGAGCAACCTGAGCCACGCTggacccaacaacccaaccgccggc                                                                                                                               >  2:30677/1‑216 (MQ=255)
                                                                                                            tAGTGTAACGTCGGGAAGTGATGAGGCTGGTGGCAGCATGAACAACCTCAACGGCTATGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCCGCATGAACGGCGTGGGGATTAACACGGGCGTTGCCGCTGGTACGGCAGAACCCCCTGTGAGTTTTGATAACGTAGCGGAAGGGAGCAAACCGACCGACCCTGAccacccctccccaaacgcccggg                                                                                                                              >  2:814131/1‑217 (MQ=255)
                                                                                                                 tAACGTCGGGAAGTGAGGGGGCTGGGTGCAGCATGAACAACCTCAACGGCTATGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGAAGCGGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATgg                                                                                                                        <  2:435997/239‑1 (MQ=255)
                                                                                                                       cGGGAAGTGATGAGGCTGGTGGCAGCATGAACAACCTCAACGGCTATGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGAAACCACTGTGAGTTTTGATAACGTAATGGAAGGTAGCAAACTGACCGACCCGGACGCCATAACCTCAACGTCCGGGGCATGGACGGcc                                                                                                                  >  1:435997/1‑239 (MQ=255)
                                                                                                                         ggAAGTGATGAGGCTGGTGGCAGCATGAACAACCTCAACGGCTATGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCa                                                                                                                                                >  1:739736/1‑208 (MQ=255)
                                                                                                                         ggAAGTGATGAGGCTGGTGGCAGCATGAACAACCTCAACGGCTATGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCa                                                                                                                                                <  2:739736/208‑1 (MQ=255)
                                                                                                                             gtggtggggtgggtggAAGCAGGAACAACCTAAAGGGTTATGGCGGCGGTACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACCAGTCCAGACTGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAgg                                                                                                             <  2:916893/227‑1 (MQ=255)
                                                                                                                             gtggtgggggttgggggagcaTTAACAACCTCAACGGCTAGGTCGTCGGTTCCAACGTTACCGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAgg                                                                                                             <  2:848961/221‑1 (MQ=255)
                                                                                                                              gtaggaGGCTGGGGGCAGTATGAACAACCTCAACGGCTATGTCGGCGGTACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGgca                                                                                                           <  1:450059/235‑1 (MQ=255)
                                                                                                                                atgaGGCTGGTGGCAGCATGAACAACCTCAACGGCTATGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGGGGTTAACACGGGCTTTACCGCTGGGACGGGAGAAACCACTGGGGGGTTTGATAACGGAGTGGAAGGGAGCAACATAGCCCACCCTGACGCCATCACCTCAAAGGCCGGGGGAAGGACGGCCAAACgcaga                                                                                                         >  2:841950/1‑238 (MQ=255)
                                                                                                                                atgaGGCTGGTGGCAGCATGAACAACCTCAACGGCTATGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAg                                                                                                                                                                                  <  2:9187/168‑1 (MQ=255)
                                                                                                                                atgaGGCTGGTGGCAGCATGAACAACCTCAACGGCTATGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAg                                                                                                                                                                                  >  1:9187/1‑168 (MQ=255)
                                                                                                                                atgaGGCTGGTGGCAGCATGAACAACCTCAACGGCTATGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGGTAACAATGCCAGCATGAACGGCGTGGGGATTAACCCGGGCGTTTACCCTGGGACGGCAGACACCCACGTGGGTTTTGGTAAAGATAGGGAAGGGAGCACACCAGCCGACCCTGCCGCCATCCCCTCAACGGCTGGGGGATAGACGCCCAAAGggaga                                                                                                         >  2:704931/1‑238 (MQ=255)
                                                                                                                                     gCTGGTGGCAGCATGAACAACCTCAACGGCTATGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGAACCGa                                                                                                    >  1:842091/1‑239 (MQ=255)
                                                                                                                                        ggtggCAGCATGAACAACCTCAACGGCTATGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGCAGCTGACGCCATCACCTCCACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATg                                                                                                  >  1:691577/1‑238 (MQ=255)
                                                                                                                                        ggtggCAGCATGAACAACCTCAACGGCTATGTCGTCGGTACCAACGTTAACGGCAACGCCGGGAAGCTGAAGGTTAAAAATGCCAGCCTGAACGGCGGGGGGATTAACACGGGCATTACCGCTGGTACGGCCGGAACCCATGTGGGGTTTGCTACGATAGTGGAAGGGAGACACCAGAGCGAACCAGCCGCCATCCCCCCCACGTCCGGGGGGAGGGACGCCAAAGgcaaggccacaa                                                                                                  >  2:173589/1‑229 (MQ=255)
                                                                                                                                         gtggCAGCATGAACAACCTCAACGGCTATGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTAACGCTGGGACGGCAGACACCACTGTGAGGTTTGATAACGTATTGGAAGGTAGCAACCTGACCGACGCTGAc                                                                                                                                                    >  1:552211/1‑188 (MQ=255)
                                                                                                                                         ggtgcgtgctgagccacCTCAACGGACTGTGCGTGGGTTCCAAAGGTAACGGCGGGGCCGGGAAGCGGAAGGTTAACAATGCCAGCAGGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACCCCACTGTGAGTTTGGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGc                                                                                                 <  2:147830/223‑1 (MQ=255)
                                                                                                                                           ggCAGTATGAACAACCTCAACGGCTAGGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGcc                                                                                                <  1:486217/238‑1 (MQ=255)
                                                                                                                                           gccagaaTGAAGAACCTCAACGGCTATGTCGACGGTACACGATTGAACGGAAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGGGTGAGGATTACCACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAg                                                                                              <  2:271915/237‑1 (MQ=255)
                                                                                                                                                   gAACAACCGCAACGGATAGTGCGTCGGTACCAACTGTAACGGGAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTCCGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACg                                                                                      <  2:539360/239‑1 (MQ=255)
                                                                                                                                                     acaccccaaaCGGCTTTGGCGGCGGTACCAAGGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGGGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGt                                                                                     <  2:159921/230‑1 (MQ=255)
                                                                                                                                                      caacCTCAACGGTTATGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTg                                                                                   <  1:219573/239‑1 (MQ=255)
                                                                                                                                                         ccgcAACGGTTATGGTGTCGGTACCAACGTTAACGGCAGGGGCGGGAAGCTTAAGGTTAACAAGGCCAGGAGGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACg                                                                                <  2:584959/236‑1 (MQ=255)
                                                                                                                                                          ccgcaccggttttgtgtgcgGTACCAACGTTAACGGGCGGGCCGGGAAGCTGAAGGTTAACAATGCCAGCAGGAACGGCGGGGAGTTTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACg                                                                                <  2:114865/221‑1 (MQ=255)
                                                                                                                                                                     aTGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGTTAACGTAGTGGAAGGTAGCAACCTGACCGCCGCTGACGCCATCACCCCAACGGCCGGGGGATGGACTGCCCAAAGCAGCACCCGAGCCACCGGGAACCGTTACGGCACCaaaggc                                                                     >  2:229358/1‑234 (MQ=255)
                                                                                                                                                                         ggtcgGTACCAACGTTAACGGCAGCGCAGGGAAGCGGAAGGTTGACCATGCCAGCATGAGCGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGGCACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCaaaa                                                                 <  1:376019/238‑1 (MQ=255)
                                                                                                                                                                         cgtcgGTACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCaaaa                                                                 <  1:834785/239‑1 (MQ=255)
                                                                                                                                                                             ggTACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGCGGAGCTTAACACGGGCTTTACCGCTGGGAGGGGAGACACCCCTGGGAGTTTTTGTAACAGATAGGAAGGGATCAACCCGACCGCACCAGACGCCCTAACCCACACGTCCGCGGCATGGACTGCCCAAGGCGGCACCGATGCCAGCGGGAACGTTGGCGACACCACGAGAcaaaaga                                                              >  2:188959/1‑236 (MQ=255)
                                                                                                                                                                                  cAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAAAGTCCGTGGTATGGACTGCCAAAGGCCGCACCGATGCCAGCGGTAACGTTGACGGCACCATGAGCAAAAAAGCCTac                                                         >  1:47939/1‑238 (MQ=255)
                                                                                                                                                                                    aCGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACAcc                                                      <  1:643453/239‑1 (MQ=255)
                                                                                                                                                                                        tAACGGAAGCGGCGGGAAGGTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGa                                                    <  1:294005/237‑1 (MQ=255)
                                                                                                                                                                                         aaCGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTACGTGGTATGGGCTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGGGCAAAAATGCCTACACCCAtgg                                                 >  1:817294/1‑238 (MQ=255)
                                                                                                                                                                                         aaCGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGGAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGCGGGATGGACTGCCAAAGACCGCCCCCAGGCACCCGGTAACCTTGCCGCTACCACGACCAAAAAAACCCACAcccaagt                                                 >  2:310572/1‑234 (MQ=255)
                                                                                                                                                                                         aaCGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGGGGTATGGACTGCCCAAGGCAGCACCGATGCCAACCGGAACGGTGCCGGCACCAGGCGCAAAAAACCCCAAACCGaatgg                                                >  2:646911/1‑235 (MQ=255)
                                                                                                                                                                                            ggCAGCGCCGGGAAGCTGAAGGTAAACAATGCCAGGAGGAACGGCGTGGAGATTAACAAGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGAAAAAATGCCTACACCGATGTGGc                                              <  1:779508/239‑1 (MQ=255)
                                                                                                                                                                                                 cgcCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCCATGCCAGCGGTAACGTTTACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACaa                                         >  1:46569/1‑238 (MQ=255)
                                                                                                                                                                                                 cgcCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAAACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCCAAGGCAGCACCGATGCCAGCGGTAACGTTTCCGTCACCATGAGCAAAAATGCCTACACCCATGGGGCAACAg                                         >  1:598286/1‑239 (MQ=255)
                                                                                                                                                                                                  gcCGGGAAGCTTAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGCCACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAg                                         <  1:188959/238‑1 (MQ=255)
                                                                                                                                                                                                     gggAAGATGATGGTTAACAGGGCCAGCAGTGAGGGGGTGGGGATTAACACGGGCTTTCCCGCTGGTTCGGCAGCCACCACTGGGGGTTGTGATAACGTAGTGGAATGTAGCAACCTGACCGACGCTGACGACATCCCCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCGCCGATGCCAGCGGTAACGTGGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGc                                     <  2:166862/239‑1 (MQ=255)
                                                                                                                                                                                                        aaGCTGAAGGTTAACAATGCCATCATGAACGGCCGGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCt                                   <  1:156696/238‑1 (MQ=255)
                                                                                                                                                                                                        aaGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCt                                   <  1:148285/238‑1 (MQ=255)
                                                                                                                                                                                                              aaGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCCAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCCATGTGGCAACAAAAGCCTCGGGGaa                            >  1:300984/1‑239 (MQ=255)
                                                                                                                                                                                                              aaGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGGGGCAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGAAGCCCGCGGTAACGTTGACGTCCCCATGAGCAAAAAAGCCCACACCCGCGTGGCAAAAGaacccccgggaga                            >  2:531892/1‑227 (MQ=255)
                                                                                                                                                                                                               agtttAAAAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGGGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAAt                           <  2:591674/236‑1 (MQ=255)
                                                                                                                                                                                                                 gttacacaTGCCAGCATGAACGGGGTGGATATTAACACGGGCTTTACCGCGGGTGCGGCAGACCCCACTGGGAGTTTTGATAACGTAGTGGAAGGGAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCt                                   <  2:916223/222‑1 (MQ=255)
                                                                                                                                                                                                                  ttAACAATGCCAGCATGAACGGCGTGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGAc                        <  1:646911/238‑1 (MQ=255)
                                                                                                                                                                                                                        atccccgCAAGAACGGGCTGGAGATTAACCCGGGTTTTACCGCTGGTACGGAAGACCCCACTGTGAGTTTTGATACCGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATcgcg                  <  2:842091/233‑1 (MQ=255)
                                                                                                                                                                                                                         tGCCAGCATGAACGGCGTGGAGATTAACCCGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGa                 <  1:353174/239‑1 (MQ=255)
                                                                                                                                                                                                                          gCCAGCATGAACGGGGGGGAAATTAACCAGGGCTTTACCGCTGGTACGGCAGAGACCCCTGTGAGTTTTGATAACGTAGTGGAATGTAGCAACCTGCACGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGaa                <  1:886778/239‑1 (MQ=255)
                                                                                                                                                                                                                           ccAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGAAGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGGGAATgc                         >  1:916223/1‑228 (MQ=255)
                                                                                                                                                                                                                                aTGAACGGCGTGGAGTTTAACGCGGGTTTTACCGCGGGTACGGCATACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTTACGCCATCCCCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACt          <  2:691577/239‑1 (MQ=255)
                                                                                                                                                                                                                                    aCGGCGTGGAGATTAACCAGGGCCTGACCGGTGGTCCGGCAGACGCCACTGTGAGTTTGGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGACATCACCTCAACGGCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGAt      <  1:462998/239‑1 (MQ=255)
                                                                                                                                                                                                                                    aCGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGAt      <  1:841950/239‑1 (MQ=255)
                                                                                                                                                                                                                                        cGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCATCGGTGAATGACATCGCGAAAGAAACGGATCCgg  >  1:167354/1‑239 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                               
GCAGGGCGACAGCATCAATGTTGAAGGTATGAACGGTAACAACGGTAATAGCAGCGAAGTGCATTATGGCGACTATACGTTGCCGGATGTGCCGAAGCCCAATACGGTTAGTGTAACGTCGGGAAGTGATGAGGCTGGTGGCAGCATGAACAACCTCAACGGCTATGTCGTCGGTACCAACGTTAACGGCAGCGCCGGGAAGCTGAAGGTTAACAATGCCAGCATGAACGGCGTGGAGATTAACACGGGCTTTACCGCTGGTACGGCAGACACCACTGTGAGTTTTGATAACGTAGTGGAAGGTAGCAACCTGACCGACGCTGACGCCATCACCTCAACGTCCGTGGTATGGACTGCCAAAGGCAGCACCGATGCCAGCGGTAACGTTGACGTCACCATGAGCAAAAATGCCTACACCGATGTGGCAACAGATGCCTCGGTGAATGACATCGCGAAAGCACTGGATGCGG  >  NC_000913/1472555‑1473024

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: