Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,424,468 A→C 6.9% intergenic (+82/+48) ulaF → / ← yjfY L‑ribulose‑5‑phosphate 4‑epimerase UlaF/DUF1471 domain‑containing protein YjfY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,424,4680AC 6.9% 165.1 / 3.9 64intergenic (+82/+48)ulaF/yjfYL‑ribulose‑5‑phosphate 4‑epimerase UlaF/DUF1471 domain‑containing protein YjfY
Reads supporting (aligned to +/- strand):  ref base A (32/26);  new base C (2/3);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.69e-01

CCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATC  >  NC_000913/4424326‑4424608
                                                                                                                                               |                                                                                                                                            
ccTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGa                                                                                                                                              >  1:571174/1‑143 (MQ=255)
  tgatgaATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAg                                                                                                                                            >  2:472145/1‑143 (MQ=255)
   gatgaATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAgg                                                                                                                                           <  1:559536/143‑1 (MQ=255)
   aacgaATCATCACATGATGCGAAAAAACTGTCAAAACGCCTAATGAGGGCAGAACGAGACCACGCGAAGAGGAAATATCCTTCCTCCAGAGATAACTGGCCAAATACCGAGACCGCTTTTAAAAAAAAAGCCCCCTGCGAAgg                                                                                                                                           <  2:417761/140‑1 (MQ=255)
         tAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGgtaa                                                                                                                                     >  1:52061/1‑143 (MQ=255)
         gaaaCCCTTCATGCGTAACCCCCCTCCCAACACTTTTTGCGTTCATAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATCACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTTCGCAGGGGgctt                                                                                                                                     <  1:29689/142‑3 (MQ=255)
         gaaaCACCTCAGGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTTCGCAGGGGgctt                                                                                                                                     <  1:615667/142‑3 (MQ=255)
          aacgagttCATACGTAAACAAGGCCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCATCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTAAAAAAAAAGCCCCCTGCGAAGGGGgcaaa                                                                                                                                    <  2:278164/137‑1 (MQ=255)
          aaaCACTTGATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGgcaaa                                                                                                                                    <  1:591315/143‑1 (MQ=255)
          aaaCACTGCCTGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACAAGCGCTGAGGAAATTTACTGCCTCGGGAGATAACTGCTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGgcaaa                                                                                                                                    <  1:399779/143‑1 (MQ=255)
          aaaCACCTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGgcaaa                                                                                                                                    <  1:135358/143‑1 (MQ=255)
            acacTTCATGCGTAAACACGGTCCTAACGCGTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCCGCGACGGGGggcaagc                                                                                                                                  >  2:269008/1‑143 (MQ=255)
              acTTCATGCGTAAACCCGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTTCGCAGGGGgcttttttt                                                                                                                                <  1:178043/143‑8 (MQ=255)
               cTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCTCCTGCGAAGGGGgcaaagcaaa                                                                                                                               <  2:319655/143‑1 (MQ=255)
                ttCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAAc                                                                                                                              >  2:654765/1‑143 (MQ=255)
                   aTGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAGCTAt                                                                                                                           >  1:468439/1‑143 (MQ=255)
                    tGCGTAAGCGCGGTCCTAACGCTTATTACGTGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATg                                                                                                                          <  2:394479/143‑1 (MQ=255)
                       gTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCGCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCa                                                                                                                       <  2:720480/143‑1 (MQ=255)
                            cacGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGtt                                                                                                                  <  2:571174/143‑1 (MQ=255)
                             acGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGctt                                                                                                                 >  1:194403/1‑143 (MQ=255)
                                  ccTAACGCTTATTACAGGCAGAAGTATAACACGCGCTGCGGAAATTTCCTTCCTCGAGAGAGAACTGGTCGAATGCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTg                                                                                                            <  2:245254/143‑1 (MQ=255)
                                    tAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATACCTGGTCTAATTCCGCAGCCGTTTTTTTTGATTTACCCCCCCGCGCAGGGGGCACAGCAACCTATGGCACTATTTCGATGGt                                                                                                          >  1:166612/1‑143 (MQ=255)
                                    tAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTACTCGGGAGATAACTGGTCTAATTCCGCAGCCGATTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGt                                                                                                          >  1:256033/1‑143 (MQ=255)
                                        gCTTATTACGGGCAGAAGTAGAACACGCGATGAGGAAATTTCCTTCGTCGGGAAGTAACTGGGCTAATTCCGCAGCCGTTTTTGAAAAAAAAGCCCGCTGCGAAGGGGGCGAAGCAGACTATGGCGGTGGGTCGTTGGTtata                                                                                                      >  2:637566/1‑143 (MQ=255)
                                           tattaCGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCt                                                                                                   <  1:540738/143‑1 (MQ=255)
                                            attaGGGGAAGACGGACAAGACGCGGGGCGGCACTTTCCTTGCTCGGGAGATAACCGGTCTAATTCCGCAGCCGTTGTTCAAACAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTg                                                                                                  <  2:265990/143‑1 (MQ=255)
                                             ttaCGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTgg                                                                                                 >  2:749367/1‑143 (MQ=255)
                                                cGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAAAACAGGCAAATGATTTCAAAAAAAAACGCCGTTGTTTGTGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGc                                                                                              <  1:654765/143‑1 (MQ=255)
                                                   gCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAg                                                                                           >  2:144017/1‑143 (MQ=255)
                                                   gCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGACCAGCAACAGATAGCAATGTTTCGTTGGTTATACCTGGTGCTAg                                                                                           >  2:391722/1‑143 (MQ=255)
                                                    cAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAAGGTTTCGTTGGTTATACCTGGTGCTAGc                                                                                          >  1:383715/1‑143 (MQ=255)
                                                      gAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGa                                                                                        >  1:881388/1‑143 (MQ=255)
                                                        aGTAGAAGCCGCGCTGCGGAAATTGCGTTCTTCGAGAGATAACTGGTTTAATTCCGCAGCGGTTTTTCAAAAAAATGCTCCCTGCGAAGGTGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATa                                                                                      <  2:152687/143‑1 (MQ=255)
                                                          tAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCGCTGCGAAGGGGGCAAAGCAACCTATGGCAAGGTTTCGTTGGTTATACCTGGTGCTAGCGATaaa                                                                                    >  2:721368/1‑143 (MQ=255)
                                                          tAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGACAATGTTTCGTTGGTTATACCTGGTGCTACCGATaaa                                                                                    >  2:651445/1‑143 (MQ=255)
                                                                    gccgCGTAAATTTCCTTCACAGGGAGATAGCTGGTCTAAGTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCt                                                                          <  2:202609/140‑1 (MQ=255)
                                                                         gcaaaTTTCCTTCCTCGGGAGATAACGGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGt                                                                     <  2:314671/141‑1 (MQ=255)
                                                                            aaTTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCa                                                                  >  1:774016/1‑143 (MQ=255)
                                                                                     cctCGGGAGTTAACTGGTATAATACCCCAGCCGTTTCTCAAAAAAAAGCCCCATGCGACGGGGGCCAAGCAAACTATCGCAATGTTTCGTTGGTTATTCCTGGCGCTAACCATAAATATCCGCGCTGGCGTGCATATTGCCGt                                                         <  2:227470/143‑1 (MQ=255)
                                                                                       ttgggAGAGAACTGGTCTAATTCGGCAGCCGGTTTTCAAAGAAAAGCCGCCTGCGAGGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTa                                                       <  2:854583/141‑1 (MQ=255)
                                                                                          ggAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTc                                                    >  2:480890/1‑143 (MQ=255)
                                                                                                aaCTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCACACTATGGCAATGTTTCGCTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCgcc                                              >  2:73536/1‑141 (MQ=255)
                                                                                                     gTCTAATCCCGCGGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCg                                         <  2:782640/143‑1 (MQ=255)
                                                                                                      tCTAATTCCGCAGCCGTTTTTCAAAAGACAGCCCCCTCCGAAGCGGGCAAAGCACACTATGGCAATGTTTCGTTGGGTATACCTGGTGCTAGCGACAAATATCCGCCCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGc                                        <  1:564941/143‑1 (MQ=255)
                                                                                                      tCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGc                                        >  1:177265/1‑143 (MQ=255)
                                                                                                      tCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTCCCCGc                                        >  1:682778/1‑143 (MQ=255)
                                                                                                       cTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTATCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCa                                       >  2:448042/1‑143 (MQ=255)
                                                                                                       cTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGCGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATCAATATCCGCGCTGGCGTGCCTATTGCCGTTACCCCCCGGTTCCCGCa                                       >  2:272849/1‑143 (MQ=255)
                                                                                                        tAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTCCTCCCCGCCTCCCGcct                                      >  1:559236/1‑141 (MQ=255)
                                                                                                        tAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCAt                                      >  2:617518/1‑143 (MQ=255)
                                                                                                         aaTTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTACATATTGCCGTTACTCCCCGGTACCCGCATc                                     >  1:523493/1‑143 (MQ=255)
                                                                                                          aTTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCa                                    <  1:108536/143‑1 (MQ=255)
                                                                                                          aTTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCa                                    >  1:137702/1‑143 (MQ=255)
                                                                                                            tCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCATACTATGGCTATTTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGa                                  <  1:651445/143‑1 (MQ=255)
                                                                                                                gCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTCTACCTGCTGCTAGCGATAAATATCCACGCTGGCGTGCATATTGCCGATACTCCCCGGTTCCCACATCAGaaat                               >  1:280230/1‑141 (MQ=255)
                                                                                                                  aGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTGATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTAc                             <  1:492883/143‑1 (MQ=255)
                                                                                                                      gTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTCCCTgg                         >  2:64718/1‑143 (MQ=255)
                                                                                                                       tttttCAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGtta                       >  1:681524/1‑143 (MQ=255)
                                                                                                                          ttCAAAAAAAAGCCCCCTGCGAAGGGGGTAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGgtagt                     <  2:52061/143‑1 (MQ=255)
                                                                                                                               aaaaaaGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGtt                <  1:756381/143‑1 (MQ=255)
                                                                                                                                   ttgCCCCCTTCGCAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTgcgc            >  1:869326/3‑143 (MQ=255)
                                                                                                                                     gCCCCCTGCGAAGTGGTCAAAGTAATCTATGGTAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAAGATCCGCGCTGGCGTGCAGATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCt          <  1:307937/143‑1 (MQ=255)
                                                                                                                                       ccccTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTg        >  1:471607/1‑143 (MQ=255)
                                                                                                                                        cccTGCGAAGGGGGCAGAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGAGCCTTgg       >  1:218640/1‑142 (MQ=255)
                                                                                                                                        cccTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTGCGTTGCGCCTTgc       >  2:657169/1‑143 (MQ=255)
                                                                                                                                             ccaaGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATc  <  2:702924/141‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                            
CCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATC  >  NC_000913/4424326‑4424608

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: