Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,756,976 G→A R273Q (CGA→CAA)  intA → integrase IntA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,756,9760GA96.2% 255.8 / ‑0.1 80R273Q (CGA→CAA) intAintegrase IntA
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (41/36);  total (43/37)
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.53e-01

ATTGATGCTGTTCCTAGTGTAAATATCGGAAAAGCTTTCGAGAAACCGCAAAAGAAAAACATGCCAAGCATCCGGCCGGATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTACGGCACTCAATGAGG  >  NC_000913/2756729‑2757206
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                      
attgatGCTGTTCCTAGTGTAAATATCGGAAAAGCTTTCGAGAAACCGCAAAAGAAAAACATGCCAAGCATCCGGCCGGATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAAt                                                                                                                                                                                                                                     <  1:739356/251‑1 (MQ=255)
attgatGCTGTTCCTAGTGTAAAGATAGGAAAAGCTTTCGAGAAACCGCAAAAGAAAAAAATGCCAAGCATCCGGCCGGATCAACTTCCGCAGCTAATGAACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCAGGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAAt                                                                                                                                                                                                                                     <  2:189028/251‑1 (MQ=255)
          ttCCTAGTGTAAATATCGGAAAAGCTTTCGAGAAACCGCAAAAGAAAAACATGCCAAGCATCCGGCCGGATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAAc                                                                                                                                                                                                                           >  1:688388/1‑251 (MQ=255)
            cctcgtgtcactaccGGAAAAGCTTTCGAGAAACCGCAAAAGAAAATCACTCCAAGCATCCGTCCGGTTCAACTTCCGCAGCTAATTCACACCATGCTTACGTCAAGTATCAGCACGGCCACAAGATGCCTGTTCATGTGGCAACTACTAACCATCACCCTCCCTGCCGAAGCTGCTGAGGCTCTATGGGATGAGAACGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAAcc                                                                                                                                                                                                                          <  2:32412/236‑1 (MQ=255)
            ccgagtttAAAAATCGGAAAAGCTTGCGAGAAACCGCAAAAGAAAAACAGGCCAACCACCCGGCCGGATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGTCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCg                                                                                                                                                                                                                         <  2:825034/244‑1 (MQ=255)
                  gtAAATATCGGAAAAGCTTTCGAGAAACCGCAAAAGAAAAACATGCCAAGCATCCGGCCGGATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGAcc                                                                                                                                                                                                                    <  1:550140/250‑1 (MQ=255)
                     aaTATCGGAAAAGCTTTCGAGAAACCGCAAAAGAAAAACATGCCAAGCATCCGGCCGGACCAACTTCCGCAGCTAATGCACCCCAGGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATaa                                                                                                                                                                                                                <  2:310118/251‑2 (MQ=255)
                        atCGGAAAAGCTTTCGAGAAACCGCAAAAGAAAAACATGCCAAGCATCCGGCCGGATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGt                                                                                                                                                                                                             >  1:869672/1‑251 (MQ=255)
                          cGGAAAAGCTTTCGAGAAACCGCAAAAGAAAAACATGCCAAGCATCCGGCCGGATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGtt                                                                                                                                                                                                            <  1:249523/250‑1 (MQ=255)
                            gAAAAGCTTTCGAGAAACCGCAAAAGAAAAACATGCCAAGCATCCGGCCGGATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTccc                                                                                                                                                                                                         >  2:704775/1‑250 (MQ=255)
                                    ttCGAGAAACCGCAAAAGAAAAACATGCCAAGCATCCGGCCGGATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTg                                                                                                                                                                                                  <  1:198150/250‑1 (MQ=255)
                                    ttCGAGAAACCGCAAAAGAAAAACATGCCAAGCATCCGGCCGGAACAACTTCCGCAGATAATGCACACCAGTCGTACGGCAAGTACCAGCAGGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTg                                                                                                                                                                                                  <  2:692202/250‑1 (MQ=255)
                                    gtCGAGAAACCGAAAAAGAAAAACATGCCAAGCACCCGGCCGGATCAACTTCCGCAGCTAATGCTCACCATGCGTACGTCAATGATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCAAACCATCCCCCGCCCTGCCGAAGCTGCTTAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTTCAGCTCAAATGAAGATGACCCGGGACCATACGGTTCCACTATCTg                                                                                                                                                                                                  <  2:692191/249‑1 (MQ=255)
                                      cGAGAAACCGCAAAAGACAAACATGCCAAGCATCCGGCCGGATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGATCAAATGAAGATGAACAGGGAcccaccggttcccttatctgctg                                                                                                                                                                                               >  2:241696/1‑230 (MQ=255)
                                          aaaCCGCAAAAGAAAAACATGCCAAGCATCCGGCCGGATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCACTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGACTGGCAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATATGATGGGGc                                                                                                                                                                                           >  1:77276/1‑251 (MQ=255)
                                                cAAAAGAAAAACATGCCAAGCATCCGGCCGGATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGc                                                                                                                                                                                     >  2:531732/1‑251 (MQ=255)
                                                   aaGAAAAACATGCCAAGCATCCGGCCGGATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGACTGGAAAATTCATGAAGCTCAAATGAAGATGAACCGGGACCCTACGGTTACACTATCTGATGGGGCTCTTGCTAt                                                                                                                                                                                  >  2:724656/1‑251 (MQ=255)
                                                    aGAAAAACATGCCAAGCATCCGGCCGGATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTAtt                                                                                                                                                                                 >  2:675395/1‑251 (MQ=255)
                                                     gAAAAACATGCCAAGCATCCGGCCGGATCAACTTACGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTc                                                                                                                                                                                >  1:278900/1‑251 (MQ=255)
                                                       aaaaCATGCCAAGCATCCGGCCGGATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTg                                                                                                                                                                              >  1:382312/1‑251 (MQ=255)
                                                          aCATGCCAAGCATCCGGCCGGATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGa                                                                                                                                                                            >  1:225459/1‑250 (MQ=255)
                                                              gCCAAGCATCCGGCCGGATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAatg                                                                                                                                                                        >  1:458173/1‑250 (MQ=255)
                                                                 aagaaTCCGGCCGGATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAAtgatga                                                                                                                                                                    <  2:209390/247‑1 (MQ=255)
                                                                 aaGCATCCGGCCGGATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAatgat                                                                                                                                                                      >  1:706764/1‑249 (MQ=255)
                                                                 aaGCATCCGGCCGGATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACGTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGTACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATGCTGGAAatgat                                                                                                                                                                      <  2:706764/249‑1 (MQ=255)
                                                                   gCATCCGGCCGGATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCCCCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGTACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAg                                                                                                                                                                  <  1:391623/251‑1 (MQ=255)
                                                                   gCATCCGGCCGGATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAAtgatga                                                                                                                                                                    >  1:635933/1‑249 (MQ=255)
                                                                    cacccggccggATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTGCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGc                                                                                                                                                                 <  2:85845/248‑1 (MQ=255)
                                                                       ccggccggATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAACTGAAGATGAACCGGGACCATACGGTTCCACTATCTGCTGGGGCTCTTGCTATTCTGGAAATGATGACGCctc                                                                                                                                                              >  2:569093/1‑251 (MQ=255)
                                                                               aTCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCTCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAgtggt                                                                                                                                                      >  1:506352/1‑251 (MQ=255)
                                                                                       ccGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCgaga                                                                                                                                              <  1:450856/251‑1 (MQ=255)
                                                                                            gctaagtcacaCCATGCGTCCGGCAAGTATCAGCATGTCCACCAGATGCCTGTTCATGTGGCAGCTTCTAACCAGCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCATTGGTGGCCGAGAATTTa                                                                                                                                         <  2:278424/244‑1 (MQ=255)
                                                                                                      acCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCt                                                                                                                                 >  1:724724/1‑249 (MQ=255)
                                                                                                         aTGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCg                                                                                                                            <  1:646626/251‑1 (MQ=255)
                                                                                                            cGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATACCCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGACTCTCCGTGGTGGCCGAGAATTTATCTTTCCGAGCCGTAt                                                                                                                         >  1:422681/1‑251 (MQ=255)
                                                                                                            cGTACGACAAGTAACAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATAGATTTCAATGCTAGCGAATGGGAAATTCCTGCAGCTCAAATGACGCTACACCGGGACCATACGGGTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGCAGCCTCTCAGTGGTGGCCGCGAATTTATCTTTCCTAGCCGTAt                                                                                                                         >  2:560682/1‑251 (MQ=255)
                                                                                                             gTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATc                                                                                                                        >  1:263163/1‑251 (MQ=255)
                                                                                                                cGGCAAGGATAAGCTTGGCCACAAGAGGCCTGTTAATGTGGCAAATCCTAACCCTCACCCGCCCTGACGAAGCTGTTGAGGCTCGATGGGATGAGAGCGATGTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACAATACGGTTCCACTATCTGATGGGGATCTTGCTATTCTGGAAATGATGAAGCCTCTTAGTGGTGGCCGAGAATTTATCTTTCCGAGCCGTATCAAg                                                                                                                     <  2:605734/251‑1 (MQ=255)
                                                                                                                  gaaattATCAGCATGCCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCCCCCTCCCGGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTTCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGcc                                                                                                                   <  2:506352/246‑1 (MQ=255)
                                                                                                                  gCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTCACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGCTCGATTTCAATGCTAGCGAATGGACAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGAAAATGCTGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTCGCCGTATGAAgac                                                                                                                   >  2:430723/1‑249 (MQ=255)
                                                                                                                    aaGTATCAGCATGTCCACAAGATGACTGTTCATGTGGAAACTTCTAAACATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAACATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTATGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCACGccc                                                                                                                  >  2:871237/1‑250 (MQ=255)
                                                                                                                        aTCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGCGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCAATATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCccacccc                                                                                                             >  2:751932/1‑247 (MQ=255)
                                                                                                                            gCATGTTCACAAGATGCCTGCTCATGTGGCAACTTCTAACCCTAACCCGCCATGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTAAATGCAAGCGAATGGACAATCCCAGCAGCTCAAATGAAGATGAACCGGGACCATAAGGTACCACTATCAGAATGGGCTCTTGCTATTCTGGAAATGATGAACCCTCTCAGTGGTGGCCGAGACTTAATCTTTccgtgccatatcagtccaacccaacc                                                                                                          >  2:56517/1‑226 (MQ=255)
                                                                                                                            gCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCAACCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGGGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCccaaccaacc                                                                                                          >  2:56563/1‑250 (MQ=255)
                                                                                                                              aTGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCTATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAAt                                                                                                       <  1:518143/251‑1 (MQ=255)
                                                                                                                                      aaGATGCCTGTTCATGGGGCAACTTCTAACCATCCCCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGcc                                                                                               <  1:789227/251‑1 (MQ=255)
                                                                                                                                      aaGATCCCTGTTCATGTGGCAACTTCTAACCACCCCCCGCCCTGCCGAAGCTGCTGAGGCACGATGGGAGGAGATCGATTTCAATGCTAGCGAATGGCAAATTCCTGCAGCTCAAATTAAGATGAACCGGGACCATACGGTTCCACTATCTGAGGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGcc                                                                                               <  2:292051/251‑1 (MQ=255)
                                                                                                                                       aGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGATCCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGGGGTGGCCGAGAATTTATCTTTCCTAACCGTCTCCAGCCCAACCAACCAATGAATAGcc                                                                                               >  1:597267/1‑250 (MQ=255)
                                                                                                                                         aTGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCaa                                                                                             >  1:481211/1‑250 (MQ=255)
                                                                                                                                               gTTCATGTGGAAACGTCTAACCACCACCCGACCGGCCGAAGCTGCTGAGGCTCGATGGTATGAGATCGATTACAAGGCTAGCGAATGGAAAATTCCGTCAGCGCAAATGAAGATGAACCGGGTCCATACGGTTCCACTATCTGATGGGGCTTTTGTTATTCTGGAAATGATTAAGCCTCTCAGTGGTGTCCGAGAATTTCTCTTACCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTg                                                                                       <  2:688388/250‑1 (MQ=255)
                                                                                                                                                 tcgtgggtcacattatAACCACCCCCCGCCCTGCCGAGGCTGCTAAGGCTGTATGGGATAAGATCGATTTCAATGCTAGCGAATGGAAAATCCCTGCAGTTCAAATGAAGATGAACCGGGACAATACGGTCCCCCTATCTGACGGGGTTCGGGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCAGAGAAGTTATCTTTCCTAGACGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGaa                                                                                     <  2:713903/235‑1 (MQ=255)
                                                                                                                                                 tCATGTGGCAACTTCTAACCATCCCCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGaa                                                                                     <  1:675395/250‑1 (MQ=255)
                                                                                                                                                        gCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCGTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCAc                                                                             <  1:172876/251‑1 (MQ=255)
                                                                                                                                                         cAACTTCTAACCATCACCCGCCCTGCCGCAGCTGCTGAGGCTCGTTGGGATGAGATCGCTTTCAATGCTAGGGAATGGAAAATTCCTGCAGCTCAATTGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGTTATTATGGAAATGCTGAAGCTTCTCAGTGGTGGCTGACATTTTATCTTTCCGAGCCGTTTCAAGCCCAACAAACCAATGAATATCCAAACAGTGAACGCAGCAcc                                                                            >  2:72309/1‑250 (MQ=255)
                                                                                                                                                              ggcgttctaccgatctccggCCGAAGTTGCGGAGGATCGATGGGATGAGAGCGATTTAAATGCTAGCGAATGGAAAATTCCTTCAGATCAAATGAAGATGAACCGGGACCATACGGTTCCACGAGCTGATGGGGCTCTTGCTATTCTGGAAATGAGGAAGCCGCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCTTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGc                                                                       <  2:874126/232‑1 (MQ=255)
                                                                                                                                                               cTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCt                                                                                                                                      >  1:154470/1‑187 (MQ=255)
                                                                                                                                                               cTAACCACCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCt                                                                                                                                      <  2:154470/187‑1 (MQ=255)
                                                                                                                                                                 aaCCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGAGTg                                                                    >  1:675803/1‑251 (MQ=255)
                                                                                                                                                                   ccATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTGCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCt                                                                  <  1:388536/251‑1 (MQ=255)
                                                                                                                                                                      ccccccGCCCTGCCGAAGGGCCTGAGGCTCTAGGTGAGGAGAATGATTTAAATGCTAACGAATGGAAAATTCCTCCAGCTCAAAGGAGGATGAACCGGGACCATACGGTTCCACTATCTGAGGGGGCTCTTGTTATTCTGGAAATAATGAAGCCTCTCAGTGGTGGCCTAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGc                                                               <  2:386288/248‑1 (MQ=255)
                                                                                                                                                                       cccccGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGc                                                               <  1:776290/248‑1 (MQ=255)
                                                                                                                                                                       cACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTAt                                                                                                                                                                                                     <  2:248442/116‑1 (MQ=255)
                                                                                                                                                                       cACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTAt                                                                                                                                                                                                     >  1:248442/1‑116 (MQ=255)
                                                                                                                                                                              ccTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGACTTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGACTGCAGCACTCAAGcagatcggaagcgcacc                                                       >  1:874126/1‑235 (MQ=255)
                                                                                                                                                                                  ccGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCtt                                                                                                                                     <  2:52569/169‑1 (MQ=255)
                                                                                                                                                                                  ccGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCtt                                                                                                                                     >  1:52569/1‑169 (MQ=255)
                                                                                                                                                                                               ggCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACg                                      >  2:29878/1‑251 (MQ=255)
                                                                                                                                                                                                   ggATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGTATCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCtt                                  <  1:546783/250‑1 (MQ=255)
                                                                                                                                                                                                    gATGGGATGAGAGCGATTTCAATGATAGCGACTGAAAAATTCATGCAGCTCAAATGAAGATGAACCGGGAACATACGGTTCCACAATCTGATGGGGCTCTTGCTTTTCAGGAAATGATGAAGCCTATCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACAAACCAATGAATAGCCAAACAGGGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTATTTGTTTAACACGGATTg                                 >  2:399517/1‑251 (MQ=255)
                                                                                                                                                                                                       gggATGAGATCGATTTCAATGCTAGCGAATGGAAACTTCCTGCAGCTCAAATGAAGCTGAACCGGGACCATACGGTTCCACTATCTGCTGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGAGGTGGCCGAGAATTTCTCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGAGAATGCCGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTccaccagcttacgt                              >  2:749778/1‑238 (MQ=255)
                                                                                                                                                                                                          aTGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACAATACGGTTCCACTATCTGCTGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCACAACCAACCAATGAATAGCCAAACAGTGATTGCAGCACTCAAGCGTGCTGGTTTAGGAGGAGTACTTGTTTCACACGGCCTGCGtttt                           >  2:52609/1‑249 (MQ=255)
                                                                                                                                                                                                            gagaTCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAAGGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTa                          <  1:569093/250‑1 (MQ=255)
                                                                                                                                                                                                              aagcGATTTCAATGCTAGCGAATGGAAAATTCCTGAAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCg                       <  2:635933/248‑1 (MQ=255)
                                                                                                                                                                                                                  gATTTCAATGCTGGCGAATGGAAAATTCCTGCAGCTCAAATGCAGATGAACCGAGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAACTGATGAAGCCTCTCAGTGGTGGCCGAGAATTATTCTTTCCTAGCCGAATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGACCTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGGTCTATCGcccg                   >  2:828006/1‑249 (MQ=255)
                                                                                                                                                                                                                  gATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCAAATGCAGATGAACCGGGACAATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCAGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTAAAGCGTGCTGGCTTAGGAGGCGTACTTGTTTCACACGGCTTGCGGTGTATCGcc                     >  2:567086/1‑249 (MQ=255)
                                                                                                                                                                                                                   aTTTCACTGCTAGCGAATGGAAACTTCCTGCAGCTCAAATGAAGATGCACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGATATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGCATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCAATCAAGCGTGCTGGCTTAGGTGGAGTACTTGTTTCCCACGGCTTGCGTTCTATCGCCAgg                  >  2:349442/1‑250 (MQ=255)
                                                                                                                                                                                                                           gCTAGCGAATGGAAAATTCCTGCAACTCAACTAAAGCTGAACCGGGACCATACAGTTCCACTATCTGAAGGGGCTGTTGCTATTCTGGAAATGATGACGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCACACAGTGAATGCAGACCTCAACCGTGCTTGCTTATGAGGGACACTCGTTTCACACGGCTAGCGGTGTATCGCCAATACGTCAcc          >  2:706242/1‑250 (MQ=255)
                                                                                                                                                                                                                                 gAATGGAAAATTCCTGCAGCTCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTACGGCACTCAATGa    <  1:704775/251‑1 (MQ=255)
                                                                                                                                                                                                                                  aaTGGAAAATTCCTGCAGATCAAATGAAGCAGCACCGGGACCATACAGATCCACTATCTGATGGGGCTCTTGCTATTATGGAAATGCTGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTAACACGCCCAACCAACCAATGAATAACCAAAAAGTGAATGAAGCAATCAAGCGTGCTGACTTATGAGATGTACTTGATCCACACGCCTAGCGTTATATcccacgaacgcaactcaatggg   >  2:82016/1‑230 (MQ=255)
                                                                                                                                                                                                                                   aTGGAAAATTCCTGAAGCGCAAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGATTCACACGGCTTGAGTTCTATCGCCAGTACGGCACTCAATGagg  >  1:238754/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                      
ATTGATGCTGTTCCTAGTGTAAATATCGGAAAAGCTTTCGAGAAACCGCAAAAGAAAAACATGCCAAGCATCCGGCCGGATCAACTTCCGCAGCTAATGCACACCATGCGTACGGCAAGTATCAGCATGTCCACAAGATGCCTGTTCATGTGGCAACTTCTAACCATCACCCGCCCTGCCGAAGCTGCTGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTACGGCACTCAATGAGG  >  NC_000913/2756729‑2757206

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