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% 1198.0 / ‑1.9 373A2A (GCA→GCG) folD5,10‑methylene‑tetrahydrofolate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (1/3);  new base C (208/161);  total (209/164)
Fisher's exact test for biased strand distribution p-value = 3.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 392 total aligned reads displayed.

GCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTT  >  NC_000913/556867‑557048
                                                                                            |                                                                                         
  ccGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTtgccgc                                                                                         >  1:829493/1‑93 (MQ=255)
  ccGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTtgccgc                                                                                         <  2:2766001/93‑1 (MQ=255)
   cGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTgccgcc                                                                                        >  1:2417449/1‑93 (MQ=255)
   cGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTgccgcc                                                                                        <  1:3637049/93‑1 (MQ=255)
    gCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCa                                                                                       <  1:6726827/93‑1 (MQ=255)
    gCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCa                                                                                       <  2:9754995/93‑1 (MQ=255)
    gCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCa                                                                                       >  1:972910/1‑93 (MQ=255)
    ccAGTCCGTCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCa                                                                                       <  1:3992808/92‑1 (MQ=255)
    ccAGTCCGGCTCCAATACGCGCCTGAACTTTTTGAGCACCTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCa                                                                                       <  2:594574/92‑1 (MQ=255)
      aGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATc                                                                                     >  2:1087393/1‑93 (MQ=255)
        tCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCag                                                                                   >  1:1619620/1‑93 (MQ=255)
          cGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCagag                                                                                 >  2:3134703/1‑93 (MQ=255)
           ggCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCagaga                                                                                >  2:589290/1‑93 (MQ=255)
           ggCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCagaga                                                                                >  2:7951403/1‑93 (MQ=255)
             cTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCagagaga                                                                              <  1:10037439/93‑1 (MQ=255)
               gCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAgg                                                                            >  2:9059231/1‑93 (MQ=255)
                  aTACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGAtt                                                                         <  1:338432/93‑1 (MQ=255)
                    aCGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTcc                                                                       <  1:8211026/93‑1 (MQ=255)
                    aCGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTcc                                                                       <  1:2282741/93‑1 (MQ=255)
                    aCGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTca                                                                       >  2:4619496/1‑92 (MQ=255)
                     cgcgCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCa                                                                      <  2:7305825/93‑1 (MQ=255)
                     cgcgCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCa                                                                      <  2:2862951/93‑1 (MQ=255)
                      gcgcCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCAt                                                                     >  1:2804409/1‑93 (MQ=255)
                       cgcCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCAtt                                                                    >  2:1752904/1‑92 (MQ=255)
                       cgcCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATc                                                                    <  2:5175096/93‑1 (MQ=255)
                         ccTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTg                                                                  <  2:804366/93‑1 (MQ=255)
                          cTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGt                                                                 >  1:7894858/1‑93 (MQ=255)
                            gAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTa                                                               <  1:4990146/93‑1 (MQ=255)
                            gAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTa                                                               <  1:5931223/93‑1 (MQ=255)
                               cTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGt                                                            <  1:3526849/93‑1 (MQ=255)
                               cTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGt                                                            >  1:8750396/1‑93 (MQ=255)
                                tttttGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTa                                                           >  1:1336099/1‑93 (MQ=255)
                                tttttGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTa                                                           <  1:3863192/93‑1 (MQ=255)
                                 ttttGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAg                                                          >  1:6119591/1‑93 (MQ=255)
                                 ttttGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAg                                                          >  2:8206494/1‑93 (MQ=255)
                                 ttttGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAg                                                          >  1:7742119/1‑93 (MQ=255)
                                   ttGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGAt                                                        <  2:4475927/93‑1 (MQ=255)
                                   ttGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGAt                                                        <  1:6488784/93‑1 (MQ=255)
                                     gAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCg                                                      >  2:4147018/1‑93 (MQ=255)
                                      aGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGa                                                     <  2:3295239/93‑1 (MQ=255)
                                      aGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGa                                                     >  2:3579343/1‑93 (MQ=255)
                                      aGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGa                                                     >  1:5495491/1‑93 (MQ=255)
                                       gCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGaa                                                    <  2:6286569/93‑1 (MQ=255)
                                       gCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGaa                                                    <  1:8548632/93‑1 (MQ=255)
                                       gCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGaa                                                    <  1:9530167/93‑1 (MQ=255)
                                         aaCTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAgg                                                  <  2:1616648/93‑1 (MQ=255)
                                         aaCTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAgg                                                  >  1:9942037/1‑93 (MQ=255)
                                         aaCTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAgg                                                  >  2:2313557/1‑93 (MQ=255)
                                         aaCTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGCAGATCGAAgg                                                  >  1:631212/1‑93 (MQ=255)
                                          aCTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAggg                                                 >  1:8113344/1‑93 (MQ=255)
                                          aCTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAggg                                                 >  1:7198713/1‑93 (MQ=255)
                                           cTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAgggg                                                >  2:1433857/1‑93 (MQ=255)
                                            ttCAGAGCGCCCCTGCTGGGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGa                                               <  2:6958373/93‑1 (MQ=255)
                                             tCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGAt                                              <  1:1086407/93‑1 (MQ=255)
                                                gagCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATgcg                                           >  2:6559132/1‑93 (MQ=255)
                                                 agCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATgcgc                                          >  2:2032334/1‑93 (MQ=255)
                                                 agCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATgcgc                                          >  2:6748329/1‑93 (MQ=255)
                                                   cgcACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCt                                        <  1:3045295/93‑1 (MQ=255)
                                                   cgcACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCt                                        <  1:9178634/93‑1 (MQ=255)
                                                     cACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTAt                                      >  2:4761779/1‑93 (MQ=255)
                                                      cccTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTAtt                                     <  2:3520740/92‑1 (MQ=255)
                                                      aCCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTAtt                                     <  1:6594104/93‑1 (MQ=255)
                                                        ctgctgCGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTAtttt                                   >  1:4636675/1‑93 (MQ=255)
                                                         tgctgcGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTAttttt                                  >  1:3159607/1‑92 (MQ=255)
                                                         tgctgcGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTg                                  >  2:9888308/1‑93 (MQ=255)
                                                         tgctgcGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTg                                  >  2:4887496/1‑93 (MQ=255)
                                                           ctgcGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTc                                >  1:4268123/1‑93 (MQ=255)
                                                            tgcGCAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCa                               <  1:6058242/93‑1 (MQ=255)
                                                             gcgcAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGTAGGGGATGCGCCTATTTTGTCAg                              >  2:5518703/1‑93 (MQ=255)
                                                             gcgcAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAg                              >  2:9373689/1‑93 (MQ=255)
                                                              cgcAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGa                             <  2:3894863/93‑1 (MQ=255)
                                                              cgcAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGa                             >  1:5787991/1‑93 (MQ=255)
                                                               gcAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGaa                            <  2:7687384/93‑1 (MQ=255)
                                                                cAATCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAg                           <  2:4052045/93‑1 (MQ=255)
                                                                 aaTCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGc                          <  1:5215548/93‑1 (MQ=255)
                                                                 aaTCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGc                          >  2:4211876/1‑93 (MQ=255)
                                                                 aaTCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCATATTTTGTCAGaatg                          >  2:2548375/1‑91 (MQ=255)
                                                                  aTCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCg                         >  2:8799867/1‑93 (MQ=255)
                                                                  aTCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCg                         >  1:5917399/1‑93 (MQ=255)
                                                                  aTCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCg                         >  2:8589677/1‑93 (MQ=255)
                                                                  aTCGTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCg                         >  1:1956508/1‑93 (MQ=255)
                                                                    tgTTTTACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCggg                       <  2:5718098/92‑1 (MQ=255)
                                                                        ttACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGgcgc                   >  1:8039146/1‑93 (MQ=255)
                                                                         tACCGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGgcgcg                  >  2:4354428/1‑93 (MQ=255)
                                                                            cGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTg               >  1:5575638/1‑93 (MQ=255)
                                                                            cGTCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTg               >  2:6125049/1‑93 (MQ=255)
                                                                              tCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTc             >  2:5343599/1‑93 (MQ=255)
                                                                              tCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTc             >  1:3298662/1‑93 (MQ=255)
                                                                              tCAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTc             >  1:4556635/1‑93 (MQ=255)
                                                                               taataatCTTGGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCa            <  1:7227381/92‑1 (MQ=255)
                                                                               cAATAATCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCa            >  1:4846512/1‑93 (MQ=255)
                                                                                 ataatCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAgg          >  1:6938477/1‑93 (MQ=255)
                                                                                  taatCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGACGCGCCTATTTTGTCATAAGCGGGGGGCGCTGTCAGgt         >  2:8395784/1‑93 (MQ=255)
                                                                                   aatCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGgtt        >  2:3795103/1‑93 (MQ=255)
                                                                                    atCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGgttt       <  2:9561067/93‑1 (MQ=255)
                                                                                    atCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGgttt       <  2:124041/93‑1 (MQ=255)
                                                                                    atCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGgttt       >  1:9358498/1‑93 (MQ=255)
                                                                                    atCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGgttt       >  2:7473804/1‑93 (MQ=255)
                                                                                    atCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGgttt       <  1:611141/93‑1 (MQ=255)
                                                                                     tCTTTGCCGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGgtttc      >  2:8681181/1‑93 (MQ=255)
                                                                                            |                                                                                         
GCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTT  >  NC_000913/556867‑557048

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