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

AGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCA  >  NC_000913/4209795‑4210268
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                   
aGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGt                                                                                                                                                                                                                                 <  1:490380/251‑1 (MQ=14)
aGGACACCGCCCTTTCACGGCGGTAACAGGGGTGCGAATCCCCTAAGGGACGCCCCTTGCTGGTTTGTGAGTCAAAGTCACCTACCTTCATCTCTCAAAACTCCTCTTCGGGTGATGTTTGAGATATTTGCTCTTTTAAAATCTGGATCAAGCTAAAAATTGAAACACTGAACAAACAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCACCCAGATGATGAATCGTAAGAAACATCTTCTGGTTGTAAGGt                                                                                                                                                                                                                                 >  1:86858/1‑251 (MQ=11)
 ggACACCGCCCTTTCATGGCGGTAACAGTGGTGCTAATCCCGTAGGGGACGCCACTTGCTGGTTTGTGATTGAAAGTCTCCTGCCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGtt                                                                                                                                                                                                                                <  2:597482/251‑1 (MQ=11)
                 cggcggTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAAACTCCTCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGCTCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTac                                                                                                                                                                                                                >  1:557405/1‑251 (MQ=14)
                  ggcggTAACAGGGGTTCGAAGCTCCTAGGGGACGCCACTTTCTGGTTTGTGAGTGAAAGTCCCCTGCCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTaca                                                                                                                                                                                                               <  2:557405/251‑1 (MQ=9)
                   gtggTAACAGGGGTTCGAATCCCTTTGGGGACGCCTCTTGGTTGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAAACTCAGCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAAGTTGCGCAGCACGATGATGAATCGTAAGAAACATCGTCGTGTTGTAAGGTTAAGCGACTAATCGTacac                                                                                                                                                                                                              <  2:55045/249‑1 (MQ=11)
                           aGGGGTTCGAATCCCCGAGGGGAAGCCACTTGCTGGTTTGTGAGTGAAAGTCAACTGCCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATg                                                                                                                                                                                                      <  2:192717/251‑1 (MQ=9)
                            ggggTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGc                                                                                                                                                                                                     >  1:225425/1‑251 (MQ=14)
                              ggttagactCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGACTCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGAAAGGTTAAGCGACTAAGCGTACACGGTGGATGccc                                                                                                                                                                                                   >  2:315258/9‑251 (MQ=255)
                                         ccTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCaga                                                                                                                                                                                        <  1:807680/251‑1 (MQ=14)
                                           tAGGGGACGCCACTTGCTGGTTTGTGAGTGAACGTCACCTTCCTTACTATCTCAAAACTCATCTTCGGGTGCTGTTTGAGATATTTGCTCTTTAACAATCTGGATCAAGCTGCAAATTGAAACACTGAACAACGACAGTTGTTCGTGAGTCTCTCAAATTTTCGCACCACGATGATGAATCGTAAGAAACCTCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCTTGGCCGTCAGAgg                                                                                                                                                                                      >  2:449900/1‑251 (MQ=11)
                                             ggggACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAAACTCATCTTCGGGGGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGAACCAGGTGGATGCCCTGGccatccgggg                                                                                                                                                                                      >  2:610854/1‑240 (MQ=11)
                                                   gCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGCGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTCAGGTTAAGCGACTAAGCGTACACGGTGGAGGCCCTGGAAGGCCGAGGCGATgccg                                                                                                                                                                              >  2:234164/1‑248 (MQ=11)
                                                   gCCAATTGCTGGTTTGTGAGTGAAAGTCACCTGGCTTAATATCACAAAACTCTTCAGCGGGTGATGATTGAGATATTTTCTTATTAAAAATCTGTATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAAATTTCGCAACAGGATGATGAATCGTAAGAAACTGCGTCGGGTTATAAGGTCAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAg                                                                                                                                                                              <  2:235283/251‑1 (MQ=11)
                                                    ccACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCTTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCTTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAgg                                                                                                                                                                             <  1:100624/251‑1 (MQ=14)
                                                    ccACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAgg                                                                                                                                                                             >  1:235207/1‑251 (MQ=14)
                                                         tGCTGGTTTGTAAGTGAAAGTCACCTGCCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTTAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTTGCAGTCAGATGCGATTAAGGACGt                                                                                                                                                                         <  2:199996/250‑1 (MQ=255)
                                                                 tgtgAGTGAAAGTCCCCTGCCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCt                                                                                                                                                                 <  2:698326/250‑1 (MQ=14)
                                                                   tgagtgaAAGTCACCTGCCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCg                                                                                                                                                              >  1:645716/1‑251 (MQ=14)
                                                                    gagtgaAAGTCACCTGCCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCg                                                                                                                                                              <  1:429058/250‑1 (MQ=14)
                                                                        aaaaGTCACCTGCCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGTAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTACGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATa                                                                                                                                                           >  2:553254/2‑249 (MQ=9)
                                                                          aaGTCACCTGCCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTATCAAATTTTCGCAACACGATGATGAATAGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGCAGGACGTGCTAATCTGCGCTAAGCg                                                                                                                                                       >  1:219940/1‑251 (MQ=9)
                                                                                ccTGCCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGAGGATGCCCTGGCAGTCAGCGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGt                                                                                                                                                  >  2:51432/1‑250 (MQ=14)
                                                                                   gactaaaTATCTCAAAACTCATATTTGGGTGATGTTGGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACCCTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAgg                                                                                                                                              <  1:288/246‑1 (MQ=255)
                                                                                        aaTATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTGTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCTACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGata                                                                                                                                         <  1:520596/251‑1 (MQ=14)
                                                                                         atatCTCAAACCTCATCTTCGTGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGACCTGGCAGTCAGAGGCGATGAAGGACGTTCTAATCTGCGATAAGCGTCGGTAAGGTGatat                                                                                                                                        <  1:635055/251‑1 (MQ=9)
                                                                                         atatCTCAAAACTCATTTTCGGGTGATGTTTGAGATATTTTCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCGCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATTCCCGGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGatat                                                                                                                                        <  2:799120/251‑1 (MQ=255)
                                                                                         acatCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGatat                                                                                                                                        <  2:21236/249‑1 (MQ=14)
                                                                                            cctcAAAACTCAGCGTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAAGCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGaa                                                                                                                                     <  2:256712/250‑1 (MQ=9)
                                                                                               cAAAACGCATCGTCGGGTGATGTTTGAGATATTTGCTCTGTAAAAATCTGGTGCAAGCTGAAAAGTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAATGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAAc                                                                                                                                    <  2:32611/249‑1 (MQ=11)
                                                                                                 aaaCTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGtt                                                                                                                                >  1:119800/1‑251 (MQ=14)
                                                                                                   aCTCAACTTCGGGTATTGTTTGAGATATTTGCTCTTTAAAACTATGGATCAAGCTGAAAATTGAAACACTGAACAACGAATGTTGTGAGTGAGATCCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCGAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTtat                                                                                                                              >  2:793706/1‑251 (MQ=11)
                                                                                                      cATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTtata                                                                                                                             <  2:755616/249‑1 (MQ=14)
                                                                                                      cATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAAGCTGGATCAAGCTGAAAATTGAAACCCTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTtata                                                                                                                             <  1:7156/249‑1 (MQ=14)
                                                                                                       aTCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAAcc                                                                                                                          >  1:731578/1‑251 (MQ=14)
                                                                                                        acgtcggtggATGTTTGAGATATTGGCTCTTGAAAAATCTGGAGCAAGATGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCg                                                                                                                         <  2:219940/242‑1 (MQ=11)
                                                                                                         cTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGGTATGAACCGTTATAACCgg                                                                                                                        >  2:418808/1‑251 (MQ=14)
                                                                                                               gTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGAttt                                                                                                                  <  2:732141/251‑1 (MQ=14)
                                                                                                                 gATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTcc                                                                                                                <  1:171513/251‑1 (MQ=14)
                                                                                                                 gATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTc                                                                                                                 <  1:353959/250‑1 (MQ=14)
                                                                                                                 gATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTc                                                                                                                 <  2:502436/250‑1 (MQ=14)
                                                                                                                  aTGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGACAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGGGATATGAACCGTTATAACCGACGATTTCCg                                                                                                               >  1:677334/1‑251 (MQ=9)
                                                                                                                       tGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACAGTTATAACCGGCGATTTCCGAATgg                                                                                                          >  1:338988/1‑251 (MQ=14)
                                                                                                                              tttGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGACTGGGGAAAcc                                                                                                   >  2:388937/1‑251 (MQ=14)
                                                                                                                                        taacaTCTGAATCAAGCTGAAAATTGAAACACTGACCAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGACTCGTAAGCAACATCTTCGAGATGTAAGGTTAAGCGACTCAGCGTACACGGTGTATGCCCTGGCAGTCAGAGGCGATGAAGAACGTGCTAATCTGCGGTAAGCGTCGGTCAGGTGATATGCACCGTTATAACCGGCGATTTCCGAATgggaggtcggcagagagtc                                                                                         >  2:788961/5‑235 (MQ=11)
                                                                                                                                        aaaaaTCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAGGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGttt                                                                                         >  1:362742/1‑251 (MQ=14)
                                                                                                                                        aaaaaTCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGtt                                                                                          <  1:570988/250‑1 (MQ=14)
                                                                                                                                        aaaaaTCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGGTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCTGGTTGTAAGGTTAAGCGACTAAGCGTACTCGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATggg                                                                                                         <  1:788961/235‑1 (MQ=9)
                                                                                                                                            aTCTGGATCAAGCTGAAAATTGAAACACTAAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCACCACGATGATGAATCGTAAGAAACATCTTCGGATTGTAAGGTTAAGCGACTAAGCGTACACGGTGGCTGCCCTGTCAGTAAGAGGCGGGGAAGGGCGAGCTAATCTGCGATAAGCGTCGGTAAGGTGCTATGAACCGTTTAAACCGGCGATTTACGAATGGGGAAAACCAGTGTGTTTCGa                                                                                      >  2:300021/1‑250 (MQ=11)
                                                                                                                                                     aagttGAAAATTGAAACACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCATTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCa                                                                            <  1:520416/247‑1 (MQ=9)
                                                                                                                                                      aGCTGAAAATTGAAACTCTGAACCACTAAAGTTGTTCGTGAGTCTCTCAAATATTCGCAACACGATGATGAATCGTAAGAAACATCTTCTGGTTGTAAGGTTAAGCGACTAAGCGTACAAGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATACGCGCCGGTAAGGTGATATGTATCGTTATAACCGTCGATTTCCGAATGGGGAAACCCAGTGTGGTGCGACACACTATCAt                                                                           <  2:99728/251‑1 (MQ=11)
                                                                                                                                                      aGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCAt                                                                           >  1:640199/1‑251 (MQ=14)
                                                                                                                                                          gAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAAc                                                                       <  1:618096/251‑1 (MQ=14)
                                                                                                                                                            aaaTTGAAACACTGAACAACGAAAGTTGTTCGTGAGGCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAGCATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTg                                                                     <  2:781946/251‑1 (MQ=14)
                                                                                                                                                               ttGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTgcat                                                                  >  1:349098/1‑248 (MQ=14)
                                                                                                                                                                 gAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATcc                                                                <  1:449217/251‑1 (MQ=14)
                                                                                                                                                                       cTGAACAACGAAAGTTGTTCGTGAGTCTATCAAATTTTCGCAACAAGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAAACAAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGt                                                          >  1:296919/1‑251 (MQ=9)
                                                                                                                                                                         gAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTa                                                        <  1:132055/251‑1 (MQ=14)
                                                                                                                                                                          aacaacGAAAGTTGTTCGTTAGTCTTTCAAATTTTCGAAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTaa                                                       <  2:726561/251‑1 (MQ=9)
                                                                                                                                                                          aacaacGAAAGGTGTTCGTGAGGTTCTCAAATTGTCGCAACACGTTGATGAACCGTAAGAAACATCTGCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGTCTATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGGTATAACCTGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCACATGTTaa                                                       <  2:327019/251‑1 (MQ=11)
                                                                                                                                                                             aacGAAAATTGTTCGTGAGTCTCTCAAATTTTCGCCACCCGAAGATACATAGTAAGACACATCTTAGGGTTGTAAGGTTAAGCGACTCAGCGTACACGAAGGATGCACTGGCAGTCAGAGGAGATGCAGGACATGCTAATCTACGATAAGCGTAGGTAAAGTGATATGAACCGTTATAACCGGCGATTTAAGACAAGAGAAACCCAGGGAGTTTTGACACACTAACATTAACTGAATCCATCGGTTaaat                                                     >  2:808709/1‑248 (MQ=2)
                                                                                                                                                                              acGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAg                                                   <  1:449900/251‑1 (MQ=14)
                                                                                                                                                                                gAAAGTTGTTCGTGAGTCTCGCAAATTTTCGCACCACGATGATGAATCGTAAGAAACATCTTCAGGTTGTAAGGTTAAGCGACTAAGAGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGAAGTGCTAATATGCGATAAGCATCAGTAGGGGGATATGAACCGTTATAACAGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACCATATCATTCCTTGAAACCATAGGTGAATgtgg                                                  >  2:432512/1‑247 (MQ=11)
                                                                                                                                                                                 aaaGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCg                                                >  1:369749/1‑251 (MQ=14)
                                                                                                                                                                                   aGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGTTTGTAAGGTTAAGCGACTAAGCGTACACGTTGGATGCCCTGGCAGTCAGATGCGATGa                                                                                                                                                                                >  1:105671/1‑121 (MQ=11)
                                                                                                                                                                                    gttgttTGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGCAAGATCTTCGGGTTGTAAGGTTCAGCGACTAAGCGTACACGGCGGATGCCCTGGCAGTCATAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGAGTTCCGAATGGTGAAACCCAGGGTGTTGCGACACACTATCATTACTTGAATCCATAGGTTAATGAGGCGaaa                                             <  2:86858/251‑2 (MQ=11)
                                                                                                                                                                                      tgttCGTGAGTCTCTCAAATTTTCGCAACACGATGATGTATCGTAAGCAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTAAGAGACGATGAAGGACGTGCTATTCTGCGATACGCGTCGGTAAGGAGATATGAACCGTTATAACCGGCGTTTTCCGACTAGGGAAACCCAGAGTGTTTTGACACACTATCATTAAATGAATCCATAGGTGAATGTGGcacaccg                                           >  2:204180/1‑245 (MQ=11)
                                                                                                                                                                                       gttCGTGAGTCTCTCAACTTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGGTAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGCTTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACCATCATTAACTGAATCCATAGGTTCATGAGGCggcccgg                                          >  2:421344/1‑245 (MQ=11)
                                                                                                                                                                                       gttCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCgg                                          >  1:580815/1‑251 (MQ=14)
                                                                                                                                                                                             gAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTACGGGTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAAAGGGGAAACCCAGTGTGTTTCGACACCCTATCATTAACTGAATCCATAGGTTAATAAGGCGAACCAGGGGAAc                                    >  2:792123/1‑251 (MQ=11)
                                                                                                                                                                                              aGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTAAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACt                                   <  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)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                   
AGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCACCTGCCTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGCGTCGGTAAGGTGATATGAACCGTTATAACCGGCGATTTCCGAATGGGGAAACCCAGTGTGTTTCGACACACTATCATTAACTGAATCCATAGGTTAATGAGGCGAACCGGGGGAACTGAAACATCTAAGTACCCCGAGGAAAAGAAATCA  >  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: