Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,480,111 C→T 11.2% L128F (CTC→TTC)  yjgN → DUF898 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,480,1110CT11.2% 246.8 / 3.6 96L128F (CTC→TTC) yjgNDUF898 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (30/53);  new base T (9/3);  total (40/56)
Fisher's exact test for biased strand distribution p-value = 2.44e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

AAATGTAAGCGTTATCTTTATGCTAATATGGAAGTTAACGGACAACGATTTTCTTATGGAATTACCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTTATATAGTCTGGTAATGAGTTTT  >  NC_000913/4479874‑4480344
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                         
aaaTGTAAGCGTTATCTTTATGCTAATATGGAAGTTAACGGTCAACGATTTTCTTATGGAATTACCGGTGGGACTGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGGGTCCGCGTTTATTCCGCTTGTTGGTTGTGTTTTGACTTTGTTATTGTTGGTTTTGCTTAATTTTATGGCAGCAAAAGGACTTCGTCTTCAGGGCTTAATGACCATtc                                                                                                                                                                                                                                           >  2:225141/1‑238 (MQ=255)
aaaTGTAAGCGTTATCTTTATGCTAATATGGAAGTTAACGGACAACGCTTTTCTTATGGAATTACCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCCATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCCTCCGGCCTTACTGACCAGtc                                                                                                                                                                                                                                           >  1:39972/1‑238 (MQ=255)
   tttAAGCGTTATCTTTATGTTAATATGGAAGTTAACGGAAAACGTTTTTTTTATTGAATTACCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCaa                                                                                                                                                                                                                                       <  2:944187/237‑1 (MQ=255)
   tGTAAGCGTTATCTTTATGCTAATATGGAAGTTAACGGACAACGATTTTCTTATGGAATTACCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTTTCa                                                                                                                                                                                                                                        >  1:523853/1‑238 (MQ=255)
     ttaGCGTTTTCTTTTTGCTAATATGGAAGTTAACGGACAACGATTTTCTTATGGAATTACCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACg                                                                                                                                                                                                                                     <  2:845484/237‑1 (MQ=255)
       aGCGTTATCTTTATGCTAATATGGAAGTTAACGGACAACGATTTTCTTATGGAATTACCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTCCTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGATTGCGTCATAAGGCTTTTATGACCAGTTTCAAAggg                                                                                                                                                                                                                                   >  2:866665/1‑238 (MQ=255)
             aTCTTTATGCTAATATGGAAGTTAACGGACAACGATTTTCTTATGGAATTACCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGCCTTTGTTACTGTTGGTTTTGCTTTTATTTATGGCAGCAAAAGGCCTTGGGCCATGGGCCTTAATGAACAATTTCAACggggaaaaa                                                                                                                                                                                                                             >  2:309440/1‑232 (MQ=255)
             aTCTTTATGCTAATATGGAAGTTAACGGACAACGATTTTCTTATGGAATTACCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTTTATTTATGGCAGGCAAAGGACTTGTTTTTCAGGGTTTATTGACCAGTTTTAACggggaaaaa                                                                                                                                                                                                                             >  2:711026/1‑232 (MQ=255)
              tttttATGCTATTTTGGATGTTTAAGGATAATGTTTTTCTTATGGAATTACCGGTGGGAATTGTTTTGTTAGTTGTCTTTTTGTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGAt                                                                                                                                                                                                                            <  2:492780/237‑1 (MQ=255)
               ccttATGCTAATATGGAATTTAACGGACAACGATTTTCTTATGGAATTTCCGGTTGGAATGTTTTTGTTAGTTGTCTTGTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGAtt                                                                                                                                                                                                                           <  1:825572/237‑1 (MQ=255)
                   aTGCTAATATGGAAGTTAACGGACAACGATTTTCTTATGGAATTACCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTACTCAGGCCTTAATTACCAATTTCAACGGCGTGAGATTTatt                                                                                                                                                                                                                       >  2:391885/1‑237 (MQ=255)
                   aTGCTAATATGGAAGTTAACGGACAACGATTTTCTTATGGAATTACCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAAGACTGCGTCTTCAGGCCTTAATTACCAGTTTCAAAGGCGAAAAATTTatt                                                                                                                                                                                                                       >  2:170188/1‑237 (MQ=255)
                    tGCTAATATGGAAGTTAACGGACAACGATTTTCTTATGGAATTACCGGTGGGAATGTTTTTGTTAGTTGACTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGTTATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCCAAAGGACTGCGCCATCAGGCCTTAATGACCATTTTCCACGGCGCAAGAtttgatt                                                                                                                                                                                                                      >  2:230475/1‑235 (MQ=255)
                        aaTATGGAAGTTAACGGACAACGATTTTCTTATGGAATTACCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCTGCAAAAGGACTGCGTCATCAGGCCTTAATTACCAGTCTCAACGGCGTAAGATTTAGTTTTaa                                                                                                                                                                                                                  >  1:762047/1‑239 (MQ=255)
                         atatGGAAGTTAACGGACAACGATTTTCTTATGGAATTACCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCTGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAAt                                                                                                                                                                                                                 >  1:884146/1‑239 (MQ=255)
                         atatGGAAGTTAACGGACAACGATTTTCTTATGGAATTACCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAAt                                                                                                                                                                                                                 <  1:801704/239‑1 (MQ=255)
                         atatGGAAGTTAACGGACAACGATTTTCTTATGGAATTACCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACAGCGTCATCAGGCCTTAATGACCAGACTAAACGGCGTACGATTTCATTTGAAt                                                                                                                                                                                                                 >  1:720775/1‑239 (MQ=255)
                            tGGAATTTAAGGGACAGTGATTTTTTTATGGAATTGCGGGTGGGAATGTTTTTGTTAGTGGTCTTTTTTTTGTTTTTGTCTATTTCGAAATCTAAATGCCAGTGCGAGCAGATTTGCCGCTTGTTGGTTGGGTTTTGCCTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGct                                                                                                                                                                                                             <  1:534864/239‑1 (MQ=255)
                             ggAAGTTAACGGACAACGATTTTCTTATGGAATTACCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAAATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGct                                                                                                                                                                                                             <  1:664413/239‑1 (MQ=255)
                                 gTTAACGGACAACGATTTTCTTATGGAATTACCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTa                                                                                                                                                                                                          <  1:309879/238‑1 (MQ=255)
                                   tAACGGACAACGATTTTCTTATGGAATTACCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGTGTCATTAGGCCTTAATGACCAGTTCCAACGGCGGAAGATTTAGTTTTAATTGCTCTCTGa                                                                                                                                                                                                       >  2:578097/1‑239 (MQ=255)
                                     acgtaaaagtgttttttttTGGAATTTCCGGTTGGAATTTTTTTGTTAGTTGTCCTTTTTTTTGTTGTTTGTATTTGGGTATCTTAATGACAGTGGCAGCAGATATGCCGGTTGTTGGTTGTGTTTTGGCTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGaaa                                                                                                                                                                                                     <  2:123379/228‑1 (MQ=255)
                                           aaCGATTTTCTTATGGAATTACCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTTTGGCAGAAAAAGGACTTCTTCTTCTGGCCTTATTGACCAGTTTCAACGGGGGAAAATTTGATTTTAATTGTCTTATAAAAAGGttt                                                                                                                                                                                               >  2:27952/1‑238 (MQ=255)
                                            aCGATTTTCTTATGGAATTACCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCCACGGCGTAAGATTTAGTTTTAATTGCTCTTTGAAAGGGt                                                                                                                                                                                                 >  1:580423/1‑236 (MQ=255)
                                                   ttttttGGATTTTCCGGTTGGAATTTTTTTGTTATTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTggtgggt                                                                                                                                                                                       <  2:310832/234‑1 (MQ=255)
                                                               aCCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCTGCTTGATGGTTGTGTTTTGCCTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCATTAATGACCAGTATCACCGGCGTAAGATTTAGTTTTCATTTCTTCATGAACGGGTTTTGGTGGGTGCCCtttttttt                                                                                                                                                                           >  1:615740/1‑239 (MQ=255)
                                                                     gggAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATTACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTTCTCTTTGAAAAGGTTTTGGTGGGGGACCTTTTTTTTGCCGAt                                                                                                                                                                     >  2:280131/1‑239 (MQ=255)
                                                                      ggAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATCTTTATGGCAGCAAAAGAACTTCTTCATTAGGCCTTAATTACCTGTCTCCACGGCTTAAGGTTTAGTGTTAATGTCTCTATAAAAGAGGTCTGGTGGGGGACCTTTTTTTTGcgtttt                                                                                                                                                                    >  2:535065/1‑234 (MQ=255)
                                                                          ttttttttttAGTTGTTTTTTTTTTTTTTTTTTTTTTTTCGCCATCTTAATGTAAGTGTCAGCATATATGCCGCTTTTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTa                                                                                                                                                                 <  2:625008/236‑1 (MQ=255)
                                                                           ttttttGTTAGTTGGCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAAt                                                                                                                                                               <  1:83038/238‑1 (MQ=255)
                                                                              ttttttAGTTGTCTTTTTTTTGTTTTTTTCTATTTCCCAATTGTAATGACAGTGTGAGCAGTTATGCCGCGTTTTGGTTGTGTTTTGGCTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGc                                                                                                                                                            <  2:662622/239‑1 (MQ=255)
                                                                                  ttagtttgtttttttttGTTTTTTTCTATTTCGCAATCTTAATGACAGGGTTAGCAGATATGCCGTTTGTTGGTTGTGTTTTGACTTTGTTTCTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCAtt                                                                                                                                                        <  1:819760/230‑1 (MQ=255)
                                                                                  ttAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCAtt                                                                                                                                                        <  1:360098/239‑1 (MQ=255)
                                                                                         tCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGGGCCCTTTTTCTTTCCGATTTTAATGGCCATTGGGTTgg                                                                                                                                                 >  1:490887/1‑239 (MQ=255)
                                                                                            tttttttGTTTTTTTCTATTTCGAATTTTTAATGTCAGTGTCTGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGa                                                                                                                                              <  2:568103/239‑1 (MQ=255)
                                                                                                   gTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGGCTCAACGGCGTAAGATTTAGGTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTTGGATtgggcagctttt                                                                                                                                       >  1:568103/1‑231 (MQ=255)
                                                                                                                gcGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACa                                                                                                                          <  1:391885/238‑1 (MQ=255)
                                                                                                                      tCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTaa                                                                                                                                                                                                                  <  1:364995/145‑1 (MQ=255)
                                                                                                                         gaaTGACAGTGTCAGCAGATATGCCGCTTGTTGGTTTTGTTTTGACTTTGTTACTTTTGGTTTGGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTAAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACTTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTa                                                                                                                 <  1:569640/238‑1 (MQ=255)
                                                                                                                          aaTGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAAGGTTCTGGGGGGGGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACGGTTTTTTTTTTTTCGAAAAAAATGATa                                                                                                                 >  2:450288/1‑238 (MQ=255)
                                                                                                                                aGTGTCTGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCa                                                                                                          <  1:866665/239‑1 (MQ=255)
                                                                                                                                aGTGTCAGGAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGATTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCa                                                                                                          <  1:809001/239‑1 (MQ=255)
                                                                                                                                aGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGTGGGGTGACCTTTTTCTTGCCGTTTTTAATGGCCATTGGGTTGGGGGCGGTTTTTTTTATTTCTACAAAAATGTTCCCTGCCa                                                                                                          >  2:569640/1‑239 (MQ=255)
                                                                                                                                aGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGCCCTTTTTCTTGCCGATTTTAATGGCCATTGGGTTGGGGACTGTTTTCTTTTTCTCGACAAAAATTCTACCttcac                                                                                                          >  2:809001/1‑235 (MQ=255)
                                                                                                                                  tttCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTa                                                                                                                                                                 <  2:1060322/180‑1 (MQ=255)
                                                                                                                                  tgtCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCaa                                                                                                                                                                                                                                       <  2:307740/112‑1 (MQ=255)
                                                                                                                                  tgtCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCaa                                                                                                                                                                                                                                       >  1:307740/1‑112 (MQ=255)
                                                                                                                                  tgtCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTggtggg                                                                                                                                                                                        <  1:543277/159‑1 (MQ=255)
                                                                                                                                  tgtCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTggtggg                                                                                                                                                                                        >  2:543277/1‑159 (MQ=255)
                                                                                                                                  tgtCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGCTTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAAt                                                                                                        >  1:410450/1‑239 (MQ=255)
                                                                                                                                  tgtCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATTGCCATTTGGATGGGGACTGTTTTCTTTATTTCGACAAAGATGCTACCTGCCaaa                                                                                                        >  2:51374/1‑238 (MQ=255)
                                                                                                                                  tgtCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAAGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTa                                                                                                                                                                 >  1:1060322/1‑182 (MQ=255)
                                                                                                                                   tttgtgagATTTGCCGCTTGTTGGGTTGGTTTTGACTTTGTTACTGTGGGTTTTGCTTATATTTATGGAAGAAAAAGGACTGCGTCATAAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAAt                                                                                                        <  1:66020/232‑1 (MQ=255)
                                                                                                                                           gtttGCCTCTTTTTGGTTGTGTTTTGGCTTTTTTCCTTTTGTTTTTGTTTTTATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAAGGACCATTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGCCCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTCCCTGCCAATAGTTCAag                                                                                                <  2:762047/235‑1 (MQ=255)
                                                                                                                                            tttGCCGTTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAag                                                                                                <  1:871817/235‑1 (MQ=255)
                                                                                                                                            tatGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAag                                                                                                <  1:170188/237‑1 (MQ=255)
                                                                                                                                            tatGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAag                                                                                                <  1:280131/237‑1 (MQ=255)
                                                                                                                                             aggctgCTTGGTGGTTGGGTTTGGACTTTGTTCCTGTTGGTTTTGATTATATTTAGGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAagtag                                                                                             <  1:873878/234‑1 (MQ=255)
                                                                                                                                                   cTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGTTGGGGCCTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGttat                                                                                       >  1:1000806/1‑239 (MQ=255)
                                                                                                                                                    ttgttgGTTGTGTTTTGCCTTTGTTCCTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGttatt                                                                                      <  1:866930/239‑1 (MQ=255)
                                                                                                                                                    aggttttttgTGTTTTGTCTTGATTTCTGTGGTTTTTCTTTTTTTTTATGGCAGAAAAAGGTCTGGGTAATTAGGGTTAAAGGACCAGTCTGAAAGGCTAAAAATTTAGTTTTAATTGCTCTATTAAAGGGTTTGGGTGGGTGACCTTTTTCGTGCGGATTTTAATGGCGATTGGGAGGGGGAATGTTTTCTTTATCTCGACAAAGATGCTCCCTGCCAATAGTTCAAGTAGTGttatt                                                                                      <  2:615740/232‑1 (MQ=255)
                                                                                                                                                     tgttgGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGGAGGGttttta                                                                                     >  1:435721/1‑239 (MQ=255)
                                                                                                                                                      tttgtttgtTTTTTGGTTTTGTTTCTGTTGTTTTTGCTTATTTTTAAGGCAGCAAAAGGATGGCTTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTtattat                                                                                    <  2:435721/238‑1 (MQ=255)
                                                                                                                                                      gtgtgttgttTTTTGACTTTGTTACTTTTGGTTTTGCTTATATTTATGGCAGCAAAAGGTCTGCGTCATCAGGCCTTAATGACCAGTCTAAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTTTtattat                                                                                    <  1:534187/235‑1 (MQ=255)
                                                                                                                                                           tttggtTTTGACTTTGTTACGTTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCa                                                                               <  1:309440/235‑1 (MQ=255)
                                                                                                                                                                ttttGGTTTTGTTACTGTTGGTTTTGCTTATATTTAGGGAAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGtt                                                                          <  1:1053209/239‑1 (MQ=255)
                                                                                                                                                                  ttGACTTTGTTACTCTTGGTTTTGATTATATTAAGGGAACAAAAAGGTCTGCGGAATCAGGCCTTAATGACCAGGTTCAACGGCTGAAGATTTAGTTTTAATTGCTATAAGAAAGGGTTCTGGTGGGTGCCCTTTTTCATTCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCCACAAAGATGGTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCt                                                                        <  1:201014/239‑1 (MQ=255)
                                                                                                                                                                  ttGACTTTGTTACCTGTGGTTTTGCTTATATGTCTGGAAGCAAAAGGTCTGCGTCATCGGGACTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGAGGTCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCt                                                                        <  1:225141/239‑1 (MQ=255)
                                                                                                                                                                  ttGACTGTGTTACTGTTGGTTTTGGTTATATTTAGGGCAGCAAAGGGATTGCGTCATCAGGCCTTAATTACCAGTCGAAACGGCGTAAGATTTAGTTTTAATTGCTCTAGGAAAGGGTTCTGGGGGGTGACTTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTGCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCt                                                                        <  1:329670/239‑1 (MQ=255)
                                                                                                                                                                     aCTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGCCATCAGGCCTTAATGACCAGTCTCCACGGAGTAAGATTTAGTCTTACTTGCTCTATGAAAGGGTTCTGGGGGGTGACCTTTTTCTTGCCGATTTTTATGGCCCTTGGGTTGGTGACTGTTTTCTTTTTCCCGCCAAAGATTATACCTGCCCATCGTTCACTAAGCGTTATTATATCCAAGGGCCTGAt                                                                     >  2:471347/1‑239 (MQ=255)
                                                                                                                                                                        ttGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGGGTTATTATATCCATGGTTCTGATGGc                                                                  >  1:361502/1‑239 (MQ=255)
                                                                                                                                                                        ttGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCACCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTg                                                                                                                                                                                                               >  1:205876/1‑98 (MQ=255)
                                                                                                                                                                        ttGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCACCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTg                                                                                                                                                                                                               <  2:205876/98‑1 (MQ=255)
                                                                                                                                                                         ttttttttttGTTTTTGCTTATTTTTAGGGCAGTAAAAGGACTTGGTCATCAGGCCTTAATGACCAGTCTCAACGGGGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACTTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGAATGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCa                                                                 <  2:1000806/231‑1 (MQ=255)
                                                                                                                                                                           ttACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTAAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAAt                                                               <  1:815404/239‑1 (MQ=255)
                                                                                                                                                                           gtACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTGAGTTTTAATGGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAAt                                                               <  1:578097/238‑1 (MQ=255)
                                                                                                                                                                              cTGTTGGTTTTTATTTTTTTTATGGCAGAAAAAGGATTGGGTAATCAGGCCTTAATGACCATTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTCTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGt                                                            <  2:580423/239‑1 (MQ=255)
                                                                                                                                                                                   ggTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATTTCGGCAAAAATGCCACCTGCCCATATGTCAGTTAGTGTTATTTTATCCATGGGTTTGATGGGAATAATTggat                                                       >  2:892179/1‑237 (MQ=255)
                                                                                                                                                                                         gCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCCATAGTTCAAGTAGGGTTATTATATCCATGGTTCTTATGGCAATAGTTGGGATTGtt                                                  >  2:33109/1‑238 (MQ=255)
                                                                                                                                                                                              tatTTATGGCAGTAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCAtt                                            <  1:711026/239‑1 (MQ=255)
                                                                                                                                                                                              tatTTATGGCAGCAAAAGGACTGCGTCATCAGGCTTTAAGGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCAGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGtt                                                                                                    <  1:201002/183‑1 (MQ=255)
                                                                                                                                                                                                    tGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGGCAAAGATGCTACCTGCCCCTAGTTCAAGTCGGGTTTTTATATCCATGGTTCTGTTGGCAATAGTTTGAATTTTTTCCATTGGTAtt                                      >  2:70379/1‑239 (MQ=255)
                                                                                                                                                                                                         gcaAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTaa                                 >  1:579191/1‑239 (MQ=255)
                                                                                                                                                                                                           aaatGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATg                               <  1:814353/239‑1 (MQ=255)
                                                                                                                                                                                                           aaaaGGACTGTGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATg                               <  2:579191/239‑1 (MQ=255)
                                                                                                                                                                                                           aaaaGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATg                               <  1:27952/239‑1 (MQ=255)
                                                                                                                                                                                                           aaaaAGAATGAGTCATCAGGCATTAATGACGAGGTTAAACGGCGTAAGATTTAGTTTTAATTGCTTTAGGAAAGGGATCGGGTGGGTACGCGTTTTCTGGCCGATTTTAAGGGCCATTGGGATGGGTACTGTTTTCTTTATATCGAAAAAGAGGTTACCTGCCAATAGTTCAAGTAGTGTTATGATATCCATGGTTCTGATGGCAATAGTTGGTATAGTTTCCATTGGTATTTTTAATg                               <  2:720775/239‑1 (MQ=255)
                                                                                                                                                                                                            aaaGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATg                               <  1:450288/238‑1 (MQ=255)
                                                                                                                                                                                                                 aCTGGGTCATCAGGCCTTAATGACCAGTCTCAACGGGGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTAc                           <  1:70379/237‑1 (MQ=255)
                                                                                                                                                                                                                      cGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTtata                     >  2:318143/1‑239 (MQ=255)
                                                                                                                                                                                                                             aGGCCTTAATGACCAGTCTCAACGGCTTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTTATATAGTCTg              <  1:318143/239‑1 (MQ=255)
                                                                                                                                                                                                                              ggCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCt                                                                        <  1:772004/180‑1 (MQ=255)
                                                                                                                                                                                                                              ggCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCt                                                                        >  2:772004/1‑180 (MQ=255)
                                                                                                                                                                                                                                ccTTAAGGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTGGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCAGTGGTATTTTTAATGGTACTTTATATAGTCTGGt            <  1:903315/238‑1 (MQ=255)
                                                                                                                                                                                                                                    aaTGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTTATATAGTCTGGTAATgg       >  1:726203/1‑238 (MQ=255)
                                                                                                                                                                                                                                      agACCAGGTTTAACGGCGTTAGTGTTAGTTTTAACGTGTTGAGGAAAGGGTTCTGGTGGGTTACCTTTTTCTTGCCGATTTTAATGGCCATGGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGTTACCTGACAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGGTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTTATATAGTCTGGTAATGAGt     <  1:471347/238‑1 (MQ=255)
                                                                                                                                                                                                                                         ccAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACACAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGCTGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTTATATTGTCTGGTAATGAGtttt  >  1:1028722/1‑239 (MQ=255)
                                                                                                                                                                                                                                          cAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGACAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCACGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTCATGGTACTTTATTTAGTCTGGTATTGAGttt   >  1:724487/1‑237 (MQ=255)
                                                                                                                                                                                                                                             tctcAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACt                          >  1:932908/1‑211 (MQ=255)
                                                                                                                                                                                                                                             tctcAACGGCGTAAGATTTAGTTTTAATTGCTCTAGGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACt                          <  2:932908/211‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                         
AAATGTAAGCGTTATCTTTATGCTAATATGGAAGTTAACGGACAACGATTTTCTTATGGAATTACCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTTATATAGTCTGGTAATGAGTTTT  >  NC_000913/4479874‑4480344

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