Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,078,229 C→T P261S (CCG→TCG)  fabY → putative fatty acid biosynthesis enzyme FabY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,078,2290CT98.8% 270.9 / NA 81P261S (CCG→TCG) fabYputative fatty acid biosynthesis enzyme FabY
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (33/47);  minor base G (1/0);  total (34/47)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

CAGGCGTGGTACGAACATATCCCGCTTAGTGAAAAAATGGGCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGCCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAGCTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGG  >  NC_000913/4077983‑4078475
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
cAGGCGTGGTACGAACATATCCCGCTTAGTGAAAAAATGGGCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGAt                                                                                                                                                                                                                                                    >  1:290767/1‑251 (MQ=255)
    cGTGGTACGAACATATCCCGCTTAGTGAAAAAATGGGCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGc                                                                                                                                                                                                                                                >  2:128876/1‑251 (MQ=255)
     gTGGTACGAACATATCCCGCTTAGTGAAAAAATGGGCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCg                                                                                                                                                                                                                                               <  1:419254/251‑1 (MQ=255)
           cGAACATATCCCGCTTAGTGAAAAAATGGGCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAAc                                                                                                                                                                                                                                         <  1:596165/251‑1 (MQ=255)
           cGAACATATCCCGCTTAGTGAAAAAATGGGCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAAc                                                                                                                                                                                                                                         <  2:179837/251‑1 (MQ=255)
           cGAACATATCCCGCTTAGTGAAAAAATGGGCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAAc                                                                                                                                                                                                                                         <  2:247490/251‑1 (MQ=255)
            gAACATATCCCGCTTAGTGAAAAAATGGGCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGCTGCTGCGTGCACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAAcc                                                                                                                                                                                                                                        >  2:258185/1‑251 (MQ=255)
                   tCCCGCTTAGTGAAAAAATGGGCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCg                                                                                                                                                                                                                                 >  2:449671/1‑251 (MQ=255)
                        cTTAGTGAAAAAATGGGCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGc                                                                                                                                                                                                                            >  1:190824/1‑251 (MQ=255)
                        cTTAGTGAAAAAATGGGCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATAACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGc                                                                                                                                                                                                                            >  2:87738/1‑251 (MQ=255)
                              gAAAAAATGGGCGTGCGCATTCAGCAATGTACCGGGCAAAAATTTATCAGTACCAATCCATAAACCGGCAATCAGAAACCGCACCATACGCGGTTTGCCGGGAGGTTTTTCTCACTGGCGACGCTGACTGGTTGGGGACTTATCTGGCTGATGCTTCGTGAACGGCACCTCTGAGGAATGATTATTCTTACGGGTTCGCATATCGGCTTCTGGGAGTCGATTAGCGGTTAACCTCAGGCGGTactggacct                                                                                                                                                                                                                      <  2:269848/251‑9 (MQ=255)
                               aaaaaaTGGGCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTc                                                                                                                                                                                                                     >  2:149119/1‑251 (MQ=255)
                                 aaaaTGGGCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCgg                                                                                                                                                                                                                   >  1:42648/1‑251 (MQ=255)
                                            gcgcATTCAACAATATACCGGGCAAAAATTTCTCAACACCAAGCCAGAAACCGGCAATCAGAATCCGCACCATTCGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTACCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAgc                                                                                                                                                                                                         <  2:39028/250‑1 (MQ=255)
                                                        atatACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGAcc                                                                                                                                                                                            <  1:258185/251‑1 (MQ=255)
                                                            aCCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCt                                                                                                                                                                                        <  1:449671/251‑1 (MQ=255)
                                                               gggCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTAGCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAGCCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGgc                                                                                                                                                                                     <  1:606604/251‑1 (MQ=255)
                                                                   aaaaaTTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGgcgcgc                                                                                                                                                                                 >  1:24163/1‑251 (MQ=255)
                                                                       agttATCACTACCATGCCAGAAACCGGCAATCATAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTTCAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGAc                                                                                                                                                                             <  2:643703/249‑1 (MQ=255)
                                                                        tttATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACg                                                                                                                                                                            >  1:478792/1‑251 (MQ=255)
                                                                           aTCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGaaa                                                                                                                                                                         <  1:31416/251‑1 (MQ=255)
                                                                           aTCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGaa                                                                                                                                                                          >  1:225153/1‑250 (MQ=255)
                                                                                    aTGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAGGCACGGGt                                                                                                                                                                <  1:134184/251‑1 (MQ=255)
                                                                                        cAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGAc                                                                                                                                                                             >  2:148644/1‑234 (MQ=255)
                                                                                        cAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGAc                                                                                                                                                                             <  1:148644/234‑1 (MQ=255)
                                                                                        cAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTAGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGGCTGGCGCGCGGACGAAAAGCACGGGTGc                                                                                                                                                              >  1:296271/1‑249 (MQ=255)
                                                                                         aGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGc                                                                                                                                                              >  1:325005/1‑248 (MQ=255)
                                                                                           aaaCCGGCAATCAGAATCCGAACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATg                                                                                                                                                         <  2:45141/251‑1 (MQ=255)
                                                                                                ggCAATCAGAATCCGCACCATACGCTGTTTGCAGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAAACTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGCTCTCGCCCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGt                                                                                                                                                    >  2:504223/1‑251 (MQ=255)
                                                                                                      cAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAAt                                                                                                                                              <  1:188599/251‑1 (MQ=255)
                                                                                                                 ccATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGcgacg                                                                                                                                   >  2:366943/1‑251 (MQ=255)
                                                                                                                    tACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGcgacg                                                                                                                                   >  1:123392/1‑248 (MQ=255)
                                                                                                                     aCGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCCCCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACgaga                                                                                                                                <  1:639318/250‑1 (MQ=255)
                                                                                                                          gTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCgg                                                                                                                          >  2:198721/1‑251 (MQ=255)
                                                                                                                             tGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGTCCTCGGTGCCTTAAGCGGCGAtc                                                                                                                                                                                                  <  2:99781/176‑1 (MQ=255)
                                                                                                                             tGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGAtc                                                                                                                                                                                                  >  1:99781/1‑176 (MQ=255)
                                                                                                                              gCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGGTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGc                                                                                                                      <  2:177463/251‑1 (MQ=255)
                                                                                                                                  ggAGTTTATTCTCACTGGCGACGCTCACTGGTTTGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGCCGAGACGCCGGGTGCAgtg                                                                                                                  >  1:312574/1‑251 (MQ=255)
                                                                                                                                     gTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGttt                                                                                                               <  1:661738/251‑1 (MQ=255)
                                                                                                                                         aTTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAg                                                                                                           <  2:407354/251‑1 (MQ=255)
                                                                                                                                         aTTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAg                                                                                                           <  1:576925/251‑1 (MQ=255)
                                                                                                                                            ctcACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGACCGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGAAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGACGGcc                                                                                                        >  2:63803/1‑250 (MQ=255)
                                                                                                                                                    cGACGCTCACTGGTTGGGGTCTTATCTGGCTGATGCTACGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATc                                                                                                <  2:405676/251‑1 (MQ=255)
                                                                                                                                                        gCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCa                                                                                                        <  2:105805/239‑1 (MQ=255)
                                                                                                                                                        gCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCa                                                                                                        >  1:105805/1‑239 (MQ=255)
                                                                                                                                                          tCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTg                                                                                          >  2:265557/1‑251 (MQ=255)
                                                                                                                                                          gccctgggttgggACTTATCGGGGTGATGGTGCGTGAACGCGAACGCGGCGGATCGGATAGTCGTGCCGATGCGCATAGCCGCTACAGCGAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCTGTGCCTTAAGCGGCGATCTCGACCGTCTTGCGCGCCTACGGAATGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGGCGAGACGCCGGGTGCAGTGTTTGAAGGcgagtttagcgttctg                                                                                          <  2:413320/241‑16 (MQ=255)
                                                                                                                                                               ggTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGccc                                                                                       <  1:481575/249‑1 (MQ=255)
                                                                                                                                                               ggTTGGGGAATTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCcg                                                                                      <  2:366393/250‑1 (MQ=255)
                                                                                                                                                                    gggACTTATCTGGCTGATGCTGCGTGAACGCCACCGCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGc                                                                                <  1:504223/251‑1 (MQ=255)
                                                                                                                                                                     ggACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCATATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGcc                                                                               <  1:367603/251‑1 (MQ=255)
                                                                                                                                                                     ggACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGGCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGcc                                                                               >  2:395865/1‑251 (MQ=255)
                                                                                                                                                                       aCTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGCTTATTCTTGCTGCTGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGCCCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCAt                                                                             >  2:64601/1‑251 (MQ=255)
                                                                                                                                                                          tttCTGGCTGATGCTGCGTGTACGCCTCCTCGGCTGAACGATTATTCTTGCGGTGGCGCATATCCGCTACATCAAGTCGATTAGCGTTAATCCTCATGCGGTAGCCTACCTCGGTGCCGTAAGCGGCGATCTCGACCGTCTGTCGCGCTGACTAAAATCACGGGTGCAGATGCAGGTCGAAATCTTGGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCGGCCCGCGAAGCCAttt                                                                           <  2:372526/248‑1 (MQ=255)
                                                                                                                                                                               ggCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTgg                                                                         <  1:43114/247‑1 (MQ=255)
                                                                                                                                                                               ggCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTgg                                                                         <  2:354179/247‑1 (MQ=255)
                                                                                                                                                                               ggCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGc                                                                        <  2:22560/248‑1 (MQ=255)
                                                                                                                                                                                  tGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTAt                                                                  >  2:218581/1‑251 (MQ=255)
                                                                                                                                                                                     tgctgcGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAgcg                                                                                                                                                                                                        >  1:110625/1‑114 (MQ=255)
                                                                                                                                                                                     tgctgcGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAgcg                                                                                                                                                                                                        <  2:110625/114‑1 (MQ=255)
                                                                                                                                                                                      gctgcGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCATAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGa                                                                >  2:516342/1‑249 (MQ=255)
                                                                                                                                                                                         gcGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAaga                                                             <  1:218581/249‑1 (MQ=255)
                                                                                                                                                                                           gTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGc                                                         <  2:186130/251‑1 (MQ=255)
                                                                                                                                                                                           gTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGc                                                         >  1:323784/1‑251 (MQ=255)
                                                                                                                                                                                             gAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCTTATGAAGAGGGCgg                                                       <  1:66827/251‑1 (MQ=255)
                                                                                                                                                                                             gAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCgg                                                       <  1:555601/251‑1 (MQ=255)
                                                                                                                                                                                                cGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGaa                                                    <  1:96029/251‑1 (MQ=255)
                                                                                                                                                                                                 gCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAAc                                                   <  1:185111/251‑1 (MQ=255)
                                                                                                                                                                                                             gAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCAGGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAg                                       <  1:510804/251‑1 (MQ=255)
                                                                                                                                                                                                               aCGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAGCt                                     <  1:396346/251‑1 (MQ=255)
                                                                                                                                                                                                               aCGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCAGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAGCt                                     <  2:25078/251‑1 (MQ=255)
                                                                                                                                                                                                                 gATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAGCTg                                    >  1:605531/1‑250 (MQ=255)
                                                                                                                                                                                                                  attattCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGAGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAGCTGAt                                  >  1:136701/1‑251 (MQ=255)
                                                                                                                                                                                                                      ttCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCg                                                                                               <  1:265442/186‑1 (MQ=255)
                                                                                                                                                                                                                      ttCTTGCGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCg                                                                                               >  2:265442/1‑186 (MQ=255)
                                                                                                                                                                                                                            cGGATGCGCATATCCGCTACAGCAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAGCTGATGGAGAGGGGc                        <  1:312846/251‑1 (MQ=255)
                                                                                                                                                                                                                            cGGATGCGCATATCAGAGACAGCAAGTCGATTAGCGGTAAACCTCATGCGGAAGCCGACCTCGGCGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAACAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACACCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAACCCCTTCGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAGCGGGTGGGGAGGGGc                        >  2:455726/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gcgcATATCCGCTACAGCACGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGAAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGACAGCGAAGCCATTTGGCCCGTATGCAGAGGGCGGGAACGAAGAAGAGTAGCTGATGGAGAGGGGCGGtgc                   >  2:71872/1‑251 (MQ=255)
                                                                                                                                                                                                                                     atatCCGCTACAGCAAGTCGAGTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCGTTGGCGACGAGACGCCGTGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATGTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAGCTGATGGAGAGGGGCGTTGCTGcc               <  2:225153/251‑1 (MQ=255)
                                                                                                                                                                                                                                              aCAGCAAGTCGATTATCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAGCTGATGGAGAGGGGCGGTGCTGCCTCTCTCAtt      <  1:71872/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  cAAGTCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAGCTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAgg  <  2:42648/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
CAGGCGTGGTACGAACATATCCCGCTTAGTGAAAAAATGGGCGTGCGCATTCAGCAATATACCGGGCAAAAATTTATCACTACCATGCCAGAAACCGGCAATCAGAATCCGCACCATACGCTGTTTGCCGGGAGTTTATTCTCACTGGCGACGCTCACTGGTTGGGGACTTATCTGGCTGATGCTGCGTGAACGCCACCTCGGCGGAACGATTATTCTTGCGGATGCGCATATCCGCTACAGCAAGCCGATTAGCGGTAAACCTCATGCGGTAGCCGACCTCGGTGCCTTAAGCGGCGATCTCGACCGTCTGGCGCGCGGACGAAAAGCACGGGTGCAGATGCAGGTCGAAATCTTTGGCGACGAGACGCCGGGTGCAGTGTTTGAAGGCACGTATATCGTTCTGCCCGCGAAGCCATTTGGCCCGTATGAAGAGGGCGGGAACGAAGAAGAGTAGCTGATGGAGAGGGGCGGTGCTGCCTCTCTCATTCAGG  >  NC_000913/4077983‑4078475

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