Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,058,486 T→C 25.8% S135P (TCG→CCG)  torC → trimethylamine N‑oxide (TMAO) reductase I, cytochrome c‑type subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,058,4860TC25.8% 152.6 / 13.0 98S135P (TCG→CCG) torCtrimethylamine N‑oxide (TMAO) reductase I, cytochrome c‑type subunit
Reads supporting (aligned to +/- strand):  ref base T (48/17);  new base C (6/21);  total (56/42)
Fisher's exact test for biased strand distribution p-value = 6.44e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATT  >  NC_000913/1058252‑1058710
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                    
gtgtATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCggg                                                                                                                                                                                                                                 >  2:579422/1‑238 (MQ=255)
gtgtATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACCCGGc                                                                                                                                                                                                                                 >  2:46321/1‑239 (MQ=255)
gtgtATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGGGAATGGGCGCGAAGGAAAAAAAAAAACCCgga                                                                                                                                                                                                                                 >  1:753114/1‑238 (MQ=255)
gtgtATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGGAAATGGGGCGCAAAGAAAGAAACAAACCCgga                                                                                                                                                                                                                                 >  1:406214/1‑238 (MQ=255)
gtgtATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAACCAAACGCACGGAACTTGCCGCGCGCGAATGGGCCGGAATGAAAGAAAACACATCGGc                                                                                                                                                                                                                                 >  2:298492/1‑239 (MQ=255)
gtgtATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGGGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACCCCCGGAACTTGCCGAGCTGGAAAGGGGCCGAATGAAAAAAAAAAAACCGGc                                                                                                                                                                                                                                 >  1:819629/1‑239 (MQ=255)
      gaagaaTATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGGCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTg                                                                                                                                                                                                                         <  2:165505/239‑1 (MQ=255)
          aaTATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCg                                                                                                                                                                                                                      >  2:810246/1‑238 (MQ=255)
           atatAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACACGGCACCCTcccca                                                                                                                                                                                                                     >  2:450706/1‑236 (MQ=255)
           atatAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAAACCTTTATTGCTCACTCCATTGATACACCTGCAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGCGAATGGGCGCGAATGAAAGACAACAAACCGGCAACCTGCCGc                                                                                                                                                                                                                     >  2:223485/1‑238 (MQ=255)
               aaaCAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGCGAAGGGGCGGGAATGAAAGAAAACAAACCGGCAACCTGCGGccccg                                                                                                                                                                                                                 >  1:403054/1‑237 (MQ=255)
                 aCAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTTCCGAGCGTGAATGGGCGCGAATGAAAAAAAACAACTCGGCAAACTGCCGcctcggcc                                                                                                                                                                                                              >  2:795235/1‑232 (MQ=255)
                    gTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCCAACGCCCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGCGAACCTGCCCCACCTGCCATa                                                                                                                                                                                                           >  2:438960/1‑239 (MQ=255)
                    gTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGACTGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATa                                                                                                                                                                                                           >  2:592657/1‑239 (MQ=255)
                      cGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCGGCCATAAc                                                                                                                                                                                                         >  2:770994/1‑239 (MQ=255)
                      cGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCCCCTGCCaaaaa                                                                                                                                                                                                         >  1:272190/1‑238 (MQ=255)
                      cGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGAGTGAATGGGCCCGAATGAAAGCAAACAACTCGGCAACCTGCCGCCCCTGCCCTAAc                                                                                                                                                                                                         >  2:550795/1‑239 (MQ=255)
                         tGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCctacctac                                                                                                                                                                                                      >  2:614943/1‑239 (MQ=255)
                          gcctttCCAGAACGCCCCCGGCGCGGGAGGTGAATGCCATGACGTTCATATCCCGCGGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGAAATGATATCTACCAGACCTTTATTGCTCACTCCAGTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGGGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTAcg                                                                                                                                                                                                     <  2:406214/236‑1 (MQ=255)
                           cctttCCAGACCGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGGGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTAcg                                                                                                                                                                                                     <  2:744385/236‑1 (MQ=255)
                               tCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATGTTCCAGGCATGGTGAAGCGCAAACTGGAAGCGCGCAATGATATCTACCAGACCTTCATTGCTCACTCCATTGATACCCCTGAAAAATTCGAAGCCAAACGCACGGAACTTGCCGAGCGTGAATGGACGCGCCTGAACGAAAAAAACTCGGCACCCTCCCGCCCCTTccctcacaccgtatag                                                                                                                                                                                                >  1:388548/1‑225 (MQ=255)
                               tCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACCCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAAATCGGCAACCTGCCGCTCCTGCCATAACTAcgacacg                                                                                                                                                                                                >  2:741359/1‑235 (MQ=255)
                                   gAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTCCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAAATCGGCAACCTGCCGACCCTGCCATAACACcgatgcaaaga                                                                                                                                                                                            >  1:166902/1‑234 (MQ=255)
                                     aCGCCTCCGGCGTGCGAGATGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTGCATTGATACACCTGAACAATTCGAAGCCAAAAGCGCGGAACTTGCCAAGCGTGACTGGGCGCGCATGAAAGCAAAACAACAGGCAACATGCCGATCCTGCCATAAATAcacgacggtggat                                                                                                                                                                                          >  1:389275/1‑227 (MQ=255)
                                                 tggggtcTGAATGCCTTGACTTTCATACCCCGCCGGATTTCCCAGGCATGGTGAAGCGAAAACTGGAGGGGAACAATGATATCTACCAGACCTTTATTGCTCACCCCATTGATACCCCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTAAATGGGCGCGAATAAAAGAAAACAACCCGGCACCCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAagca                                                                                                                                                                               <  1:116506/232‑1 (MQ=255)
                                                  gCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAagcag                                                                                                                                                                              <  2:753114/238‑1 (MQ=255)
                                                   gcggtgtaaatccatgcatctCATATCCCGCGCGGATATCCCGGGCTGGGGGAGGCGAAAACGGGAAGGGAGCAATAATACCTACCAGCCTTTTATTGCTCACTCCATGGATACACCGGAAAAATTCGAAGCAAAACGCGGGGAACTTGCCGAGGGGAAAGGGGCGCGAATGAAAAAAAACAACCCGGCACCCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAagcag                                                                                                                                                                              <  1:685091/218‑1 (MQ=255)
                                                   cGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAACCAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGACCATCGGCagcagc                                                                                                                                                                             >  2:816045/1‑238 (MQ=255)
                                                         gaatgccggttcTGTCAAACCCCGCGGGATTTTCAAGGCGTGGTGAAGCGAAAACTGGAAGCGAGAAAGAAGATCCACCAGGCTTTTATTGCTGCCTCCATTGATACACCTGAAAAATGCGAGCCAAAACGCGGGGAACTTGCCGGGCGAGAATGGGCGCGAATGAAAGAAACCAACCCGGCAACCGGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGa                                                                                                                                                                      <  1:96742/228‑1 (MQ=255)
                                                                acgACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAAAAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGGTCAGGCGAAGCAGCAACCTGAAacaaca                                                                                                                                                               >  1:139436/3‑239 (MQ=255)
                                                                  gacgttctgtccgccccgGATTTCCCACGCTTGGGGAAGCGCAAATGGGAAAGAAGAAAAAACATCGCCCAGTCCTTTATTGCTCACACAGTTGTTCCCCCGCAAAATTTCGAGTCCAACCCCGCGGACTTGGCCGGAGTGAAAGGGGGGCGAAGGAAAGAAACCAACCCGGGACCCTGCAGCGCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACg                                                                                                                                                             <  1:509113/224‑1 (MQ=255)
                                                                  gACTGGCATAGCCCGCCGGATATTCCAGGGATGGTGAAGCGGAAAGTGGAAGCGAACAATGATATCTACCAGACCTTTATTGCTCACTCCAGTGATCCACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGGGTGAATGGGCGGGAATTAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACg                                                                                                                                                             <  2:381553/239‑1 (MQ=255)
                                                                       tCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGa                                                                                                                                                       >  2:550296/1‑239 (MQ=255)
                                                                       tCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGCCAGa                                                                                                                                                       >  2:54270/1‑239 (MQ=255)
                                                                               cgccgGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCCAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGACTGAAAGAAAACAACTCGGCAACCTGCCGCCCCTGCCATAACTACGATGCGAGGGCACAAGCGAACGAGACTCCTAAAACACAACGTAAGAGGaaaggg                                                                                                                                               >  1:116369/1‑235 (MQ=255)
                                                                               cgccgGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGAGGGATCATGGGAAGCACAACCCTGAAGCAGCCCGCCAGAGGAAgggg                                                                                                                                               >  1:714701/1‑237 (MQ=255)
                                                                                 ccgGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGCTCATGCGAAGCAGCATCCTGAAGCAGCACGCCAGCGAAAGGTgg                                                                                                                                              >  2:192100/1‑238 (MQ=255)
                                                                                 ccgGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGGGGc                                                                                                                                             >  2:202391/1‑239 (MQ=255)
                                                                                 ccgGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAACCAGCAGTCAGGCTGAAGGGGGc                                                                                                                                             >  2:771310/1‑239 (MQ=255)
                                                                                 ccgGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAACAAGCACGCCAGATGAAGGGGGc                                                                                                                                             >  2:761282/1‑239 (MQ=255)
                                                                                 ccgGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTCCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCGAAACTACGATCCGATGGATAATGCGAAGAAGCATCCTGAACCACAACGTCAGAACAAGGGGGc                                                                                                                                             >  1:234833/1‑239 (MQ=255)
                                                                                 ccgGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAAAGCGCGGAACTTGCCGAGCGTTAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACCACGATGCGATGGATCATGCAAACCACCACCCCGAACCAGCACCCCAAATGAAGGTgga                                                                                                                                             >  1:786038/1‑238 (MQ=255)
                                                                                         ttCCCGGTATGGTGAAGCACAACTTGGAGGGAAAAAAATATATTTACCAGCCTTTTATTCCTCACCCCTGTGATCACCCGAAAAAATTCGAACCAAACCGCGGGGAATTTGCCGGGCGTAAAGGGGCGGGAATAAAAAAAAACACCCCGGCACCCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGaaa                                                                                                                                      <  1:750836/239‑1 (MQ=255)
                                                                                               gCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGTTaat                                                                                                                                >  2:252174/1‑239 (MQ=255)
                                                                                               gCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGAGGCAGGTCAGATGAAGGGGGCAGCGAAAGAAaat                                                                                                                                >  1:765208/1‑239 (MQ=255)
                                                                                                cATGTTGAAGCGCAATCTGGAAGCGACCATTTATTTCTTCCAGACCTTTATTGCTCACTCCATTGATACCCCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTTAATGGCCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATaat                                                                                                                                <  2:765208/238‑1 (MQ=255)
                                                                                                    gTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATaatcaat                                                                                                                            >  2:223934/1‑238 (MQ=255)
                                                                                                             aaaCTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACCCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGa                                                                                                                  <  2:139436/239‑1 (MQ=255)
                                                                                                              aaCTGGAAGGGAGCAATTATATCTCCAAGACCTTTTTGCCTCCCCCAATGGACAACCCTAAAAATTTAGAACCAAAACGCGCGAAACGTGCCGGGCGTGAATGGGGGAGAAAGAAAAAAAACAACCCGGAAACCTGCCCCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGAc                                                                                                                 <  1:592657/239‑1 (MQ=255)
                                                                                                                    aaggaggacttgttatCTACCAACCCTTTGTTGGCTACTCCTTGGAACCCCCGAAAAAATTGGAAGCAAACCGCGGGGACTTGGGCGGGCGGAAAGGGGCGGAAATAAAAGAAAACAACCCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCa                                                                                                            <  1:395758/225‑1 (MQ=255)
                                                                                                                     acggaaagataatatCCACCGGGCTTTTATCTCCCCCTCCATTGATATCCCAAAAAAATCCGAAGCCAAACCCGGGGACTTGGCCCGGGGGGAATGGCGGAAATGAAAAAAAAACATCCCGGCACCCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATa                                                                                                          <  1:539775/228‑1 (MQ=255)
                                                                                                                         aGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGACAGATAACCAATCCTGCATCGATCGTCCTAAAgg                                                                                                      >  1:179476/1‑239 (MQ=255)
                                                                                                                            aaTGATATCTACCAGACCTTTATTGCTCACCCCATTGATCCACCTAAAAAATGCGAAGCCAACCGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTa                                                                                                    <  2:403054/238‑1 (MQ=255)
                                                                                                                                 tatCTACCAGACTTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGGGCGTGAAGGGGCGCGAATGAAAGAAAACAACCCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGcc                                                                                               <  2:655490/238‑1 (MQ=255)
                                                                                                                                   ttcacacctctttttttGCTCACTCCAGTGAAACACATAAAAAATTGGAACCCAAACGCGGGGGAGTGGCCGGGGGGAAATGGGGGCAAATGAAAGAAACACACTGGGCACCCTGCCCCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCccac                                                                                            <  1:653490/228‑1 (MQ=255)
                                                                                                                                   tCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGCAGGCGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGAATTGCccac                                                                                            >  2:398497/1‑239 (MQ=255)
                                                                                                                                        cAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCCGtt                                                                                       >  2:37472/1‑239 (MQ=255)
                                                                                                                                          gACCTTTATTGCTCACTCCATTGATACACCTGAAAAATCCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGACGGCAGCGAAAGATAATCAATCCTGCATCGCCGTACATAAAGGAATTGCCCACCAGttc                                                                                      >  2:630675/1‑237 (MQ=255)
                                                                                                                                           aCCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGAATTGCCCACCAGTTAcc                                                                                    >  2:542582/1‑239 (MQ=255)
                                                                                                                                                     cTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGGGAAAGGTAATCAATCCTGCACCGACTGTCATACAGGTATTGCCCCCCAGTTACCGGatatagct                                                                          >  2:156203/1‑235 (MQ=255)
                                                                                                                                                       cACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGAt                                                                                                                                                      >  2:183484/1‑161 (MQ=255)
                                                                                                                                                       cACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTAc                                                                                     <  2:8008/226‑1 (MQ=255)
                                                                                                                                                       cACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTAAc                                                                                     >  1:8008/1‑226 (MQ=255)
                                                                                                                                                       cACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCACATGAAGGGGGCAGCGAACGATAATCAATCCTGCACCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGAAAGGAgaac                                                                        >  1:717423/1‑236 (MQ=255)
                                                                                                                                                          tCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCgg                                                                     >  2:253258/1‑239 (MQ=255)
                                                                                                                                                          tCCATTGAAACACCTAAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGGGCGTGAATGGGGGCGAATGAAAGAAAACAACCCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCgg                                                                     <  2:272190/239‑1 (MQ=255)
                                                                                                                                                          ccctttttatacccTAAAAAATTGAACACAAACCGCCCGGAATTGGCCGAGGGTAAAGGGGCGCAAAGAAAAAAAACAACCTGGGACACCTCCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCgg                                                                     <  1:761282/227‑1 (MQ=255)
                                                                                                                                                          acttttGATACCCTTAAAAATTTGAACGCCAAACGCGGGAAATTGGCCAGGGGGGAAGGGGCGGAAATGAAAAAAAACCACCCGGCACCCGGCCCCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCgg                                                                     <  1:46321/235‑1 (MQ=255)
                                                                                                                                                                    cacCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAagcag                                                                                                                                                                              <  2:39909/125‑1 (MQ=255)
                                                                                                                                                                    cacCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAagcag                                                                                                                                                                              >  1:39909/1‑125 (MQ=255)
                                                                                                                                                                    cacCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAg                                                           >  2:328048/1‑239 (MQ=255)
                                                                                                                                                                    cacCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTCAATCACCTTTCATAAAGGAATTGCCCACCAGTTACCGGAAATGATAAGGGGGTTCCGTAAg                                                           >  1:25247/1‑239 (MQ=255)
                                                                                                                                                                    cacCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCGGCATCGCCTGTCATAAAGGAATAGCCCCCCAGTTACCGGGTATAAGAAACGAGTTCCGAAAg                                                           >  1:339701/1‑239 (MQ=255)
                                                                                                                                                                     cacattaaaaaTTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACCACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGAAGGCAGCGTAAGATAATCACTCCTTCATCGACTGTCATAACGGTATTCCCCACCCATTACCGGATATAAGAAACGGGTTGCGTaaa                                                           >  1:862047/7‑238 (MQ=255)
                                                                                                                                                                      ccTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATCAAGGAATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCa                                                         >  2:897609/1‑239 (MQ=255)
                                                                                                                                                                      ccTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCa                                                         >  2:482192/1‑239 (MQ=255)
                                                                                                                                                                      ccTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGGATTGcc                                                                                               >  2:40456/1‑201 (MQ=255)
                                                                                                                                                                      ccTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGGGGCAGCGAAAGATAATCAACCCGGCATCGACTGTCATAAAGGTATTGCCCACCAATAACCGGACATGAGTAGCGGCCTCCggacaca                                                         >  1:346107/1‑233 (MQ=255)
                                                                                                                                                                                      acccccaccgccgtgAACTTGCCGGGGGTGAGTGGGGGGGAAAAAAAAAAACCAGCTCGGCACCCTGCCGCCCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGGATTGCCCACCAGTTGCCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTg                                            <  1:770994/225‑1 (MQ=255)
                                                                                                                                                                                           aaggtgagggagtttgcagtggttGATTGGGCGAGAAGAAAGGAACCTAACTCGCAACCTGGCCTCTCCGGCCATAACGACGATGTGATGGATCATGCGAAGCAGCATCCTGAAGCAGCAAGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTgcgcgc                                       <  1:833568/216‑1 (MQ=255)
                                                                                                                                                                                           aaacggcgggAACTTGCCGGGGTGTAATGGGCCCAAAATAAAAAAAACCACCCGGCAACCTGCCGCCCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTgcgcgc                                       <  2:679444/231‑1 (MQ=255)
                                                                                                                                                                                             ccgcgcgGAACTGGCCGAGCGTGAAGGGGCGGGAAGGAAAGAAAACACCGCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCa                                     <  2:346107/238‑1 (MQ=255)
                                                                                                                                                                                                   ggAACTTGCCTGGCGTGAATGGGCACAAATGAAAGAAAACAACCCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCt                                <  2:166902/238‑1 (MQ=255)
                                                                                                                                                                                                   gaaactgccgggggttGAAGGGCCGCGAATGAAAGAAACAACCTCGTCACCCTGCCGCTCCTTCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCt                                <  2:211330/223‑1 (MQ=255)
                                                                                                                                                                                                       cTTGCCGAGGGGAAATGGGCGAGAATGAAAAAAAACAACCGGGCAACAGCCGGCCCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGa                           <  2:193628/239‑1 (MQ=255)
                                                                                                                                                                                                          gCCGGGCGGGAAGGGGCGGGAAGGAAAAAAACCAACCCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAg                        <  2:685050/239‑1 (MQ=255)
                                                                                                                                                                                                          gCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGAATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAACTGCGCGCCCGGGTTAATGAAAg                        >  1:685050/1‑239 (MQ=255)
                                                                                                                                                                                                                   ggaagggggcagAATGAAAGAAAACAACGCGGCAACCTGGCGCTCCTGCCATAAGTACGATGCGATGGATCAGGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGac                 <  2:388548/227‑1 (MQ=255)
                                                                                                                                                                                                                     ggaaggggcgaaaagtaaagaaaCCACCCCGGCACCCGGACGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGac                 <  2:849145/222‑1 (MQ=255)
                                                                                                                                                                                                                      aTGGGCCGAAAGAAAAAAAACACCCCGGGCACCCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACTGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACg              <  1:328048/238‑1 (MQ=255)
                                                                                                                                                                                                                           tggggcgaaaagaaaaaaaCCACCCGGCCACACGGCCGCCCCGGCCATAAATACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGGACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCt            <  1:771310/224‑1 (MQ=255)
                                                                                                                                                                                                                           cgcgAATGAAAGAAAACACCTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTa         <  2:25247/238‑1 (MQ=255)
                                                                                                                                                                                                                           cgcgAATAAAAAAAACCAACCGGGCAACCTGCCGACCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTa         <  2:786038/238‑1 (MQ=255)
                                                                                                                                                                                                                           cgcaAATAAAAGAAAACACCTCGGAACCCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTa         <  2:714701/238‑1 (MQ=255)
                                                                                                                                                                                                                             gccgaaacaaaaaaaaCACTGCGGCCCCCTGCCCCTCCTGCCACAACTACGATGCGATGGAACAGGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTa         <  1:54270/230‑1 (MQ=255)
                                                                                                                                                                                                                              gcgacgaaaaaaaaaCTACTCGTCACCCTCCCGCTCCTGCCCACACTACGAGGCGTTGGATCATGCGAAGCAGCATCCTGAAACAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTAct       <  1:220607/232‑1 (MQ=255)
                                                                                                                                                                                                                                   taaaaaaaaaCAACCCGGCCACCCGCCGCCCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTAtt  <  1:202391/237‑1 (MQ=255)
                                                                                                                                                                                                                                   aaaaaaaaaaCCACGCGGCACCCCGCCGCCCCGGCCAAAACTACGGTGCGAGGGAGCATGCGAAGAAGCATCCCTAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTAtt  <  1:706536/237‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                    
GTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTATT  >  NC_000913/1058252‑1058710

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