Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,420,658 A→C N227H (AAC→CAC)  ulaA → L‑ascorbate‑specific enzyme IIC permease component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,420,6580AC96.7% 296.2 / ‑1.1 90N227H (AAC→CAC) ulaAL‑ascorbate‑specific enzyme IIC permease component of PTS
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base C (49/38);  total (51/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

CTCCATGTGGACCACCATTATCTGTACCGCGATTCTGGTTTCGCTCTACTGGGGCATCACTTCCAACATGATGTACAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGAACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGACACCGTGCAGGCGATGGCAGGCAAAGTGCACTGGACGGTGTACATCCTGCAAACTGGTTTCTCCTTTGCGGTGGCGATCTTCATCATCACGCAGGGTGTGCGCATGTTTGTGGCGGAA  >  NC_000913/4420438‑4420858
                                                                                                                                                                                                                            |                                                                                                                                                                                                        
cTCCATGTGGACCACCATTATCTGTACCGCGATTCTGGTTTCGCTCTACTGGGGCATCACTTCCAACATGATGTACAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGac                                                                                                                                                                                            >  2:19903/1‑235 (MQ=255)
accaattgggACCACCTTTATCTGTACCGCGATTCTGGTTTCGCTCTACTGGGGCATCACTTCCAACATGATGTACAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGac                                                                                                                                                                                            <  2:298279/227‑1 (MQ=255)
          accaccATTATCTGTACCGCGATTCTGGTTTCGCTCTACTGGGGCATCACTTCCAACATGATGTACAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGtctc                                                                                                                                                                                 >  1:593157/1‑236 (MQ=255)
          accaccATTATCTGTACCGCGATTCTGGTTTCGCTCTACTGGGGCATCACTTCCAACATGATGTACAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGtctc                                                                                                                                                                                 >  1:465370/1‑236 (MQ=255)
           ccaccaTTATCTGTACCGCGATTCTGGTTTCGCTCTACTGGGGCATCACTTCCAACATGATGTACAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTcc                                                                                                                                                                                >  1:458502/1‑236 (MQ=255)
            caccaTTATCTGTACCGCGATTCTGGTTTCGCTCTACTGGGGCATCACTTCCAACATGATGTACAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAAAGCGCTGAAGACCTCAAATTGCCGTGCTGGCTGCCAAATTTCCACGACAACAGCGtc                                                                                                                                                                                   >  2:209709/1‑232 (MQ=255)
               caTTATCTGTACCGCGATTCTGGTTTCGCTCTACTGGGGCATCACTTCCAACATGATGTACAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACgg                                                                                                                                                                            >  1:196995/1‑236 (MQ=255)
               caTTATCTGTACCGCGATTCTGGTTTCGCTCTACTGGGGCATCACTTCCAACATGATGTACAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACgg                                                                                                                                                                            >  1:469820/1‑236 (MQ=255)
                 ttACCTGTACCGGGATTCTGGTTTCGCTTTACTGGGGCATCACTTCCAACATGATGTACAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCg                                                                                                                                                                          <  2:619285/236‑1 (MQ=255)
                   aTCTGTACCGCGATTCTGGTTTCGCTCTACTGGGGCATCACTTCCAACTTGATGTCCAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGAt                                                                                                                                                                        <  2:367632/236‑1 (MQ=255)
                      tgactcgggtttCGGGTTTCGCTCTACAGGGGCATCACTTCCAACATGATGTACAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATtgt                                                                                                                                                                     <  2:391312/226‑1 (MQ=255)
                      tGTACCGCGATTCTGGTTTCGCTCTACTGGGGCATCACTTCCAACATGATGTACAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGGCTCCACGGCGATtgt                                                                                                                                                                     >  1:425459/1‑236 (MQ=255)
                        tACCGCGATTCTGGTTTCGCTCTACTGGGGCATCACTTCCAACATGATGTACAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGACCCTTTTCCACGACAACATCGTATCCACGGCGAAtggg                                                                                                                                                                    >  2:660322/1‑233 (MQ=255)
                         aCCGCGATTCTGGTTTCGCTCTACTGGGGCATCACTTCCAACATGATGTACAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAGGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGtga                                                                                                                                                                   >  1:76761/1‑235 (MQ=255)
                         aCCGCGATTCTGGTTTCGCTCTACTGGGGCATCACTTCCAACATGATGTACAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGtga                                                                                                                                                                   >  1:391312/1‑235 (MQ=255)
                         aCCGCGATTCTGGTTTCGCTCTACTGGGGCATCACTTCCAACATGATGTACAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCAACGACAACATCGTTCTCACGGCGATTGtga                                                                                                                                                                   >  2:742135/1‑235 (MQ=255)
                           cgcgATTCTGGTTTCGCTCTACTGGGGCATCACTTCCAACATGATGTACAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACTAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGtgatga                                                                                                                                                                >  2:629733/1‑236 (MQ=255)
                             cgATTCTGGTTTCGCTCTACTGGGGCATCACTTCCAACATGATGTACAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAACAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACACCATCGTATCCACGGCGACTGtgatga                                                                                                                                                                >  2:149513/1‑234 (MQ=255)
                                ttCTGGTTTCGCTCTACTGGGGCATCACTTCCAACATGATGTACAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCAACTTTTCCACGACAACATCGTCGCCACGGCGCTTGTGATGACCa                                                                                                                                                             >  2:598047/1‑234 (MQ=255)
                                     gtttccttatACTGGGCAATCCCTTCAAACTTGATGTACAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATcttctt                                                                                                                                                      <  2:196995/227‑1 (MQ=255)
                                          gCTCTACTGGGGCATCACTTCCAACATGATGTACAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTg                                                                                                                                                 <  1:442994/236‑1 (MQ=255)
                                          gCTCTACTGGGGCATCACTTCCAACATGATGTACAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTg                                                                                                                                                 >  1:351524/1‑236 (MQ=255)
                                          gCTCTACTGGGGCATCACTTCCAACATGATGTACAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCGTCTTTGGTg                                                                                                                                                 >  2:442994/1‑236 (MQ=255)
                                          gCTCTACTGGGGAATCGCTTCCAACATGATGTACAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTg                                                                                                                                                 <  1:273307/236‑1 (MQ=255)
                                                    ggCATCACTTCCAACATGATGTACAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGAACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGc                                                                                                                                        >  2:151767/1‑235 (MQ=255)
                                                        tCACTTCCAACATGATGTACAAGCCGACTCAGGAAGTGACGGATGGTTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACAttt                                                                                                                                                                                                    <  1:537738/171‑1 (MQ=255)
                                                        tCACTTCCAACATGATGTACAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCTTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTACACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTg                                                                                                                                                 >  1:427220/1‑222 (MQ=255)
                                                        tCACTTCCAACATGATGTACAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCTTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTACACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTg                                                                                                                                                 <  2:427220/222‑1 (MQ=255)
                                                        tCACTTCCAACATGATGTACAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGcaca                                                                                                                                                                                                       >  1:749716/1‑168 (MQ=255)
                                                        tCACTTCCAACATGATGTACAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGcaca                                                                                                                                                                                                       <  2:749716/168‑1 (MQ=255)
                                                        tCACTTCCAACATGATGTACAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACAttt                                                                                                                                                                                                    >  2:537738/1‑171 (MQ=255)
                                                        tCACTTCCAACATGATGTACAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGCCCATCTTCTTTGGTGCCATTCTGCTCTcc                                                                                                                                   >  2:467351/1‑236 (MQ=255)
                                                        tCACTTCCAACATGATGTACAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTcc                                                                                                                                   >  1:132068/1‑236 (MQ=255)
                                                        tCACTTCCAACATGATGTACAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCGTCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCACCTTCTTTGGTGCCATTCTGCTCTcc                                                                                                                                   >  2:646100/1‑236 (MQ=255)
                                                          aCTTCCAACATGATGTACAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCt                                                                                                                                  <  1:660322/235‑1 (MQ=255)
                                                           cctCCCACTTGGTGTACAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTc                                                                                                                                <  2:564300/234‑1 (MQ=255)
                                                                     gagtaaCAGGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTATCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGaca                                                                                                                      <  2:469820/231‑1 (MQ=255)
                                                                      atgTACAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACAttt                                                                                                                                                                                                    <  2:85274/157‑1 (MQ=255)
                                                                      atgTACAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACAttt                                                                                                                                                                                                    >  1:85274/1‑157 (MQ=255)
                                                                      atgTACAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGacac                                                                                                                     >  1:88176/1‑236 (MQ=255)
                                                                        gTACAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGACACCg                                                                                                                   >  1:315711/1‑236 (MQ=255)
                                                                                cGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGACACCGTGCAGGCg                                                                                                           <  1:19903/236‑1 (MQ=255)
                                                                                 gACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATcttctt                                                                                                                                                      >  1:664774/1‑192 (MQ=255)
                                                                                 gACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATcttctt                                                                                                                                                      <  2:664774/192‑1 (MQ=255)
                                                                                 gACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTc                                                                                                                                <  1:323834/214‑1 (MQ=255)
                                                                                 gACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTc                                                                                                                                >  2:323834/1‑214 (MQ=255)
                                                                                 gACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCggtacccgtca                                                                                                               >  2:158071/1‑221 (MQ=255)
                                                                                    tCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGACACCGTGCAGGCGATgg                                                                                                       <  1:629733/236‑1 (MQ=255)
                                                                                                  aTGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGCTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGACACCGTGAAGTCTATGGCAGGCCAAGTGCAc                                                                                         >  1:284959/1‑236 (MQ=255)
                                                                                                   gagCTGTGGTTTCGCCAGCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGACACCGTGCAGGCGATGGCAGGCAAAGTGCACt                                                                                        <  2:593157/234‑1 (MQ=255)
                                                                                                      cTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGACACCGTGCAGGCGATGGCAGGCAAAGTGCACTGGa                                                                                     <  1:630719/236‑1 (MQ=255)
                                                                                                      cTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCGCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGACACCGTGCAGGCGATGGCAGGCAAAGTGCACTgg                                                                                      <  1:742135/235‑1 (MQ=255)
                                                                                                       tgtgGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAGGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGACACCGTGCAGGCGATGGCAGGCAAAGTGCACTGGa                                                                                     <  2:76761/235‑1 (MQ=255)
                                                                                                       tgtgGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGACACCGTGCAGGCGATGGCAGGCAAAGTGCACTGGa                                                                                     <  2:458502/235‑1 (MQ=255)
                                                                                                       tgtgGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGACACCGTGCAGGCGATGGCAGGCAAAGTGCACTGGa                                                                                     <  1:684143/235‑1 (MQ=255)
                                                                                                       tgtgGTTTCTCCATCGGCACCCCGCAGCAGTTTGCATCATGGATTGCCTAAAACGTCGCGCCGTTACTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGACACCGTGCAGGCGATGGCAGGCAAAGTGCACTGGa                                                                                     <  2:465370/235‑1 (MQ=255)
                                                                                                       tgtgGTTTCTCCAGCGGTCCACAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGACACCGTGCAGGCGATGGCAGGCAAAGTGCACTGGa                                                                                     <  2:229953/235‑1 (MQ=255)
                                                                                                           gTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGACACCGTGCAGGCGATGGCAGGCAAAGTGCACTGGACGgtg                                                                                >  1:76791/1‑236 (MQ=255)
                                                                                                            tttCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTg                                                                                                                                                 >  1:660685/1‑170 (MQ=255)
                                                                                                            tttCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTg                                                                                                                                                 <  2:660685/170‑1 (MQ=255)
                                                                                                                tcCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGACACCGTGCAGGCGATGGCAGGCAAATGGCACTGGACGGGGTa                                                                              >  2:142166/1‑233 (MQ=255)
                                                                                                                tcCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGACACCGTGCAGGCGATGGCAGGCAAAGTGCACTGGACGGTGTa                                                                              >  1:709624/1‑233 (MQ=255)
                                                                                                                 ccATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGac                                                                                                                                                                                            >  1:394410/1‑122 (MQ=255)
                                                                                                                 ccATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGac                                                                                                                                                                                            <  2:394410/122‑1 (MQ=255)
                                                                                                                     cGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGACACCGTGCAGGCGATGGCAGGCAAAGTGCACTGGACGGTGTACATCCTg                                                                       >  2:260324/1‑235 (MQ=255)
                                                                                                                          aCCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACgg                                                                                                                                                                            <  2:757499/129‑1 (MQ=255)
                                                                                                                          aCCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACgg                                                                                                                                                                            >  1:757499/1‑129 (MQ=255)
                                                                                                                                  cagTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGACACCGTGCAGGCGATGGCAGGCAAAGTGCACTGGACGGTGTACATCCTGCAAACTgggt                                                             >  1:364479/1‑230 (MQ=255)
                                                                                                                                      ttGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGACACCGTGCAGGCGATGGCAGGCAAAGTGCACTGGACGGGGTAAATCCGGCAAATCGGTTTTTCCttt                                                     >  2:99933/1‑236 (MQ=255)
                                                                                                                                             atGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGACACCGTGCAGGCGATGGCAGGCAAAGTGCACTGGACGGTGTACATCCTGCAAACTGGTTTCTCCTTTGCGGTgg                                              <  1:191046/236‑1 (MQ=255)
                                                                                                                                             atGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGAACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGACACCGTGCAGGCGATGGCAGGCAAAGTGCACTGGACGGTGTACATCCTGCAAACTGGTTTCTCCTTTGCGGTgg                                              <  1:151767/236‑1 (MQ=255)
                                                                                                                                                   tGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGACACCGTGCAGGCGATGGCAGGCAAAGTGCACTGGACGGTGTACATCCTGCAAACTGGTTTCTCCTTTGCGGTGGCGATCt                                        >  1:351222/1‑236 (MQ=255)
                                                                                                                                                        ataAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGGGCCATTCTGCTCTCCTTCGGTATCGACACCGTGCAGGCGATGGCAGGCAAAGTGCACTGGACGGTGTACATCCTGCAAACTGGTTTCTCCTTTGCGGTGGCGATCTtcatc                                   <  2:132068/236‑1 (MQ=255)
                                                                                                                                                         taAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGACACCGTGCAGGCGATGGCAGGCAAATTGCACTGGACGGGGTACATCCTGCAAATTGGTTTCTCCTTTGCGATGCCGACCTtaatc                                   >  2:242619/1‑235 (MQ=255)
                                                                                                                                                          aaaGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGACACCGTGCAGGCGATGGCAGGCAAAGTGCACTGGACGGTGTACATCCTGCAAACTGGTTTCTCCTTTGCGGTGGCGATCTtcatcat                                 >  2:452110/1‑236 (MQ=255)
                                                                                                                                                               cgcgCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGACACCGTGCAGGCGATGGCAGGCAAAGTGCACTGGACGGTGTACATCCTGCAAACTGGTTTCTCCTTTTCGGTGGCGATCTtcatcatcaacc                            >  2:597187/1‑233 (MQ=255)
                                                                                                                                                               cgcgCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGACACCGTGCAGGCGATGGCAGGCAAAGTGCACTGGACGGTGTACATCCTGCAAACTGGTTTCTCCTTTGCGGTGGCGATCTTCATCATCACGc                            >  1:691897/1‑236 (MQ=255)
                                                                                                                                                               cgcgCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGACACCGTGCAGGCGATGGCAGGCAAAGTGCACTGGACGGTGTACATCCTGCAAACTGGTTTCTCCTTTGCGGTGGCGATCTTCATCATCACGc                            >  1:111843/1‑236 (MQ=255)
                                                                                                                                                                 cgcCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGtga                                                                                                                                                                   >  1:466553/1‑99 (MQ=255)
                                                                                                                                                                 cgcCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGtga                                                                                                                                                                   <  2:466553/99‑1 (MQ=255)
                                                                                                                                                                  gcCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGACACCGTGCAGGCGATGGCAGGCAAAGTGCACTGGACGGTGTACATCCTGCAAACTGGTTTCTCCTTTGCGGTGGCGATCTTCATCATCACg                             <  2:709624/232‑1 (MQ=255)
                                                                                                                                                                  gcCGTTCCGGGGCAAAAAGGTGGAGAGCTTTAAAGGCCTCAAATTGCGAGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGTTATCGACACCGTGCAGGCGATGGCAGGCAAAGTGCACTGGACGGTGTACATCCTGCAAACTGGTTTCTCCTTTGCGGTGGCGATCTTCATCATCACGCa                           <  2:560135/234‑1 (MQ=255)
                                                                                                                                                                  gcCGTTCCGCGGCAAAAAAGATGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGACACCGTGCAGGCGATGGCAGGCAAAGTGCACTGGACGGTGTACATCCTGCAAACTGGTTTCTCCTTTGCGGTGGCGATCTTCATCATCACGCa                           <  2:199674/234‑1 (MQ=255)
                                                                                                                                                                  gcCGTTCATGGGCAAAAAAGAGGCGGACGTTGAAGCCCTCAAATTGCCGGGCTGGCTGCCCATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGACACCGTGCAGGCGATGGCAGGCAAAGTGCACTGGACGGTGTACATCCTGCAAACTGGTTTCTCCTTTGCGGTGGCGATCTTCATCATCACg                             <  2:76791/232‑1 (MQ=255)
                                                                                                                                                                   ccGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGACACCGTGCAGGCGATGGCAGGCAAAGTGCACTGGACGGTGTACATCCTGCAAACTGGTTTCTCCTTTGCGGTGGCGATCTTCATCATCACGCAggg                        >  2:654664/1‑236 (MQ=255)
                                                                                                                                                                         cTCGACAAAAAAGGGGCGAGCGTTGAAGGCCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGACACCGTGCAGGCGATGGCAGGCAAAGTGCACTGGACGGTGTACATCCTGCAAACTGGTTTCTCCTTTGCGGTGGCGATCTTCATCATCACGCAGGGTGTgcg                  <  2:351524/236‑1 (MQ=255)
                                                                                                                                                                                        gagaGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGACACCGTGCAGGCGATGGCAGGCAAAGTGCACTGGACGGTGTACATCCTGCAAACTGGTTTCTCCTTTGCGGTGGCGATCTTCATCATCACGCAGGGTGTGCGCATGTTTGTGGCg     >  1:508490/1‑234 (MQ=255)
                                                                                                                                                                                             cGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGACACCGTGCAGGCGATGGCAGGCAAAGTGCACTgg                                                                                      >  1:536528/1‑148 (MQ=255)
                                                                                                                                                                                             cGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGACACCGTGCAGGCGATGGCAGGCAAAGTGCACTgg                                                                                      <  2:536528/148‑1 (MQ=255)
                                                                                                                                                                                             cGTTGAAGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGACACCGTGCAGGCGATGGCAGGCAAAGTGCACTGGACGGTGTACATCCTGCAAACTGGTTTCTCCTTTGCGGTGGCGATCTTCATCATCACGCAGGGTGTGCGCATGTTTGTGGCGGaa  >  1:71107/1‑232 (MQ=255)
                                                                                                                                                                                                  aaGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGACACCGTGCAGGCGATGGCAGGCAAAGTGCACTGGACGGTGTACATCCTGCAAACTGGTTTCTCCTTTGCGGTGGCGATCTTCATCATCACGCa                           >  1:569781/1‑202 (MQ=255)
                                                                                                                                                                                                  aaGACCTCAAATTGCCGGGCTGGCTGCACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGACACCGTGCAGGCGATGGCAGGCAAAGTGCACTGGACGGTGTACATCCTGCAAACTGGTTTCTCCTTTGCGGTGGCGATCTTCATCATCACGCa                           <  2:569781/202‑1 (MQ=255)
                                                                                                                                                                                                                            |                                                                                                                                                                                                        
CTCCATGTGGACCACCATTATCTGTACCGCGATTCTGGTTTCGCTCTACTGGGGCATCACTTCCAACATGATGTACAAGCCGACTCAGGAAGTGACGGATGGCTGTGGTTTCTCCATCGGTCACCAGCAGCAGTTTGCATCATGGATTGCCTATAAAGTCGCGCCGTTCCTCGGCAAAAAAGAGGAGAGCGTTGAAGACCTCAAATTGCCGGGCTGGCTGAACATTTTCCACGACAACATCGTCTCCACGGCGATTGTGATGACCATCTTCTTTGGTGCCATTCTGCTCTCCTTCGGTATCGACACCGTGCAGGCGATGGCAGGCAAAGTGCACTGGACGGTGTACATCCTGCAAACTGGTTTCTCCTTTGCGGTGGCGATCTTCATCATCACGCAGGGTGTGCGCATGTTTGTGGCGGAA  >  NC_000913/4420438‑4420858

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