Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,285,409 A→G S40S (TCT→TCC)  unknown ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,285,4090AG100.0% 129.5 / NA 44S40S (TCT→TCC) unknownhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/18);  total (26/18)

GCTAATTACGTAAGAGCGGTAATTATATAAAACGAATTTAAATAAACAGAAAATTATATGAAATGGGGAAGTATCCCACTCGATTATTTATCGGCGATATAATGCAATGATTTATAAATGAATAAAAATGATATTGTTTATTCTCTACCGTCATTTTTACCGTCACCAGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGGCCCATAGACGCTTTAACACTTTTGGGGAATTTACCTTGTTTAATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGCCCATTTCTGGGCTTGTTATCTAACCAATATCGCGATTTAATATTTCGCGGATCTTGTCGCCTGATTTCTGCTCAATGAAATCATTTATCACTTTATACA  >  1/1285168‑1285653
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
gCTAATTACGTAAGAGCGGTAATTATATAAAACGAATTTAAATAAACAGAAAATTATATGAAATGGGGAAGTATCCCACTCGATTATTTATCGGCGACATAATGCAATGATTTATCAATGAATAAAATTGATATTGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTa                                                                                                                                                                                                                                             >  1:267202/1‑251 (MQ=255)
                                                   aaTTATATGAAATGGGGAAGTATCCCACTCGATTATTTATCGGCGACATAATGCAATGATTTATCAATGAATAAAATTGATATTGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCg                                                                                                                                                                                          >  1:353600/1‑251 (MQ=255)
                                                                   gAAGTATCCCACTCGATTATTTATCGGCGACATAATGCAATGATTTATCAATGAATAAAATTGATATTGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACtt                                                                                                                                                                          >  1:534059/1‑251 (MQ=255)
                                                                    aaGTATCCCACTCGATTATTTATCGGCGACATAATGCAATGATTTATCAATGAATAAAATTGATATTGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTcatcat                                                                                                                                                                                                             >  1:573166/1‑215 (MQ=255)
                                                                                cGATTATTTATCGGCGACATAATGCAATGATTTATCAATGAATAAAATTGATATTGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTg                                                                                                                                                              >  1:811968/1‑250 (MQ=255)
                                                                                     atttatCGGCGACATAATGCAATGATTTATCAATGAATAAAATTGATATTGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAg                                                                                                                                                                                    <  1:186895/223‑1 (MQ=255)
                                                                                     atttatCGGCGACATAATGCAATGATTTATCAATGAATAAAATTGATATTGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTa                                                                                                                                                         >  1:113098/1‑250 (MQ=255)
                                                                                           cGGCGACATAATGCAATGATTTATCAATGAATAAAATTGATATTGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTtagta                                                                                                                                                  >  1:787709/1‑251 (MQ=255)
                                                                                                                       gAATAAAATTGATATTGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGtctcatca                                                                                                                                                                                                              <  1:558209/163‑8 (MQ=255)
                                                                                                                        aaTAAAATTGATATTGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTtctc                                                                                                                     <  1:791741/251‑1 (MQ=255)
                                                                                                                               ttGATATTGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAAttt                                                                                                                                                                                                                            >  1:787692/2‑141 (MQ=255)
                                                                                                                                      tGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACGCGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAAtt                                                                                                                                                                                                   <  1:844659/159‑1 (MQ=255)
                                                                                                                                      tGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTAcc                                                                                                                                                                                                                         >  1:374983/1‑137 (MQ=255)
                                                                                                                                      tGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCtt                                                                                                                                                           <  1:78960/199‑1 (MQ=255)
                                                                                                                                       gTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTtagtag                                                                                                                                                 <  1:208337/208‑1 (MQ=255)
                                                                                                                                       gTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAg                                                                                                      >  1:825922/1‑251 (MQ=255)
                                                                                                                                           aTTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGTTTTTATCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGCCCa                                                                                                  >  1:412077/1‑251 (MQ=255)
                                                                                                                                           aTTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCtt                                                                                                                                                           <  1:730566/194‑1 (MQ=255)
                                                                                                                                           aTTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGCCCa                                                                                                  >  1:340052/1‑251 (MQ=255)
                                                                                                                                           aTTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCAATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCt                                                                                                                                                                                  >  1:361858/1‑171 (MQ=255)
                                                                                                                                            ttCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGCCCAt                                                                                                 >  1:488936/1‑251 (MQ=255)
                                                                                                                                                   ccGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTgggg                                                                                                                                                                                                                                 <  1:468394/116‑1 (MQ=255)
                                                                                                                                                    cGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGCCCATTTCTGGGc                                                                                         <  1:412717/251‑1 (MQ=255)
                                                                                                                                                         tttttaccgtctttggaagaaAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGCCCATTTCTGGGCTTGtt                                                                                    >  1:640421/16‑251 (MQ=255)
                                                                                                                                                         tttttaccgtctttggaagaaAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGCCCATTTCTGGGCTTGtt                                                                                    >  1:142387/16‑251 (MQ=255)
                                                                                                                                                                  tctttggaagaaAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTGCCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTGCTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGCCCATTTCTGGGCTTGTTATCTAACCa                                                                           <  1:249440/245‑1 (MQ=255)
                                                                                                                                                                                   cATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGCCCATTTCTGGGCTTGTTATCTAACCAATATCGCGATTTAatat                                                          >  1:336520/1‑251 (MQ=255)
                                                                                                                                                                                    aTCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACtt                                                                                                                                                      <  1:884412/158‑1 (MQ=255)
                                                                                                                                                                                    aTCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGCCCATTTCTGGGCTTGTTATCTAACCAATATCGCGATTTAATAtt                                                         >  1:287794/1‑251 (MQ=255)
                                                                                                                                                                                      cAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAAttgttg                                                                                                                                                                                               >  1:683844/1‑115 (MQ=255)
                                                                                                                                                                                         gAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGCCCATTTCTGGGCTTGTTATCTAACCAATATCGCGATTTAATAttt                                                        <  1:15111/247‑1 (MQ=255)
                                                                                                                                                                                                    aaGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGCCCATTTCTGGGCTTGTTATCTAACCAATATCGCGATTTAATATTTCGCGGATCTTGTCGc                                         <  1:417985/251‑1 (MQ=255)
                                                                                                                                                                                                      gCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAAt                                                                                                                                                                                                    <  1:846544/94‑1 (MQ=255)
                                                                                                                                                                                                           gTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGCCCATTTc                                                                                              <  1:397612/191‑1 (MQ=255)
                                                                                                                                                                                                                 atatCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGCCCATTTCTGGGCTTGTTATCTAACCAATATCACGATTTAATATTTCGCGGATCTTGTCGCCTGATTTCTGCt                             >  1:246589/1‑250 (MQ=255)
                                                                                                                                                                                                                                    cTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGCCCATTTCTGGGCTTGTTATCTAACCAATATCGCGATTTAATATTTCGCGGATCTTGTCGCCTGATTTCTGCTc                            <  1:671454/232‑1 (MQ=255)
                                                                                                                                                                                                                                    cTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTG‑‑‑‑‑CTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGCCCATTTCTGGGCTTGTTATCTAACCAATATCGCGATTTAATATTTCGCGGATCTTGTCGCCTGATTTCTGCTCAATGAAATCATTTATCACTTtata    >  1:68068/1‑251 (MQ=37)
                                                                                                                                                                                                                                         cgccccATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGCCCATTTCTGGGCTTGTTATCTAACCAATATCGCGATTTAATATTTCGCGGATCTTGTCGCCTGATTTCTGCTCAATGAAATCATTTATCACTTtata    >  1:387462/4‑251 (MQ=255)
                                                                                                                                                                                                                                           ccccATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGcc                                                                                                    <  1:832729/152‑1 (MQ=255)
                                                                                                                                                                                                                                           ccccATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGCCCATTTCTGGGCTTGTTATCTAACCAATATCGCGATTTAATATTTCGCGGATCTTGTCGCCTGATTTCTGCTCAATGAAATCATTTATCACTTTATAc   >  1:82326/2‑250 (MQ=255)
                                                                                                                                                                                                                                           ccccATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGCCCATTTCTGGGCTTGTTATCTAACCAATATCGCGATTTAATATTTCGCGGATCTTGTCGCCTGATTTCTGCTCAATGAAATCATTTATCACTTTATAc   <  1:67171/249‑1 (MQ=255)
                                                                                                                                                                                                                                           ccccATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGCCCATTTCTGGGCTTGTTATCTAACCAATATCGCGATTTAATATTTCGCGGATCTTGTCGCCTGATTTCTGCTCAATGAAATCATTTATCACTTTATACa  >  1:593200/2‑251 (MQ=255)
                                                                                                                                                                                                                                           ccccATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGCCCATTTCTGGGCTTGTTATCTAACCAATATCGCGATTTAATATTTCGCGGATCTTGTCGCCTGATTTCTGCTCAATGAAATCATTTATCACTTTATACa  >  1:525486/2‑251 (MQ=255)
                                                                                                                                                                                                                                             ccATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGCCCATTTCTGGGCTTGTTATCTAACCAATATCGCGATTTAATATTTCGCGGATCTTGTCGCCTGATTTCTGCTCAATGAAATCATTTATCACTTTATAc   >  1:661920/1‑248 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
GCTAATTACGTAAGAGCGGTAATTATATAAAACGAATTTAAATAAACAGAAAATTATATGAAATGGGGAAGTATCCCACTCGATTATTTATCGGCGATATAATGCAATGATTTATAAATGAATAAAAATGATATTGTTTATTCTCTACCGTCATTTTTACCGTCACCAGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGGCCCATAGACGCTTTAACACTTTTGGGGAATTTACCTTGTTTAATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGCCCATTTCTGGGCTTGTTATCTAACCAATATCGCGATTTAATATTTCGCGGATCTTGTCGCCTGATTTCTGCTCAATGAAATCATTTATCACTTTATACA  >  1/1285168‑1285653

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