Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 556,636 A→G V110A (GTG→GCG)  folD ← 5,10‑methylene‑tetrahydrofolate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913556,6360AG97.2% 790.2 / NA 256V110A (GTG→GCG) folD5,10‑methylene‑tetrahydrofolate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  major base G (174/73);  minor base T (0/7);  total (174/82)
Fisher's exact test for biased strand distribution p-value = 2.55e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 339 total aligned reads displayed.

GGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCC  >  NC_000913/556546‑556725
                                                                                          |                                                                                         
ggACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCt                                                                                         <  1:2721123/93‑1 (MQ=255)
 gACGCAGACGCGGCGCGCGCTGGGACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCtt                                                                                        <  2:4735377/93‑1 (MQ=255)
 gACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCtt                                                                                        >  1:7029177/1‑93 (MQ=255)
  aCGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCttt                                                                                       <  2:1104368/93‑1 (MQ=255)
   cGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTg                                                                                      <  2:423544/93‑1 (MQ=255)
        gcgcgGCGCGCCCTGGCTCAGACGACCTACCTTGTAAGTATGCAACCCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCTCTTTGACGtt                                                                                 <  1:1909033/92‑1 (MQ=255)
        aCGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGtt                                                                                 <  1:4364692/93‑1 (MQ=255)
           cggcgCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATc                                                                              >  2:2203262/1‑93 (MQ=255)
           cggcgCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAAACGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATc                                                                              >  2:7694922/1‑93 (MQ=255)
               gcgcgcTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATa                                                                          <  2:4488101/93‑1 (MQ=255)
               gcgcgcGGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATa                                                                          <  1:1759967/93‑1 (MQ=255)
                cgcgcTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATAc                                                                         >  2:6857133/1‑93 (MQ=255)
                  cgcTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATAccc                                                                       >  2:5412806/1‑93 (MQ=255)
                     tGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACccgcc                                                                    <  1:8092047/93‑1 (MQ=255)
                      ggCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACccgccg                                                                   <  1:3389163/93‑1 (MQ=255)
                      ggCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACccgccg                                                                   >  1:4444716/1‑93 (MQ=255)
                      ggCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACccgccg                                                                   >  2:2520142/1‑93 (MQ=255)
                       gCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCgg                                                                  >  1:2206538/1‑93 (MQ=255)
                       gCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCgg                                                                  >  1:833086/1‑93 (MQ=255)
                         acaGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTa                                                                <  1:4301789/93‑1 (MQ=255)
                          caGCCGACAGACGTTGTAAGGATGGAAACCGCCCACGTCTGTGCCCGAATGAATCCGTCCCAGCCCTTTGACGTTATCAATCCCCGCCGGTaa                                                               <  2:1010132/93‑1 (MQ=255)
                          caGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTaa                                                               >  1:928757/1‑93 (MQ=255)
                           aGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTCCCAGCCCTTTGACGTTATCAATCCCCGCCGGTAAc                                                              <  1:1008530/93‑1 (MQ=255)
                            gacgacCGACGTTGTAAGGATGGAACCCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACg                                                             <  2:2828944/93‑1 (MQ=255)
                            gacgacCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACg                                                             >  2:7453682/1‑93 (MQ=255)
                             acgacCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACgg                                                            >  2:1688798/1‑93 (MQ=255)
                               gaccgacGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGcg                                                          >  1:6597035/1‑92 (MQ=255)
                               gaccgacGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCa                                                          <  2:5819984/93‑1 (MQ=255)
                                 ccgacGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGt                                                        >  2:3373981/1‑93 (MQ=255)
                                 ccgacGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGt                                                        <  1:5471890/93‑1 (MQ=255)
                                   gacGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTg                                                      >  2:2716218/1‑93 (MQ=255)
                                   gacGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTg                                                      >  1:3851575/1‑93 (MQ=255)
                                   gacGTTGTAAGGATGAAACCCTTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTg                                                      <  1:5925026/93‑1 (MQ=255)
                                    acGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGa                                                     >  1:2511953/1‑93 (MQ=255)
                                       ttGTAAGGAGGGAACCCGTCCACGTCTTTGTCCGGATAAATACGTCCCAGCGTTTTGACTTTATCAATCCCCGCCGGTAACGGCAGTTGAAcc                                                  <  2:3601141/93‑1 (MQ=255)
                                        tGTAAGGATGAAACCCGTCCACGTTTTTGTCCGGATGAATCCTTTCCAGCTTTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCa                                                 <  2:8109244/93‑1 (MQ=255)
                                        tGTAAGGATGAAAACCGTCCACGTCTTTGTCGGGATAAATACTTTCCAGCGTTTTGACTTTATCAATACCCGCCGGTAACGGCAGTTGAACCa                                                 <  2:6403442/93‑1 (MQ=255)
                                         gTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAg                                                <  1:5827398/93‑1 (MQ=255)
                                          tAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGa                                               >  1:50049/1‑93 (MQ=255)
                                          tAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGa                                               >  1:3172936/1‑93 (MQ=255)
                                            atgATGGAAACCGTCCACGTCGTTGTCCGCATGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGGCGGTAACGGCAGTTGAACCCGAAt                                             >  2:4834521/3‑93 (MQ=255)
                                            aGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAAt                                             >  1:4572762/1‑93 (MQ=255)
                                             ggATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATg                                            >  1:6305043/1‑93 (MQ=255)
                                              gATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCCGAATGc                                           >  2:3232185/1‑93 (MQ=255)
                                                tGGAACCCGTCCGCGTCTTTGTCGGGATGAATCCTTTCCAGCGCTTTGACGTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCa                                         <  2:8605079/93‑1 (MQ=255)
                                                tGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCTCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCa                                         <  2:1449928/93‑1 (MQ=255)
                                                 ggAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCAt                                        <  1:7661022/93‑1 (MQ=255)
                                                  gAAACCGTCCCCGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATc                                       <  1:2312835/93‑1 (MQ=255)
                                                  gAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATc                                       >  1:2280135/1‑93 (MQ=255)
                                                   aaaCCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCg                                      >  2:3794165/1‑93 (MQ=255)
                                                    acccGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGa                                     <  1:2571381/91‑1 (MQ=255)
                                                    aaCCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGa                                     >  1:6254053/1‑93 (MQ=255)
                                                     aCCGTCCACGTCTTTGCCGGGATGAATGCGTTCCAGCGCTTTGACGTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGAt                                    <  1:2953888/93‑1 (MQ=255)
                                                        gTCCACTTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGgt                                 <  2:5004207/93‑1 (MQ=255)
                                                         tCCACGTCTTTGTCCGGATGAATACGTTCCAGCTCTTGGACGTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGgtg                                <  1:120009/93‑1 (MQ=255)
                                                         tCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGgtg                                <  1:8704980/93‑1 (MQ=255)
                                                         tCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGgtg                                >  2:7356483/1‑93 (MQ=255)
                                                         tCCACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGgtg                                >  2:759230/1‑93 (MQ=255)
                                                         cccccGTTTTTGCCGGGATCAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGgtg                                <  1:2012978/89‑1 (MQ=255)
                                                           cACGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGtgtt                              >  1:3269054/1‑93 (MQ=255)
                                                           cACGTCTTTGTCCGGATAAATACTTTCCAGCGTTTTGACTTTATAAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGtgtt                              <  2:2024869/93‑1 (MQ=255)
                                                            aCGTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGtgttg                             >  2:3424138/1‑93 (MQ=255)
                                                              gTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTc                           >  1:1580929/1‑93 (MQ=255)
                                                              gTCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTc                           >  1:6958529/1‑93 (MQ=255)
                                                               tCTTTGTCGGGATGAATACGTTCCAGCGTTTTGACGTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCg                          <  2:165551/93‑1 (MQ=255)
                                                               tCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCg                          >  2:4370361/1‑93 (MQ=255)
                                                               tCTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCg                          >  2:8345826/1‑93 (MQ=255)
                                                                tgttGTCCGGATAAATACGTTCCAGCGTTTTGACGTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCgg                         <  1:7230030/91‑1 (MQ=255)
                                                                cTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCgg                         >  1:3094111/1‑93 (MQ=255)
                                                                cTTTGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGAGTTGGCgg                         >  1:4188005/1‑93 (MQ=255)
                                                                 tttGTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCgga                        >  2:3968356/1‑92 (MQ=255)
                                                                    gcccGGATGAATACGTTCCAGCGCTTTGACTTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCAtt                     <  1:1998051/91‑1 (MQ=255)
                                                                    gTCCGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCAtt                     >  1:3711542/1‑93 (MQ=255)
                                                                     tCCGGATAAATACTTTCCAGCGCTTTGACGTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTc                    <  2:6028336/93‑1 (MQ=255)
                                                                      ccGGATGAATACGTTCCAGCTCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCa                   <  2:4593442/93‑1 (MQ=255)
                                                                      ccGGATGAATACGTTCCAGCGCTTTGACGTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCa                   <  1:1365374/93‑1 (MQ=255)
                                                                      ccGGATGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCa                   >  1:7516335/1‑93 (MQ=255)
                                                                      ccGAATGAATACGTTCCAGCTTTTTGATGTTATCAACCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCa                   <  2:1898758/93‑1 (MQ=255)
                                                                       cGGATGAATACGTTCCAGCGCTTTGACGTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAg                  <  1:6155947/93‑1 (MQ=255)
                                                                        ggATGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGc                 >  1:6805862/1‑93 (MQ=255)
                                                                        ggATAAATACGTTCCAGCGCTTTGACGTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGc                 <  1:8219805/93‑1 (MQ=255)
                                                                         gATGAATACTTTCCAGCGTTTTGACGTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCg                <  2:6578285/93‑1 (MQ=255)
                                                                         gATGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCg                <  1:3982366/93‑1 (MQ=255)
                                                                          aTGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGt               >  2:7238360/1‑93 (MQ=255)
                                                                          aTGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGt               >  2:1761272/1‑93 (MQ=255)
                                                                          aTGAATACGTGCCAGCGTTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGt               <  2:5677822/93‑1 (MQ=255)
                                                                           tGAATACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTa              >  2:1425708/1‑93 (MQ=255)
                                                                             aaTACTTCCCAGCTCTTGGACTTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATc            <  1:3531869/93‑1 (MQ=255)
                                                                             aaTACGTTCCAGCGTTTTGACGTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCAATCGATGGTGTTGTCGGCATTCAGCGTATc            <  2:8746417/93‑1 (MQ=255)
                                                                             aaTACGTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATc            >  2:3892612/1‑93 (MQ=255)
                                                                               tccGTTCCGCCGCTTTGACGTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGa          <  2:1297122/91‑1 (MQ=255)
                                                                                 cGTTCCAGCGCTTTGACGTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATa        <  1:6503437/93‑1 (MQ=255)
                                                                                  gTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATaa       >  1:5187509/1‑93 (MQ=255)
                                                                                  gTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATaa       >  2:1075090/1‑93 (MQ=255)
                                                                                  gTTCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATaa       >  2:1126410/1‑93 (MQ=255)
                                                                                  ctcccAGCGCTTTGCGTTTAACAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATaa       <  1:7906443/90‑1 (MQ=255)
                                                                                   ttCCAGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAg      >  1:463832/1‑93 (MQ=255)
                                                                                      cAGCGCTTTGACGTTATCAATACACGCCGGTAACGGCAGGTCAAACATAAAGCCCTCGATGGTGTTGTCGGCATTCAGCGCATCGATAAGCta   >  2:7934846/1‑92 (MQ=255)
                                                                                       aGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTcc  >  2:3032231/1‑93 (MQ=255)
                                                                                       aGCGCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTcc  >  1:8524853/1‑93 (MQ=255)
                                                                                          |                                                                                         
GGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCC  >  NC_000913/556546‑556725

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑