Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,249,670 A→G V131A (GTC→GCC)  b1200 ← putative dihydroxyacetone kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,249,6700AG99.1% 339.4 / NA 111V131A (GTC→GCC) b1200putative dihydroxyacetone kinase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (51/59);  minor base T (1/0);  total (52/59)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 112 total aligned reads displayed.

CTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGACCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATC  >  NC_000913/1249581‑1249759
                                                                                         |                                                                                         
cTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCca                                                                                        >  1:4136981/1‑93 (MQ=255)
cTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCca                                                                                        <  1:2350142/93‑1 (MQ=255)
 tACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCcac                                                                                       >  2:1121075/1‑93 (MQ=255)
 tACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCcac                                                                                       >  1:4567501/1‑93 (MQ=255)
 tACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCcac                                                                                       >  1:4524879/1‑93 (MQ=255)
  aCGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCcaca                                                                                      >  1:3236114/1‑93 (MQ=255)
  aCGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCcaca                                                                                      >  1:2584449/1‑93 (MQ=255)
   cGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAg                                                                                     <  2:3161863/93‑1 (MQ=255)
   cGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAg                                                                                     <  1:147655/93‑1 (MQ=255)
   cGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAg                                                                                     <  2:4204514/93‑1 (MQ=255)
   cGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTGCCAGCAACGTGGTCATCAATGGCCACAg                                                                                     >  1:1865199/1‑93 (MQ=255)
    gAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGt                                                                                    <  1:4041421/93‑1 (MQ=255)
     aGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTg                                                                                   <  1:657555/93‑1 (MQ=255)
      gTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTgg                                                                                  >  2:767854/1‑93 (MQ=255)
      gTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTgg                                                                                  <  1:654723/93‑1 (MQ=255)
      gTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTgg                                                                                  >  1:3809143/1‑93 (MQ=255)
       tttttCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGt                                                                                 <  1:1517207/93‑1 (MQ=255)
       tttttCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGt                                                                                 >  2:3985719/1‑93 (MQ=255)
        ttttCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTc                                                                                <  1:1501366/93‑1 (MQ=255)
        ttttCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTc                                                                                <  1:3067318/93‑1 (MQ=255)
          ttCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTCAc                                                                              <  1:568103/93‑1 (MQ=255)
               ttAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTa                                                                         <  1:793789/93‑1 (MQ=255)
                tAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTAc                                                                        <  2:2600228/93‑1 (MQ=255)
                tAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTAc                                                                        <  1:996579/93‑1 (MQ=255)
                 aaTACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTAcg                                                                       >  2:3035359/1‑93 (MQ=255)
                  aTACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTAcgc                                                                      <  1:2203885/93‑1 (MQ=255)
                   tACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTAcgcc                                                                     >  2:2902233/1‑93 (MQ=255)
                      ggTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCt                                                                  >  1:801072/1‑93 (MQ=255)
                      ggTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCt                                                                  <  2:3994134/93‑1 (MQ=255)
                       gtgtTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTa                                                                 <  1:1641114/93‑1 (MQ=255)
                        tgtTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCCCAGTGGTCACTTTTACGCGGCTAt                                                                >  2:3677554/1‑93 (MQ=255)
                         gtTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATc                                                               <  1:2775680/93‑1 (MQ=255)
                         gtTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATc                                                               <  1:19262/93‑1 (MQ=255)
                         gtTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATc                                                               >  2:3197331/1‑93 (MQ=255)
                         gtTGGCAACGCCGCGTCGCCAGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATc                                                               <  1:4153757/93‑1 (MQ=255)
                             gCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGc                                                           >  2:654171/1‑93 (MQ=255)
                                  gccgcGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTaa                                                      >  1:2870262/1‑93 (MQ=255)
                                   ccgcGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAAc                                                     <  1:1833632/93‑1 (MQ=255)
                                   ccgcGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAAc                                                     <  1:2683346/93‑1 (MQ=255)
                                     gcgTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTc                                                   <  1:219317/93‑1 (MQ=255)
                                       gTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCgg                                                 >  2:589497/1‑93 (MQ=255)
                                        tCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGt                                                <  2:3103423/93‑1 (MQ=255)
                                        tCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGt                                                <  2:2214353/93‑1 (MQ=255)
                                         cGCCCGGCATTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTc                                               <  1:321857/93‑1 (MQ=255)
                                           cccGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGc                                             <  2:1405955/93‑1 (MQ=255)
                                              ggCAGTATAAAGGCTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGt                                          <  2:2258139/93‑1 (MQ=255)
                                              ggCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGt                                          <  2:4242965/93‑1 (MQ=255)
                                                 aGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTc                                       <  1:2962760/93‑1 (MQ=255)
                                                  gTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCa                                      <  2:2411566/93‑1 (MQ=255)
                                                   tataAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCaa                                     >  2:2027297/1‑93 (MQ=255)
                                                   tataAAGACTGGCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTCTCaa                                     >  2:233823/1‑93 (MQ=255)
                                                     taAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAg                                   >  1:47833/1‑93 (MQ=255)
                                                     taAAGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAg                                   <  1:451506/93‑1 (MQ=255)
                                                      aaaGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGt                                  >  1:3944443/1‑93 (MQ=255)
                                                        aGACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTa                                >  1:4184587/1‑93 (MQ=255)
                                                         gACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTaa                               >  1:4076011/1‑93 (MQ=255)
                                                         gACTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTaa                               >  1:2011986/1‑93 (MQ=255)
                                                           cTGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGa                             >  2:430359/1‑93 (MQ=255)
                                                            tGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGaa                            <  2:3217458/93‑1 (MQ=255)
                                                            tGTCTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGaa                            >  2:1469097/1‑93 (MQ=255)
                                                               cTTTTACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAatat                         <  2:1453064/93‑1 (MQ=255)
                                                                ttttACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATc                        <  1:4631825/93‑1 (MQ=255)
                                                                  ttACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATCGc                      <  1:3004913/93‑1 (MQ=255)
                                                                   tACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATCGcc                     <  1:3694314/93‑1 (MQ=255)
                                                                   tACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATCGcc                     <  2:3908115/93‑1 (MQ=255)
                                                                   tACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATCGcc                     <  1:2093352/93‑1 (MQ=255)
                                                                   tACCGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATCGcc                     <  1:1575510/93‑1 (MQ=255)
                                                                     ccGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATCGCCgg                   >  2:1492826/1‑93 (MQ=255)
                                                                     ccGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATCGCCgg                   <  2:1613736/93‑1 (MQ=255)
                                                                      cGCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATCGCCGgt                  >  2:55654/1‑93 (MQ=255)
                                                                       gCAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATCGCCGgtg                 >  2:4616675/1‑93 (MQ=255)
                                                                        cAACGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATCGCCGgtgt                >  2:3439544/1‑93 (MQ=255)
                                                                         aaCGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTa               >  1:2523981/1‑93 (MQ=255)
                                                                          aCGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTaa              <  1:1236262/93‑1 (MQ=255)
                                                                          aCGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTaa              >  1:893938/1‑93 (MQ=255)
                                                                          aCGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTaa              <  1:1878850/93‑1 (MQ=255)
                                                                          aCGTCGTCATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTaa              >  1:325146/1‑93 (MQ=255)
                                                                             tcgtcATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAAttt           >  2:851869/1‑93 (MQ=255)
                                                                             tcgtcATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAAttt           >  2:1266743/1‑93 (MQ=255)
                                                                              cgtcATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAAtttt          <  2:3360717/93‑1 (MQ=255)
                                                                               gtcATCAATGTCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAAttttt         >  1:1215723/1‑93 (MQ=255)
                                                                               gtcATCAATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAAttttt         >  2:2416413/1‑93 (MQ=255)
                                                                                tcatcaATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTg        >  2:1994696/1‑93 (MQ=255)
                                                                                tcatcaATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTg        >  1:1318267/1‑93 (MQ=255)
                                                                                tcatcaATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTg        >  1:1688408/1‑93 (MQ=255)
                                                                                tcatcaATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTg        >  2:477483/1‑93 (MQ=255)
                                                                                tcatcaATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTg        >  1:1875264/1‑93 (MQ=255)
                                                                                tcatcaATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTg        >  2:3368853/1‑93 (MQ=255)
                                                                                 catcaATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGa       >  1:3377663/1‑93 (MQ=255)
                                                                                 catcaATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGa       >  1:708977/1‑93 (MQ=255)
                                                                                 catcaATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGa       <  2:333885/93‑1 (MQ=255)
                                                                                  atcaATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGAt      <  2:1606999/93‑1 (MQ=255)
                                                                                   tcaATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATa     >  2:4197900/1‑93 (MQ=255)
                                                                                   tcaATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATa     >  2:4590900/1‑93 (MQ=255)
                                                                                    caATGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATaa    <  1:2285440/93‑1 (MQ=255)
                                                                                     aaTGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAAt   >  1:334349/1‑93 (MQ=255)
                                                                                      aTGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATc  <  1:4082/93‑1 (MQ=255)
                                                                                      aTGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATc  >  1:4523423/1‑93 (MQ=255)
                                                                                      aTGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATc  <  1:1820538/93‑1 (MQ=255)
                                                                                      aTGGCCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATc  <  2:3151357/93‑1 (MQ=255)
                                                                                         |                                                                                         
CTACGAGTTTTTCAATTAATACGGTGTTGGCAACGCCGCGTCGCCCGGCAGTATAAAGACTGTCTTTTACCGCAACGTCGTCATCAATGACCACAGTGGTCACTTTTACGCCGCTATCGTGCAGTAACTCGGTCGCTGTTTCAAAGTTAAGAATATCGCCGGTGTAATTTTTGATAATC  >  NC_000913/1249581‑1249759

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: