Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,418,862 C→T G255D (GGT→GAT)  ulaG ← L‑ascorbate‑6‑phosphate lactonase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,418,8620CT98.6% 228.8 / NA 73G255D (GGT→GAT) ulaGL‑ascorbate‑6‑phosphate lactonase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (38/34);  total (39/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

GATAGTGGTACTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATACCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAGGGTGATCAGTGCAGTACGGTCGAAAGCATCAAGCGCATGAATTTCAATG  >  NC_000913/4418623‑4419104
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                  
gATAGTGGTACTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGttt                                                                                                                                                                                                                                         <  2:218036/251‑1 (MQ=255)
     tGGTACTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATAGCCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGcc                                                                                                                                                                                                                                    >  1:405153/1‑251 (MQ=255)
      ggTACTCTAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATACCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGcc                                                                                                                                                                                                                                    <  2:489776/250‑1 (MQ=255)
        tACTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGt                                                                                                                                                                                                                                  >  1:192115/1‑250 (MQ=255)
             gAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATc                                                                                                                                                                                                                            >  1:276020/1‑251 (MQ=255)
                gTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCt                                                                                                                                                                                                                          >  2:381167/1‑250 (MQ=255)
                   gTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACg                                                                                                                                                                                                                      >  1:431255/1‑251 (MQ=255)
                     cTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCCCAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGcc                                                                                                                                                                                                                    >  2:527979/1‑251 (MQ=255)
                        tATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTTACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACg                                                                                                                                                                                                                 >  2:158354/1‑251 (MQ=255)
                           ccAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTc                                                                                                                                                                                                               <  1:381167/250‑1 (MQ=255)
                            cAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGa                                                                                                                                                                                                             <  2:165234/251‑1 (MQ=255)
                                 gCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGGTTTCGCCGTACGATCCTAACGCCACGTCGATCTg                                                                                                                                                                                                         >  1:37789/1‑250 (MQ=255)
                                 gCCAGGTAAAGTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGAACACTACGTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGa                                                                                                                                                                                                        <  2:320305/251‑1 (MQ=255)
                                  ccAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGAt                                                                                                                                                                                                       <  1:457857/251‑1 (MQ=255)
                                  ccAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGAt                                                                                                                                                                                                       >  1:99461/1‑251 (MQ=255)
                                     ggTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGCCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGtt                                                                                                                                                                                                    >  2:168887/1‑251 (MQ=255)
                                        aaaTTTGCCACCCACATGCCAGATAAACGGCTTGACGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCCCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGt                                                                                                                                                                                                 >  2:92119/1‑251 (MQ=255)
                                               ccacccacCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTg                                                                                                                                                                                          >  1:399040/1‑251 (MQ=255)
                                                         gCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGACGATCTGATGTTCGTTACCGTGCTTCGCataa                                                                                                                                                                                <  2:101020/251‑1 (MQ=255)
                                                          ccAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTAGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCctaat                                                                                                                                                                               >  1:262476/1‑251 (MQ=255)
                                                                      cTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGagtagt                                                                                                                                                                   >  2:555160/1‑251 (MQ=255)
                                                                         gAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGGCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTggg                                                                                                                                                                <  1:285269/251‑1 (MQ=255)
                                                                              cATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTcg                                                                                                                                                           <  1:122416/251‑1 (MQ=255)
                                                                                tACTTCAGGCGATCTTTTTAAATCTCCCACAGCACGCGAATATCTTGCGGCTCGGCCTGGAAGTTTGCCCAGATATCGTGGTGGACCGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCCCTGTTCCTTTTGTAGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGACACGTCGATCTGCTGTTCGTTTCCGTGCTTCGCCTAATAGTTAGAGTCGTGGGAGTcgcc                                                                                                                                                         >  2:155830/1‑251 (MQ=255)
                                                                                    tCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCt                                                                                                                                                      <  1:85943/250‑1 (MQ=255)
                                                                                         cgctCTTTTTTCATCTCCCACAGCTCTAGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGCTATCGTGGTGGAACGGGATCACTCCTTTCGCATTCCGCGCTTCACCCATACGCAGAATATCGGCGCTGGTCATTTTGTCGGTGATCTCGCGCGGGTTTTCGCCGTATGATCCTAACGCCACGTCGCTCTGATGTTCGTTCCCGTGCTTCGCATATTAGTTAGTGTAGTTGTGGTCGCCGCTGTGAt                                                                                                                                                 >  2:194839/4‑250 (MQ=255)
                                                                                          gATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGGTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATAc                                                                                                                                               >  1:38846/1‑251 (MQ=255)
                                                                                            tCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATAc                                                                                                                                               >  1:262446/1‑249 (MQ=255)
                                                                                             cTTTTTTCAACTAACACAGCCCGCGAATCAATTGCGGACCGGACTGGAAGTTTGACCAGAAATCGTGGTGGAACCGGCTCAATAATTTCGCATTCAGCGCTTAACCCATACGCAGCCTATCGGCGCTGGTCATTTTGTCGGTGATCTCACGCGGGTTTTCGCCGCACGATCCTCACGCCACGTCGATCTGATGTTCGTTACCGTGCTTAGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAgg                                                                                                                                            >  2:273046/1‑251 (MQ=255)
                                                                                               tttttCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATa                                                                                                                                                <  1:528402/245‑1 (MQ=255)
                                                                                                         cccACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAgccgcc                                                                                                                                    <  1:558167/247‑1 (MQ=255)
                                                                                                           cacaGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGt                                                                                                                              >  1:327249/1‑251 (MQ=255)
                                                                                                            acaGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGtt                                                                                                                             <  2:192115/251‑1 (MQ=255)
                                                                                                             caGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGttt                                                                                                                            >  1:237118/1‑251 (MQ=255)
                                                                                                                cACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGCTCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTg                                                                                                                          >  2:328387/1‑250 (MQ=255)
                                                                                                                     ggATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGCTACAGGGAGCCGCCAGGCGTTTTGAACAgg                                                                                                                    >  2:478505/1‑251 (MQ=255)
                                                                                                                        tctcTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGt                                                                                                                   >  2:529018/1‑249 (MQ=255)
                                                                                                                              gCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCAcc                                                                                                           <  1:50625/251‑1 (MQ=255)
                                                                                                                              gCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCAcc                                                                                                           <  1:210320/251‑1 (MQ=255)
                                                                                                                                     ggCCTGGAAGTTTGACCAGAGATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGAt                                                                                                    <  2:15960/251‑1 (MQ=255)
                                                                                                                                       ccTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCg                                                                                                  >  1:366640/1‑251 (MQ=255)
                                                                                                                                          ggAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTc                                                                                                <  2:237118/250‑1 (MQ=255)
                                                                                                                                             aGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGCGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATg                                                                                            >  1:559934/1‑251 (MQ=255)
                                                                                                                                             aGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATg                                                                                            <  1:387757/251‑1 (MQ=255)
                                                                                                                                                  gACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGCTCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTAGCCGCTGCGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATc                                                                                       >  1:281934/1‑251 (MQ=255)
                                                                                                                                                   aCCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCt                                                                                      <  1:545202/251‑1 (MQ=255)
                                                                                                                                                     cAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTgg                                                                                    >  1:438725/1‑251 (MQ=255)
                                                                                                                                                      aGATATCGTGGTGGAAAGGGATCACTACTTTCGCATTCAGAGCTTAACCCATACGCAGCATATCGGCGCTGCTCATTTTGTCGGTGATACCGCGCGGGTTTTCGCCGTACGCTCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGACCAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGc                                                                                   >  2:570101/1‑251 (MQ=255)
                                                                                                                                                            cGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGGAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACg                                                                             <  2:366640/251‑1 (MQ=255)
                                                                                                                                                             gtggtgGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCAGGCAGTACGc                                                                            >  2:395962/1‑251 (MQ=255)
                                                                                                                                                                           tCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTAGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTg                                                               >  2:291321/1‑250 (MQ=255)
                                                                                                                                                                                 ctTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCggc                                                         >  2:532557/1‑250 (MQ=255)
                                                                                                                                                                                     cGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCggcaggc                                                     >  2:284093/1‑250 (MQ=255)
                                                                                                                                                                                     cGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCggcaggc                                                     <  1:291321/250‑1 (MQ=255)
                                                                                                                                                                                        aTTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTTATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTCCGCCAGCCGCTTTTTGATCGGCAGGCAggg                                                 <  1:231195/251‑1 (MQ=255)
                                                                                                                                                                                           cAGCGCTTCACCCATCCGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAGGGTGa                                              <  2:37789/251‑1 (MQ=255)
                                                                                                                                                                                                 ttCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAGGGTGATCAGTg                                        >  2:314957/1‑251 (MQ=255)
                                                                                                                                                                                                      ccATACGCAGCATAGCGGCGCTGGTCATTTTGGCGGTGATATCGCGCGGGTTTTCGCCGTACGAGCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGAGCGGCAGGCAGGGTGATCAGTGCAGTa                                   <  2:276020/251‑1 (MQ=255)
                                                                                                                                                                                                       cATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAGGGTGATCAGTGCAGTAc                                  <  2:99461/251‑1 (MQ=255)
                                                                                                                                                                                                        agtcGCAGCTTATCGACGCTGGTCATTTTGTCGTTGATATCGCGCGGGTTTTCGCCGTACTATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGATTCGCATAATAGTTAGGGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAGGGTGATCAGTGCAGTACg                                 <  2:313290/248‑1 (MQ=255)
                                                                                                                                                                                                                caTATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAGGGTGATCAGTGCAGTACGGTCGaa                           >  1:285527/1‑249 (MQ=255)
                                                                                                                                                                                                                 atatCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAGGGTGATCAGTGCAGTACGGTCGAAAGc                        <  1:555160/251‑1 (MQ=255)
                                                                                                                                                                                                                   atCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCttt                                                                  <  2:366660/207‑1 (MQ=255)
                                                                                                                                                                                                                   atCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCttt                                                                  >  1:366660/1‑207 (MQ=255)
                                                                                                                                                                                                                   atCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAGGGTGATCAGTGCAGTACGGTCGAAAGCa                       <  1:521957/250‑1 (MQ=255)
                                                                                                                                                                                                                   atCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCAGCTGGCAGTACGCCAGCCGCTTTGTGATCGGCAGGCAGGGTGATCAGTGCAGTACGGTCGAAAGCa                       <  1:475900/250‑1 (MQ=255)
                                                                                                                                                                                                                   atCGGCGATGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGACCCTAACGCCACGTCGATCTGATGTTCGTTCCAGCGCTTCGCATCCTAGCTAGAGTAGTGGGTGTCCCCGCTGTCATCCAGGGAGCCGCCCGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCAATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCCAGGTGATTAGTGCAGTACGGTCGAAAGCAt                      >  2:574962/1‑251 (MQ=255)
                                                                                                                                                                                                                    ccGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTCCCGTGCTTCGCATAATAGTTAGAGTTGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGTCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAGGGTGATCAGTGCAGTACGGTCGAAAGCATc                     <  2:121550/250‑1 (MQ=255)
                                                                                                                                                                                                                       gcgcTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGCTCGGCAGGCAGGGGGCTCAGTGCAGTACGGTCGCAAGCATCa                    >  2:543342/1‑249 (MQ=255)
                                                                                                                                                                                                                         gcTGGTCATTTTTGCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAGGGTGATCAGTGCAGTACGGTCGAAAGCATCAAgcg                <  2:431255/251‑1 (MQ=255)
                                                                                                                                                                                                                                    tGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAGGGTGATCAGTGCAGTACGGTCGAAAGCATCAAGCGCATGAATTTc      <  2:399040/250‑1 (MQ=255)
                                                                                                                                                                                                                                    tGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAGGGTGATCAGTGCAGTACGGTCGAAAGCATCAAGCGCATGAATTTc      <  2:262446/250‑1 (MQ=255)
                                                                                                                                                                                                                                       cGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAGGGTGATCAGTGCAGTACGGTCGAAAGCATCAAGCGCATGAATTTCAATg  >  1:447618/1‑251 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                  
GATAGTGGTACTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATACCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAGGGTGATCAGTGCAGTACGGTCGAAAGCATCAAGCGCATGAATTTCAATG  >  NC_000913/4418623‑4419104

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: