Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 556,959 T→C A2A (GCA→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,9590TC98.9% 571.0 / ‑2.9 179A2A (GCA→GCG) folD5,10‑methylene‑tetrahydrofolate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (91/86);  total (92/87)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.20e-01

Only 100 of 189 total aligned reads displayed.

CCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTC  >  NC_000913/556869‑557044
                                                                                          |                                                                                     
ccGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTtgccgc                                                                                     >  2:3020425/1‑93 (MQ=255)
ccGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTtgccgc                                                                                     >  1:2898165/1‑93 (MQ=255)
 cGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTgccgcc                                                                                    >  1:362353/1‑93 (MQ=255)
   cAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCAt                                                                                  <  2:1632369/93‑1 (MQ=255)
    cgTCCGACTGCAATGCGCACCTGAGCTTTTTGAGCACCTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATc                                                                                 <  1:4197758/92‑1 (MQ=255)
      tCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCag                                                                               <  1:4167876/93‑1 (MQ=255)
       ccGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCATAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCaga                                                                              <  1:70282/93‑1 (MQ=255)
       ccGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCaga                                                                              >  1:1430149/1‑93 (MQ=255)
        cGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGGCGCAATAagag                                                                             >  1:1505520/1‑93 (MQ=255)
           cTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCagagaga                                                                          >  1:4535888/1‑93 (MQ=255)
            tGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCagagagag                                                                         >  2:1058544/1‑93 (MQ=255)
            tGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCagagagag                                                                         <  1:2435739/93‑1 (MQ=255)
             gCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAgg                                                                        >  1:843346/1‑93 (MQ=255)
              cAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGa                                                                       >  1:2646028/1‑93 (MQ=255)
               aaTACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGAt                                                                      >  2:4492273/1‑93 (MQ=255)
               aaTACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGAt                                                                      <  1:338247/93‑1 (MQ=255)
               aaTACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGAt                                                                      >  2:1169411/1‑93 (MQ=255)
               aaTACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGAt                                                                      <  1:1249046/93‑1 (MQ=255)
                aTACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGAtt                                                                     <  2:2969552/93‑1 (MQ=255)
                 tACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTc                                                                    >  1:2376163/1‑93 (MQ=255)
                 tACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTc                                                                    >  1:4368718/1‑93 (MQ=255)
                 tACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGATATCAGAGAGAGCATTc                                                                    >  2:3429782/1‑93 (MQ=255)
                   cgcgCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCa                                                                  <  2:1602500/93‑1 (MQ=255)
                    gcgcCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCAt                                                                 >  1:1728500/1‑93 (MQ=255)
                      gcCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCt                                                               >  2:2888331/1‑93 (MQ=255)
                        cTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGt                                                             <  1:1837547/93‑1 (MQ=255)
                        cTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGt                                                             >  1:1330495/1‑93 (MQ=255)
                         tGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGtt                                                            >  1:930186/1‑93 (MQ=255)
                         tGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGtt                                                            >  1:1794657/1‑93 (MQ=255)
                         tGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGtt                                                            <  1:4509083/93‑1 (MQ=255)
                          gAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTa                                                           >  1:1236620/1‑93 (MQ=255)
                           aaCTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTAc                                                          <  1:884737/93‑1 (MQ=255)
                            aCTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCTATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACg                                                         >  2:637689/1‑93 (MQ=255)
                              tttttGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTa                                                       >  2:4043638/1‑93 (MQ=255)
                              tttttGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTa                                                       <  2:888894/93‑1 (MQ=255)
                              tttttGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTa                                                       <  1:1064359/93‑1 (MQ=255)
                              tttttGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTa                                                       >  2:3612942/1‑93 (MQ=255)
                              tttttGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTa                                                       <  2:2808453/93‑1 (MQ=255)
                              tttgtGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTa                                                       <  2:2808454/93‑1 (MQ=255)
                               ttttGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAg                                                      <  2:2232839/93‑1 (MQ=255)
                               ttttGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAg                                                      >  2:4166594/1‑93 (MQ=255)
                               ttttGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAg                                                      >  1:478008/1‑93 (MQ=255)
                                  tGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATc                                                   <  1:1014775/93‑1 (MQ=255)
                                     gCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGaa                                                <  1:4414744/93‑1 (MQ=255)
                                      cAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAg                                               <  2:534482/93‑1 (MQ=255)
                                      cAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAg                                               >  2:1116023/1‑93 (MQ=255)
                                       aaCTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAgg                                              >  2:805461/1‑93 (MQ=255)
                                       aaCTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAgg                                              <  1:3317919/93‑1 (MQ=255)
                                        aCTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAggg                                             >  2:4590031/1‑93 (MQ=255)
                                          ttCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGa                                           >  2:572419/1‑93 (MQ=255)
                                          ttCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGa                                           <  1:1912360/93‑1 (MQ=255)
                                           tCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGAt                                          <  1:350204/93‑1 (MQ=255)
                                           tCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGAt                                          <  1:1309743/93‑1 (MQ=255)
                                           tCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGAt                                          >  2:4325668/1‑93 (MQ=255)
                                            cAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATg                                         >  1:4001505/1‑93 (MQ=255)
                                            cAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATg                                         <  2:2752710/93‑1 (MQ=255)
                                            cAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATg                                         <  1:4677309/93‑1 (MQ=255)
                                             agagCGCACCTGCTGCGCAATCGTTTTCCCGTCAATAATCTTTGCCGCCATCACAGAGAGGATTCCATCTGTTACGTACATCGTAGGGGATgg                                        >  2:1938909/1‑92 (MQ=255)
                                             agagCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATgc                                        <  1:3816635/93‑1 (MQ=255)
                                             agagCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATgc                                        >  1:2550472/1‑93 (MQ=255)
                                              gagCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTNTTACGTAGATCGAAGGGGATgcg                                       >  1:2041387/1‑93 (MQ=255)
                                              gagCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATgcg                                       >  2:4264234/1‑93 (MQ=255)
                                                   cACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTAt                                  <  1:4058280/93‑1 (MQ=255)
                                                    cccTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTAtt                                 <  1:681613/92‑1 (MQ=255)
                                                    aCCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTAtt                                 >  1:583629/1‑93 (MQ=255)
                                                    aCCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTAtt                                 >  1:3015112/1‑93 (MQ=255)
                                                     ccTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTAttt                                >  2:2505318/1‑93 (MQ=255)
                                                     ccTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAAGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTAttt                                <  2:3148150/93‑1 (MQ=255)
                                                       tgctgcGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTg                              <  1:3243423/93‑1 (MQ=255)
                                                          tgcGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCa                           >  2:3406566/1‑93 (MQ=255)
                                                           gcgcAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTCTCac                          >  1:3268700/1‑92 (MQ=255)
                                                            cgcAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGa                         <  2:3587139/93‑1 (MQ=255)
                                                               aaTCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGc                      >  2:3868386/1‑93 (MQ=255)
                                                               aaTCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGc                      <  1:2120757/93‑1 (MQ=255)
                                                               aaTCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGc                      >  1:341292/1‑93 (MQ=255)
                                                               aaTCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGc                      >  1:2769849/1‑93 (MQ=255)
                                                                aTCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGAGCCTATTTTGTCAGAAGCg                     >  2:2105714/1‑93 (MQ=255)
                                                                 tCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCgg                    >  1:4560816/1‑93 (MQ=255)
                                                                 tCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCgg                    >  1:1977513/1‑93 (MQ=255)
                                                                  cGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCggg                   <  1:188304/93‑1 (MQ=255)
                                                                   gTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCgggg                  >  2:1495582/1‑93 (MQ=255)
                                                                    ttttACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGgc                 <  2:4002047/93‑1 (MQ=255)
                                                                      ttACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGgcgc               <  1:3794318/93‑1 (MQ=255)
                                                                      ttACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGgcgc               <  2:1373543/93‑1 (MQ=255)
                                                                        aCCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATAGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGgcgcgt             >  2:3324636/1‑92 (MQ=255)
                                                                        aCCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGgcgcgc             >  2:3493391/1‑93 (MQ=255)
                                                                        aCCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGgcgcgc             <  1:3720295/93‑1 (MQ=255)
                                                                          cGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTg           <  2:1833711/93‑1 (MQ=255)
                                                                              aataatCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAg       <  1:2841374/93‑1 (MQ=255)
                                                                              aataatCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAg       >  1:2073081/1‑93 (MQ=255)
                                                                              aataatCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAg       <  2:3220410/93‑1 (MQ=255)
                                                                                taatCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGgt     <  2:132564/93‑1 (MQ=255)
                                                                                taatCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGgt     >  1:524721/1‑93 (MQ=255)
                                                                                taatCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGgt     >  2:735640/1‑93 (MQ=255)
                                                                                 aatCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGgtt    >  1:2353792/1‑93 (MQ=255)
                                                                                 aatCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGgtt    >  2:2748533/1‑93 (MQ=255)
                                                                                  atCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGgttt   >  1:3577241/1‑93 (MQ=255)
                                                                                   tCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGgtttc  >  2:3664613/1‑93 (MQ=255)
                                                                                   tCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGgtttc  <  1:1454337/93‑1 (MQ=255)
                                                                                   tCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGgtttc  <  1:4432689/93‑1 (MQ=255)
                                                                                          |                                                                                     
CCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTC  >  NC_000913/556869‑557044

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