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

AAATGTAAGCGTTATCTTTATGCTAATATGGAAGTTAACGGACAACGATTTTCTTATGGAATTACCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTTATATAGTCTGGTAATGAGTTTT  >  NC_000913/4479874‑4480344
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                         
aaaTGTAAGCGTTATCTTTATGCTAATATGGAAGTTAACGGTCAACGATTTTCTTATGGAATTACCGGTGGGACTGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGGGTCCGCGTTTATTCCGCTTGTTGGTTGTGTTTTGACTTTGTTATTGTTGGTTTTGCTTAATTTTATGGCAGCAAAAGGAC‑TTCGTCTTCAGGGCTTAATGACCATtc                                                                                                                                                                                                                                           >  2:225141/1‑238 (MQ=255)
aaaTGTAAGCGTTATCTTTATGCTAATATGGAAGTTAACGGACAACGCTTTTCTTATGGAATTACCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCCATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCCTCCGGCCTTACTGACCAGtc                                                                                                                                                                                                                                           >  1:39972/1‑238 (MQ=255)
   tttAAGCGTTATCTTTATGTTAATATGGAAGTTAACGGAAAACGTTTTTTTTATTGAATTACCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCaa                                                                                                                                                                                                                                       <  2:944187/237‑1 (MQ=255)
   tGTAAGCGTTATCTTTATGCTAATATGGAAGTTAACGGACAACGATTTTCTTATGGAATTACCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTTTCa                                                                                                                                                                                                                                        >  1:523853/1‑238 (MQ=255)
     ttaGCGTTTTCTTTTTGCTAATATGGAAGTTAACGGACAACGATTTTCTTATGGAATTACCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACg                                                                                                                                                                                                                                     <  2:845484/237‑1 (MQ=255)
       aGCGTTATCTTTATGCTAATATGGAAGTTAACGGACAACGATTTTCTTATGGAATTACCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTCCTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAT‑TGCGTCATAAGGCTTTTATGACCAGTTTCAAAggg                                                                                                                                                                                                                                   >  2:866665/1‑238 (MQ=255)
             aTCTTTATGCTAATATGGAAGTTAACGGACAACGATTTTCTTATGGAATTACCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGCCTTTGTTACTGTTGGTTTTGCTTTTATTTATGGCAGCAAAAGGCCTTGGGCCAT‑GGGCCTTAATGAACAATTTCAACggggaaaaa                                                                                                                                                                                                                             >  2:309440/1‑232 (MQ=255)
             aTCTTTATGCTAATATGGAAGTTAACGGACAACGATTTTCTTATGGAATTACCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTTTATTTATGGCAGGCAAAGGAC‑TTGTTTTTCAGGGTTTATTGACCAGTTTTAACggggaaaaa                                                                                                                                                                                                                             >  2:711026/1‑232 (MQ=255)
              tttttATGCTATTTTGGATGTTTAAGGATAATGTTTTTCTTATGGAATTACCGGTGGGAATTGTTTTGTTAGTTGTCTTTTTGTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGAt                                                                                                                                                                                                                            <  2:492780/237‑1 (MQ=255)
               ccttATGCTAATATGGAATTTAACGGACAACGATTTTCTTATGGAATTTCCGGTTGGAATGTTTTTGTTAGTTGTCTTGTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGAtt                                                                                                                                                                                                                           <  1:825572/237‑1 (MQ=255)
                   aTGCTAATATGGAAGTTAACGGACAACGATTTTCTTATGGAATTACCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTACTCAGGCCTTAATTACCAATTTCAACGGCGTGAGATTTatt                                                                                                                                                                                                                       >  2:391885/1‑237 (MQ=255)
                   aTGCTAATATGGAAGTTAACGGACAACGATTTTCTTATGGAATTACCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAAGAC‑TGCGTCTTCAGGCCTTAATTACCAGTTTCAAAGGCGAAAAATTTatt                                                                                                                                                                                                                       >  2:170188/1‑237 (MQ=255)
                    tGCTAATATGGAAGTTAACGGACAACGATTTTCTTATGGAATTACCGGTGGGAATGTTTTTGTTAGTTGACTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGTTATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCCAAAGGAC‑TGCGCCATCAGGCCTTAATGACCATTTTCCACGGCGCAAGAtttgatt                                                                                                                                                                                                                      >  2:230475/1‑235 (MQ=255)
                        aaTATGGAAGTTAACGGACAACGATTTTCTTATGGAATTACCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCTGCAAAAGGAC‑TGCGTCATCAGGCCTTAATTACCAGTCTCAACGGCGTAAGATTTAGTTTTaa                                                                                                                                                                                                                  >  1:762047/1‑239 (MQ=255)
                         atatGGAAGTTAACGGACAACGATTTTCTTATGGAATTACCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCTGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAAt                                                                                                                                                                                                                 >  1:884146/1‑239 (MQ=255)
                         atatGGAAGTTAACGGACAACGATTTTCTTATGGAATTACCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAAt                                                                                                                                                                                                                 <  1:801704/239‑1 (MQ=255)
                         atatGGAAGTTAACGGACAACGATTTTCTTATGGAATTACCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑AGCGTCATCAGGCCTTAATGACCAGACTAAACGGCGTACGATTTCATTTGAAt                                                                                                                                                                                                                 >  1:720775/1‑239 (MQ=255)
                            tGGAATTTAAGGGACAGTGATTTTTTTATGGAATTGCGGGTGGGAATGTTTTTGTTAGTGGTCTTTTTTTTGTTTTTGTCTATTTCGAAATCTAAATGCCAGTGCGAGCAGATTTGCCGCTTGTTGGTTGGGTTTTGCCTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAG‑AAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGct                                                                                                                                                                                                             <  1:534864/239‑1 (MQ=255)
                             ggAAGTTAACGGACAACGATTTTCTTATGGAATTACCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAAATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGct                                                                                                                                                                                                             <  1:664413/239‑1 (MQ=255)
                                 gTTAACGGACAACGATTTTCTTATGGAATTACCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTa                                                                                                                                                                                                          <  1:309879/238‑1 (MQ=255)
                                   tAACGGACAACGATTTTCTTATGGAATTACCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGTGTCATTAGGCCTTAATGACCAGTTCCAACGGCGGAAGATTTAGTTTTAATTGCTCTCTGa                                                                                                                                                                                                       >  2:578097/1‑239 (MQ=255)
                                     acgtaaaagtgttttttttTGGAATTTCCGGTTGGAATTTTTTTGTTAGTTGTCCTTTTTTTTGTTGTTTGTATTTGGGTATCTTAATGACAGTGGCAGCAGATATGCCGGTTGTTGGTTGTGTTTTGGCTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGaaa                                                                                                                                                                                                     <  2:123379/228‑1 (MQ=255)
                                           aaCGATTTTCTTATGGAATTACCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTTTGGCAGAAAAAGGAC‑TTCTTCTTCTGGCCTTATTGACCAGTTTCAACGGGGGAAAATTTGATTTTAATTGTCTTATAAAAAGGttt                                                                                                                                                                                               >  2:27952/1‑238 (MQ=255)
                                            aCGATTTTCTTATGGAATTACCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCCACGGCGTAAGATTTAGTTTTAATTGCTCTTTGAAAGGGt                                                                                                                                                                                                 >  1:580423/1‑236 (MQ=255)
                                                   ttttttGGATTTTCCGGTTGGAATTTTTTTGTTATTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTggtgggt                                                                                                                                                                                       <  2:310832/234‑1 (MQ=255)
                                                               aCCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCTGCTTGATGGTTGTGTTTTGCCTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCATTAATGACCAGTATCACCGGCGTAAGATTTAGTTTTCATTTCTTCATGAACGGGTTTTGGTGGGTGCCCtttttttt                                                                                                                                                                           >  1:615740/1‑239 (MQ=255)
                                                                     gggAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATTACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTTCTCTTTGAAAAGGTTTTGGTGGGGGACCTTTTTTTTGCCGAt                                                                                                                                                                     >  2:280131/1‑239 (MQ=255)
                                                                      ggAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATCTTTATGGCAGCAAAAGAAC‑TTCTTCATTAGGCCTTAATTACCTGTCTCCACGGCTTAAGGTTTAGTGTTAATGTCTCTATAAAAGAGGTCTGGTGGGGGACCTTTTTTTTGcgtttt                                                                                                                                                                    >  2:535065/1‑234 (MQ=255)
                                                                          ttttttttttAGTTGTTTTTTTTTTTTTTTTTTTTTTTTCGCCATCTTAATGTAAGTGTCAGCATATATGCCGCTTTTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTa                                                                                                                                                                 <  2:625008/236‑1 (MQ=255)
                                                                           ttttttGTTAGTTGGCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAAt                                                                                                                                                               <  1:83038/238‑1 (MQ=255)
                                                                              ttttttAGTTGTCTTTTTTTTGTTTTTTTCTATTTCCCAATTGTAATGACAGTGTGAGCAGTTATGCCGCGTTTTGGTTGTGTTTTGGCTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGc                                                                                                                                                            <  2:662622/239‑1 (MQ=255)
                                                                                  ttagtttgtttttttttGTTTTTTTCTATTTCGCAATCTTAATGACAGGGTTAGCAGATATGCCGTTTGTTGGTTGTGTTTTGACTTTGTTTCTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCAtt                                                                                                                                                        <  1:819760/230‑1 (MQ=255)
                                                                                  ttAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCAtt                                                                                                                                                        <  1:360098/239‑1 (MQ=255)
                                                                                         tCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGGGCCCTTTTTCTTTCCGATTTTAATGGCCATTGGGTTgg                                                                                                                                                 >  1:490887/1‑239 (MQ=255)
                                                                                            tttttttGTTTTTTTCTATTTCGAATTTTTAATGTCAGTGTCTGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGa                                                                                                                                              <  2:568103/239‑1 (MQ=255)
                                                                                                   gTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGGCTCAACGGCGTAAGATTTAGGTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTTGGATtgggcagctttt                                                                                                                                       >  1:568103/1‑231 (MQ=255)
                                                                                                                gcGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACa                                                                                                                          <  1:391885/238‑1 (MQ=255)
                                                                                                                      tCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTaa                                                                                                                                                                                                                  <  1:364995/145‑1 (MQ=255)
                                                                                                                         gaaTGACAGTGTCAGCAGATATGCCGCTTGTTGGTTTTGTTTTGACTTTGTTACTTTTGGTTTGGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTAAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACTTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTa                                                                                                                 <  1:569640/238‑1 (MQ=255)
                                                                                                                          aaTGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAAGGTTCTGGGGGGGGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACGGTTTTTTTTTTTTCGAAAAAAATGATa                                                                                                                 >  2:450288/1‑238 (MQ=255)
                                                                                                                                aGTGTCTGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCa                                                                                                          <  1:866665/239‑1 (MQ=255)
                                                                                                                                aGTGTCAGGAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAT‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCa                                                                                                          <  1:809001/239‑1 (MQ=255)
                                                                                                                                aGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGTGGGGTGACCTTTTTCTTGCCGTTTTTAATGGCCATTGGGTTGGGGGCGGTTTTTTTTATTTCTACAAAAATGTTCCCTGCCa                                                                                                          >  2:569640/1‑239 (MQ=255)
                                                                                                                                aGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGCCCTTTTTCTTGCCGATTTTAATGGCCATTGGGTTGGGGACTGTTTTCTTTTTCTCGACAAAAATTCTACCttcac                                                                                                          >  2:809001/1‑235 (MQ=255)
                                                                                                                                  tttCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTa                                                                                                                                                                 <  2:1060322/180‑1 (MQ=255)
                                                                                                                                  tgtCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCaa                                                                                                                                                                                                                                       <  2:307740/112‑1 (MQ=255)
                                                                                                                                  tgtCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCaa                                                                                                                                                                                                                                       >  1:307740/1‑112 (MQ=255)
                                                                                                                                  tgtCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTggtggg                                                                                                                                                                                        <  1:543277/159‑1 (MQ=255)
                                                                                                                                  tgtCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTggtggg                                                                                                                                                                                        >  2:543277/1‑159 (MQ=255)
                                                                                                                                  tgtCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGCTTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAAt                                                                                                        >  1:410450/1‑239 (MQ=255)
                                                                                                                                  tgtCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATTGCCATTTGGATGGGGACTGTTTTCTTTATTTCGACAAAGATGCTACCTGCCaaa                                                                                                        >  2:51374/1‑238 (MQ=255)
                                                                                                                                  tgtCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAAGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTa                                                                                                                                                                 >  1:1060322/1‑182 (MQ=255)
                                                                                                                                   tttgtgagATTTGCCGCTTGTTGGGTTGGTTTTGACTTTGTTACTGTGGGTTTTGCTTATATTTATGGAAGAAAAAGGAC‑TGCGTCATAAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAAt                                                                                                        <  1:66020/232‑1 (MQ=255)
                                                                                                                                           gtttGCCTCTTTTTGGTTGTGTTTTGGCTTTTTTCCTTTTGTTTTTGTTTTTATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAAGGACCATTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGCCCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTCCCTGCCAATAGTTCAag                                                                                                <  2:762047/235‑1 (MQ=255)
                                                                                                                                            tttGCCGTTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAag                                                                                                <  1:871817/235‑1 (MQ=255)
                                                                                                                                            tatGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAag                                                                                                <  1:170188/237‑1 (MQ=255)
                                                                                                                                            tatGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAag                                                                                                <  1:280131/237‑1 (MQ=255)
                                                                                                                                             aggctgCTTGGTGGTTGGGTTTGGACTTTGTTCCTGTTGGTTTTGATTATATTTAGGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAagtag                                                                                             <  1:873878/234‑1 (MQ=255)
                                                                                                                                                   cTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGTTGGGGCCTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGttat                                                                                       >  1:1000806/1‑239 (MQ=255)
                                                                                                                                                    ttgttgGTTGTGTTTTGCCTTTGTTCCTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGttatt                                                                                      <  1:866930/239‑1 (MQ=255)
                                                                                                                                                    aggttttttgTGTTTTGTCTTGATTTCTGTGGTTTTTCTTTTTTTTTATGGCAGAAAAAGGTC‑TGGGTAATTAGGGTTAAAGGACCAGTCTGAAAGGCTAAAAATTTAGTTTTAATTGCTCTATTAAAGGGTTTGGGTGGGTGACCTTTTTCGTGCGGATTTTAATGGCGATTGGGAGGGGGAATGTTTTCTTTATCTCGACAAAGATGCTCCCTGCCAATAGTTCAAGTAGTGttatt                                                                                      <  2:615740/232‑1 (MQ=255)
                                                                                                                                                     tgttgGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGGAGGGttttta                                                                                     >  1:435721/1‑239 (MQ=255)
                                                                                                                                                      tttgtttgtTTTTTGGTTTTGTTTCTGTTGTTTTTGCTTATTTTTAAGGCAGCAAAAGGAT‑GGCTTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTtattat                                                                                    <  2:435721/238‑1 (MQ=255)
                                                                                                                                                      gtgtgttgttTTTTGACTTTGTTACTTTTGGTTTTGCTTATATTTATGGCAGCAAAAGGTC‑TGCGTCATCAGGCCTTAATGACCAGTCTAAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTTTtattat                                                                                    <  1:534187/235‑1 (MQ=255)
                                                                                                                                                           tttggtTTTGACTTTGTTACGTTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCa                                                                               <  1:309440/235‑1 (MQ=255)
                                                                                                                                                                ttttGGTTTTGTTACTGTTGGTTTTGCTTATATTTAGGGAAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGtt                                                                          <  1:1053209/239‑1 (MQ=255)
                                                                                                                                                                  ttGACTTTGTTACTCTTGGTTTTGATTATATTAAGGGAACAAAAAGGTC‑TGCGGAATCAGGCCTTAATGACCAGGTTCAACGGCTGAAGATTTAGTTTTAATTGCTATAAGAAAGGGTTCTGGTGGGTGCCCTTTTTCATTCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCCACAAAGATGGTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCt                                                                        <  1:201014/239‑1 (MQ=255)
                                                                                                                                                                  ttGACTTTGTTACCTGTGGTTTTGCTTATATGTCTGGAAGCAAAAGGTC‑TGCGTCATCGGGACTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGAGGTCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCt                                                                        <  1:225141/239‑1 (MQ=255)
                                                                                                                                                                  ttGACTGTGTTACTGTTGGTTTTGGTTATATTTAGGGCAGCAAAGGGAT‑TGCGTCATCAGGCCTTAATTACCAGTCGAAACGGCGTAAGATTTAGTTTTAATTGCTCTAGGAAAGGGTTCTGGGGGGTGACTTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTGCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCt                                                                        <  1:329670/239‑1 (MQ=255)
                                                                                                                                                                     aCTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGCCATCAGGCCTTAATGACCAGTCTCCACGGAGTAAGATTTAGTCTTACTTGCTCTATGAAAGGGTTCTGGGGGGTGACCTTTTTCTTGCCGATTTTTATGGCCCTTGGGTTGGTGACTGTTTTCTTTTTCCCGCCAAAGATTATACCTGCCCATCGTTCACTAAGCGTTATTATATCCAAGGGCCTGAt                                                                     >  2:471347/1‑239 (MQ=255)
                                                                                                                                                                        ttGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGGGTTATTATATCCATGGTTCTGATGGc                                                                  >  1:361502/1‑239 (MQ=255)
                                                                                                                                                                        ttGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCACCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTg                                                                                                                                                                                                               >  1:205876/1‑98 (MQ=255)
                                                                                                                                                                        ttGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCACCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTg                                                                                                                                                                                                               <  2:205876/98‑1 (MQ=255)
                                                                                                                                                                         ttttttttttGTTTTTGCTTATTTTTAGGGCAGTAAAAGGAC‑TTGGTCATCAGGCCTTAATGACCAGTCTCAACGGGGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACTTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGAATGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCa                                                                 <  2:1000806/231‑1 (MQ=255)
                                                                                                                                                                           ttACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTAAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAAt                                                               <  1:815404/239‑1 (MQ=255)
                                                                                                                                                                           gtACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTGAGTTTTAATGGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAAt                                                               <  1:578097/238‑1 (MQ=255)
                                                                                                                                                                              cTGTTGGTTTTTATTTTTTTTATGGCAGAAAAAGGAT‑TGGGTAATCAGGCCTTAATGACCATTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTCTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGt                                                            <  2:580423/239‑1 (MQ=255)
                                                                                                                                                                                   ggTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATTTCGGCAAAAATGCCACCTGCCCATATGTCAGTTAGTGTTATTTTATCCATGGGTTTGATGGGAATAATTggat                                                       >  2:892179/1‑237 (MQ=255)
                                                                                                                                                                                         gCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCCATAGTTCAAGTAGGGTTATTATATCCATGGTTCTTATGGCAATAGTTGGGATTGtt                                                  >  2:33109/1‑238 (MQ=255)
                                                                                                                                                                                              tatTTATGGCAGTAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCAtt                                            <  1:711026/239‑1 (MQ=255)
                                                                                                                                                                                              tatTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCTTTAAGGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCAGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGtt                                                                                                    <  1:201002/183‑1 (MQ=255)
                                                                                                                                                                                                    tGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGGCAAAGATGCTACCTGCCCCTAGTTCAAGTCGGGTTTTTATATCCATGGTTCTGTTGGCAATAGTTTGAATTTTTTCCATTGGTAtt                                      >  2:70379/1‑239 (MQ=255)
                                                                                                                                                                                                         gcaAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTaa                                 >  1:579191/1‑239 (MQ=255)
                                                                                                                                                                                                           aaatGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATg                               <  1:814353/239‑1 (MQ=255)
                                                                                                                                                                                                           aaaaGGAC‑TGTGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATg                               <  2:579191/239‑1 (MQ=255)
                                                                                                                                                                                                           aaaaGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATg                               <  1:27952/239‑1 (MQ=255)
                                                                                                                                                                                                           aaaaAGAA‑TGAGTCATCAGGCATTAATGACGAGGTTAAACGGCGTAAGATTTAGTTTTAATTGCTTTAGGAAAGGGATCGGGTGGGTACGCGTTTTCTGGCCGATTTTAAGGGCCATTGGGATGGGTACTGTTTTCTTTATATCGAAAAAGAGGTTACCTGCCAATAGTTCAAGTAGTGTTATGATATCCATGGTTCTGATGGCAATAGTTGGTATAGTTTCCATTGGTATTTTTAATg                               <  2:720775/239‑1 (MQ=255)
                                                                                                                                                                                                            aaaGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATg                               <  1:450288/238‑1 (MQ=255)
                                                                                                                                                                                                                 aC‑TGGGTCATCAGGCCTTAATGACCAGTCTCAACGGGGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTAc                           <  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)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                         
AAATGTAAGCGTTATCTTTATGCTAATATGGAAGTTAACGGACAACGATTTTCTTATGGAATTACCGGTGGGAATGTTTTTGTTAGTTGTCTTTTTTTTGTTTTTTTCTATTTCGCAATCTTAATGACAGTGTCAGCAGATATGCCGCTTGTTGGTTGTGTTTTGACTTTGTTACTGTTGGTTTTGCTTATATTTATGGCAGCAAAAGGAC‑TGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTTATATAGTCTGGTAATGAGTTTT  >  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: ‑