Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,210,041 G→A intergenic (+192/‑2) gltV → / → rrlE tRNA‑Glu/23S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,210,0410GA97.6% 142.2 / ‑4.3 84intergenic (+192/‑2)gltV/rrlEtRNA‑Glu/23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (36/45);  total (39/45)
Fisher's exact test for biased strand distribution p-value = 2.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCA  >  NC_000913/4209795‑4210268
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                   
aGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGt                                                                                                                                                                                                                                 <  1:490380/251‑1 (MQ=14)
aGGACACCGCCCTTTCACGGCGGTAACAGGGGTGCGAATCCCCTAAGGGACGCCCCTTGCTGGTTTGTGAGTCAAAGTCACCTACCTTCATCTCTCAAAACTCCTCTTCGGGTGATGTTTGAGATATTTGCTCTTTTAAAATCTGGATCAAGCTAAAAATTGAAACACTGAACAAACAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCACCCAGATGATGAATCGTAAGAAACATCTTCTGGTTGTAAGGt                                                                                                                                                                                                                                 >  1:86858/1‑251 (MQ=11)
 ggACACCGCCCTTTCATGGCGGTAACAGTGGTGCTAATCCCGTAGGGGACGCCACTTGCTGGTTTGTGATTGAAAGTCTCCTGCCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGtt                                                                                                                                                                                                                                <  2:597482/251‑1 (MQ=11)
                 cggcggTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAAACTCCTCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGCTCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTac                                                                                                                                                                                                                >  1:557405/1‑251 (MQ=14)
                  ggcggTAACAGGGGTTCGAAGCTCCTAGGGGACGCCACTTTCTGGTTTGTGAGTGAAAGTCCCCTGCCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTaca                                                                                                                                                                                                               <  2:557405/251‑1 (MQ=9)
                   gtggTAACAGGGGTTCGAATCCCTTTGGGGACGCCTCTTGGTTGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAAACTCAGCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAAGTTGCGCAGCACGATGATGAATCGTAAGAAACATCGTCGTGTTGTAAGGTTAAGCGACTAATCGTacac                                                                                                                                                                                                              <  2:55045/249‑1 (MQ=11)
                           aGGGGTTCGAATCCCCGAGGGGAAGCCACTTGCTGGTTTGTGAGTGAAAGTCAACTGCCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATg                                                                                                                                                                                                      <  2:192717/251‑1 (MQ=9)
                            ggggTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGc                                                                                                                                                                                                     >  1:225425/1‑251 (MQ=14)
                              ggttagactCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGACTCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGAAAGGTTAAGCGACTAAGCGTACACGGTGGATGccc                                                                                                                                                                                                   >  2:315258/9‑251 (MQ=255)
                                         ccTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCaga                                                                                                                                                                                        <  1:807680/251‑1 (MQ=14)
                                           tAGGGGACGCCACTTGCTGGTTTGTGAGTGAACGTCACCTTCCTTACTATCTCAAAACTCATCTTCGGGTGCTGTTTGAGATATTTGCTCTTTAACAATCTGGATCAAGCTGCAAATTGAAACACTGAACAACGACAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCACCACGATGATGAATCGTAAGAAACCTCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCTTGGCCGTCAGAgg                                                                                                                                                                                      >  2:449900/1‑251 (MQ=11)
                                             ggggACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAAACTCATCTTCGGGGGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGAACCAGGTGGATGCCCTGGccatccgggg                                                                                                                                                                                      >  2:610854/1‑240 (MQ=11)
                                                   gCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGCGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTCAGGTTAAGCGACTAAGCGTACACGGTGGAGGCCCTGGAAGGCCGAGGCGATgccg                                                                                                                                                                              >  2:234164/1‑248 (MQ=11)
                                                   gCCAATTGCTGGTTTGTGAGTGAAAGTCACCTGGCTTAATATCACAAAACTCTTCAGCGGGTGATGATTGAGATATTTTCTTATTAAAAATCTGTATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAAATTTCGCAACAGGATGATGAATCGTAAGAAACTGCGTCGGGTTATAAGGTCAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAg                                                                                                                                                                              <  2:235283/251‑1 (MQ=11)
                                                    ccACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCTTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCTTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAgg                                                                                                                                                                             <  1:100624/251‑1 (MQ=14)
                                                    ccACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAgg                                                                                                                                                                             >  1:235207/1‑251 (MQ=14)
                                                         tGCTGGTTTGTAAGTGAAAGTCACCTGCCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTTAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTTGCAGTCAGATGCGATTAAGGACGt                                                                                                                                                                         <  2:199996/250‑1 (MQ=255)
                                                                 tgtgAGTGAAAGTCCCCTGCCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCt                                                                                                                                                                 <  2:698326/250‑1 (MQ=14)
                                                                   tgagtgaAAGTCACCTGCCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCg                                                                                                                                                              >  1:645716/1‑251 (MQ=14)
                                                                    gagtgaAAGTCACCTGCCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCg                                                                                                                                                              <  1:429058/250‑1 (MQ=14)
                                                                        aaaaGTCACCTGCCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGTAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTACGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATa                                                                                                                                                           >  2:553254/2‑249 (MQ=9)
                                                                          aaGTCACCTGCCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTATCAAATTTTCGCAACACGATGATGAATAGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGCAGGACGTGCTAATCTGCGCTAAGCg                                                                                                                                                       >  1:219940/1‑251 (MQ=9)
                                                                                ccTGCCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGAGGATGCCCTGGCAGTCAGCGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGt                                                                                                                                                  >  2:51432/1‑250 (MQ=14)
                                                                                   gactaaaTATCTCAAAACTCATATTTGGGTGATGTTGGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACCCTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAgg                                                                                                                                              <  1:288/246‑1 (MQ=255)
                                                                                        aaTATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTGTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCTACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGata                                                                                                                                         <  1:520596/251‑1 (MQ=14)
                                                                                         atatCTCAAACCTCATCTTCGTGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGACCTGGCAGTCAGAGGCGATGAAGGACGTTCTAATCTGCGATAAGCGTCGGTAAGGTGatat                                                                                                                                        <  1:635055/251‑1 (MQ=9)
                                                                                         atatCTCAAAACTCATTTTCGGGTGATGTTTGAGATATTTTCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCGCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATTCCCGGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGatat                                                                                                                                        <  2:799120/251‑1 (MQ=255)
                                                                                         acatCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGatat                                                                                                                                        <  2:21236/249‑1 (MQ=14)
                                                                                            cctcAAAACTCAGCGTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAAGCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGaa                                                                                                                                     <  2:256712/250‑1 (MQ=9)
                                                                                               cAAAACGCATCGTCGGGTGATGTTTGAGATATTTGCTCTGTAAAAATCTGGTGCAAGCTGAAAAGTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAATGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAAc                                                                                                                                    <  2:32611/249‑1 (MQ=11)
                                                                                                 aaaCTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGtt                                                                                                                                >  1:119800/1‑251 (MQ=14)
                                                                                                   aCTCAACTTCGGGTATTGTTTGAGATATTTGCTCTTTAAAACTATGGATCAAGCTGAAAATTGAAACACTGAACAACGAATGTTGTGAGTGAGATC‑CTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCGAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTtat                                                                                                                              >  2:793706/1‑251 (MQ=11)
                                                                                                      cATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTtata                                                                                                                             <  2:755616/249‑1 (MQ=14)
                                                                                                      cATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAAGCTGGATCAAGCTGAAAATTGAAACCCTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTtata                                                                                                                             <  1:7156/249‑1 (MQ=14)
                                                                                                       aTCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAAcc                                                                                                                          >  1:731578/1‑251 (MQ=14)
                                                                                                        acgtcggtggATGTTTGAGATATTGGCTCTTGAAAAATCTGGAGCAAGATGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCg                                                                                                                         <  2:219940/242‑1 (MQ=11)
                                                                                                         cTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGGTATGAACCGTTATAACCgg                                                                                                                        >  2:418808/1‑251 (MQ=14)
                                                                                                               gTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGAttt                                                                                                                  <  2:732141/251‑1 (MQ=14)
                                                                                                                 gATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTcc                                                                                                                <  1:171513/251‑1 (MQ=14)
                                                                                                                 gATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTc                                                                                                                 <  1:353959/250‑1 (MQ=14)
                                                                                                                 gATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTc                                                                                                                 <  2:502436/250‑1 (MQ=14)
                                                                                                                  aTGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGACAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGGGATATGAACCGTTATAACCGACGATTTCCg                                                                                                               >  1:677334/1‑251 (MQ=9)
                                                                                                                       tGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACAGTTATAACCGGCGATTTCCGAATgg                                                                                                          >  1:338988/1‑251 (MQ=14)
                                                                                                                              tttGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGACTGGGGAAAcc                                                                                                   >  2:388937/1‑251 (MQ=14)
                                                                                                                                        taacaTCTGAATCAAGCTGAAAATTGAAACACTGACCAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGACTCGTAAGCAACATCTTCGAGATGTAAGGTTAAGCGACTCAGCGTACACGGTGTATGCCCTGGCAGTCAGAGGCGATGAAGAACGTGCTAATCTGCGGTAAGCGTCGGTCAGGTGATATGCACCGTTATAACCGGCGATTTCCGAATgggaggtcggcagagagtc                                                                                         >  2:788961/5‑235 (MQ=11)
                                                                                                                                        aaaaaTCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAGGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGttt                                                                                         >  1:362742/1‑251 (MQ=14)
                                                                                                                                        aaaaaTCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGtt                                                                                          <  1:570988/250‑1 (MQ=14)
                                                                                                                                        aaaaaTCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGGTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCTGGTTGTAAGGTTAAGCGACTAAGCGTACTCGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATggg                                                                                                         <  1:788961/235‑1 (MQ=9)
                                                                                                                                            aTCTGGATCAAGCTGAAAATTGAAACACTAAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCACCACGATGATGAATCGTAAGAAACATCTTCGGATTGTAAGGTTAAGCGACTAAGCGTACACGGTGGCTGCCCTGTCAGTAAGAGGCGGGGAAGGGCGAGCTAATCTGCGATAAGCGTCGGTAAGGTGCTATGAACCGTTTAAACCGGCGATTTACGAATGGGGAAAACCAGTGTGTTTCGa                                                                                      >  2:300021/1‑250 (MQ=11)
                                                                                                                                                     aagttGAAAATTGAAACACTGAACAATGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCATTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCa                                                                            <  1:520416/247‑1 (MQ=9)
                                                                                                                                                      aGCTGAAAATTGAAACTCTGAACCACTAAAGTTGTTCGTGAG‑TCTCTCAAATATTCGCAACACGATGATGAATCGTAAGAAACATCTTCTGGTTGTAAGGTTAAGCGACTAAGCGTACAAGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATACGCGCCGGTAAGGTGATATGTATCGTTATAACCGTCGATTTCCGAATGGGGAAACCCAGTGTGGTGCGACACACTATCAt                                                                           <  2:99728/251‑1 (MQ=11)
                                                                                                                                                      aGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCAt                                                                           >  1:640199/1‑251 (MQ=14)
                                                                                                                                                          gAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAAc                                                                       <  1:618096/251‑1 (MQ=14)
                                                                                                                                                            aaaTTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑GCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAGCATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTg                                                                     <  2:781946/251‑1 (MQ=14)
                                                                                                                                                               ttGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTgcat                                                                  >  1:349098/1‑248 (MQ=14)
                                                                                                                                                                 gAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATcc                                                                <  1:449217/251‑1 (MQ=14)
                                                                                                                                                                       cTGAACAACGAAAGTTGTTCGTGAG‑TCTATCAAATTTTCGCAACAAGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAAACAAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGt                                                          >  1:296919/1‑251 (MQ=9)
                                                                                                                                                                         gAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTa                                                        <  1:132055/251‑1 (MQ=14)
                                                                                                                                                                          aacaacGAAAGTTGTTCGTTAG‑TCTTTCAAATTTTCGAAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTaa                                                       <  2:726561/251‑1 (MQ=9)
                                                                                                                                                                          aacaacGAAAGGTGTTCGTGAG‑GTTCTCAAATTGTCGCAACACGTTGATGAACCGTAAGAAACATCTGCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGTCTATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGGTATAACCTGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCACATGTTaa                                                       <  2:327019/251‑1 (MQ=11)
                                                                                                                                                                             aacGAAAATTGTTCGTGAG‑TCTCTCAAATTTTCGCCACCCGAAGATACATAGTAAGACACATCTTAGGGTTGTAAGGTTAAGCGACTCAGCGTACACGAAGGATGCACTGGCAGTCAGAGGAGATGCAGGACATGCTAATCTACGATAAGCGTAGGTAAAGTGATATGAACCGTTATAACCGGCGATTTAAGACAAGAGAAACCCAGGGAGTTTTGACACACTAACATTAACTGAATCCATCGGTTaaat                                                     >  2:808709/1‑248 (MQ=2)
                                                                                                                                                                              acGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAg                                                   <  1:449900/251‑1 (MQ=14)
                                                                                                                                                                                gAAAGTTGTTCGTGAG‑TCTCGCAAATTTTCGCACCACGATGATGAATCGTAAGAAACATCTTCAGGTTGTAAGGTTAAGCGACTAAGAGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGAAGTGCTAATATGCGATAAGCATCAGTAGGGGGATATGAACCGTTATAACAGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACCATATCATTCCTTGAAACCATAGGTGAATgtgg                                                  >  2:432512/1‑247 (MQ=11)
                                                                                                                                                                                 aaaGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCg                                                >  1:369749/1‑251 (MQ=14)
                                                                                                                                                                                   aGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGTTTGTAAGGTTAAGCGACTAAGCGTACACGTTGGATGCCCTGGCAGTCAGATGCGATGa                                                                                                                                                                                >  1:105671/1‑121 (MQ=11)
                                                                                                                                                                                    gttgttTGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGCAAGATCTTCGGGTTGTAAGGTTCAGCGACTAAGCGTACACGGCGGATGCCCTGGCAGTCATAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGAGTTCCGAATGGTGAAACCCAGGGTGTTGCGACACACTATCATTACTTGAATCCATAGGTTAATGAGGCGaaa                                             <  2:86858/251‑2 (MQ=11)
                                                                                                                                                                                      tgttCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGTATCGTAAGCAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTAAGAGACGATGAAGGACGTGCTATTCTGCGATACGCGTCGGTAAGGAGATATGAACCGTTATAACCGGCGTTTTCCGACTAGGGAAACCCAGAGTGTTTTGACACACTATCATTAAATGAATCCATAGGTGAATGTGGcacaccg                                           >  2:204180/1‑245 (MQ=11)
                                                                                                                                                                                       gttCGTGAG‑TCTCTCAACTTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGGTAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGCTTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACCATCATTAACTGAATCCATAGGTTCATGAGGCggcccgg                                          >  2:421344/1‑245 (MQ=11)
                                                                                                                                                                                       gttCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCgg                                          >  1:580815/1‑251 (MQ=14)
                                                                                                                                                                                             gAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTACGGGTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAAAGGGGAAACCCAGTGTGTTTCGACACCCTATCATTAACTGAATCCATAGGTTAATAAGGCGAACCAGGGGAAc                                    >  2:792123/1‑251 (MQ=11)
                                                                                                                                                                                              aG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACt                                   <  2:743159/251‑1 (MQ=14)
                                                                                                                                                                                                ttctctcAAATTTTCGCATAACGATGATGATTCGTAAGAAACATCTTCGGGTTGTACGTTTAAGTGACTAAGCTTACACGGTCGATGCCCTGTCATTCACAGGCGATGAAGGACGTGCTACTCTGCGATAATCGTCGGTAAGGTTATATTAACCGTTATAACCGGCGATTTCCGAATGTGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACt                                   <  2:381823/249‑1 (MQ=11)
                                                                                                                                                                                                   tctcAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTTCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGt                                                          >  1:401359/1‑224 (MQ=14)
                                                                                                                                                                                                   tctcAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTTCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGt                                                          <  2:401359/224‑1 (MQ=14)
                                                                                                                                                                                                      cAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAAc                              <  1:782715/249‑1 (MQ=14)
                                                                                                                                                                                                      aaaaTTTTCGCAAAACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAGGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACAt                            <  2:362742/250‑1 (MQ=9)
                                                                                                                                                                                                       aaaTTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCACTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACAGTTATAACCGGCGCTTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCAGGGGAACTGAAACATc                           >  1:819524/1‑251 (MQ=9)
                                                                                                                                                                                                         aTTTTCGCAAAACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTa                         <  2:629239/251‑1 (MQ=14)
                                                                                                                                                                                                           tttCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGTTAGGCGTCGGTAAGGGGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAATTGAATCCATAGGTTAATGGGGCGAACCGGGGGAACTGAAACATcaa                         >  1:441918/1‑247 (MQ=11)
                                                                                                                                                                                                           ttgtgCAACACGATGATGAATCGTAAGAAACAGCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTaa                        <  2:225425/246‑1 (MQ=14)
                                                                                                                                                                                                             tCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGAATAAGCGTACACGGCGGATGCCCTGGCAGTCAGAGGCGATGAAGAACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACACAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACTGGGGGAACTGAAACATCTAAga                      >  1:404531/1‑249 (MQ=255)
                                                                                                                                                                                                                      gatgatGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAgg             >  1:733665/1‑250 (MQ=14)
                                                                                                                                                                                                                            gaATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGaa      <  1:553254/251‑1 (MQ=14)
                                                                                                                                                                                                                                 gTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAAAATCTAAGTACCCCGAGGAAAAGAAAtaa  >  1:37673/1‑248 (MQ=11)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                   
AGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAG‑TCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCA  >  NC_000913/4209795‑4210268

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