Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,989,853 T→G E305A (GAA→GCA)  hemC ← hydroxymethylbilane synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,989,8530TG98.5% 249.3 / NA 65E305A (GAA→GCA) hemChydroxymethylbilane synthase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (28/36);  total (28/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTAATTGTTGACCCGGAGAAAACTCAATCAGCGGAAAATGCCAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAG  >  NC_000913/3989707‑3989996
                                                                                                                                                  |                                                                                                                                               
ggTAATTGTTGACCCGGAGAAAACTCAATCAGCGGAAAATGCCAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGc                                                                                                                                             <  1:323416/151‑1 (MQ=255)
 gTAATTGTTGACCCGGAGAAAACTCAATCAGCGGAAAATGCCAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCg                                                                                                                                            <  1:731712/151‑1 (MQ=255)
 gTAATTGTTGACCCGGAGAAAACTCAATCAGCGGAAAATGCCAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCg                                                                                                                                            <  2:715731/151‑1 (MQ=255)
    aTTGTTGACCCGGAGAAAACTCAATCAGCGGAAAATGCCAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAgg                                                                                                                                         >  1:887167/1‑151 (MQ=255)
      tgttgACCCGGAGAAAACTCAATCAGCGGAAAATGCCAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGa                                                                                                                                        <  2:379589/150‑1 (MQ=255)
       gttgACCCGGAGAAAACTCAATCAGCGGAAAATGCCAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGAtc                                                                                                                                      <  1:135879/151‑1 (MQ=255)
       gttgACCCGGAGAAAACTCAATCAGCGGAAAATGCCAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGAtc                                                                                                                                      <  1:523721/151‑1 (MQ=255)
        ttgACCCGGAGAAAACTCAATCAGCGGAAAATGCCAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGAtc                                                                                                                                      <  1:277201/150‑1 (MQ=255)
         tgACCCGGAGAAAACTCAATCAGCGGAAAATGCCAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGAtct                                                                                                                                     <  1:251595/150‑1 (MQ=255)
         tgACCCGGAGAAAACTCAATCAGCGGAAAATGCCAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGAtct                                                                                                                                     <  1:831086/150‑1 (MQ=255)
          gACCCGGAGAAAACTCAATCAGCGGAAAATGCCAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGAtctc                                                                                                                                    <  2:257636/150‑1 (MQ=255)
            cccGGAGAAAACTCAATCAGCGGAAAATGCCAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTcgcg                                                                                                                                 >  1:817123/1‑151 (MQ=255)
              cGGAGAAAACTCAATCAGCGGAAAATGCCAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTcgcgcg                                                                                                                               <  2:182658/151‑1 (MQ=255)
                gagaAAACTCAATCAGCGGAAAATGCCAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTcgcgcgcg                                                                                                                             >  2:924132/1‑151 (MQ=255)
                    aaaCTCAATCAGCGGAAAATGCCAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGt                                                                                                                         >  2:133884/1‑151 (MQ=255)
                     aaCTCAATCAGCGGAAAATGCCAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGtt                                                                                                                        <  2:318210/151‑1 (MQ=255)
                           aTCAGCGGAAAATGCCAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAg                                                                                                                  <  2:597232/151‑1 (MQ=255)
                             cAGCGGAAAATGCCAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGCTCTCGCGCGCGCCGTTATTCAGt                                                                                                                 >  2:651211/1‑150 (MQ=255)
                              aGCGGAAAATGCCAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGt                                                                                                                 >  1:885556/1‑149 (MQ=255)
                              aGCGGAAAATGCCAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGt                                                                                                                 >  1:1059972/1‑149 (MQ=255)
                               gCGGAAAATGCCAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGc                                                                                                              >  2:6616/1‑151 (MQ=255)
                               gCGGAAAATGCCAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGc                                                                                                              <  1:178163/151‑1 (MQ=255)
                                cGGAAAATGCCAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGct                                                                                                             >  2:409158/1‑151 (MQ=255)
                                  gAAAATGCCAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGctct                                                                                                           >  1:911300/1‑151 (MQ=255)
                                      aTGCCAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGCTACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCtcttct                                                                                                        >  2:554094/1‑150 (MQ=255)
                                      aTGCCAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTg                                                                                                       >  2:1049537/1‑151 (MQ=255)
                                         ccAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGcc                                                                                                     >  2:428915/1‑150 (MQ=255)
                                                 ccTGCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATc                                                                                            <  2:72610/151‑1 (MQ=255)
                                                    gCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATcccc                                                                                         <  1:508395/151‑1 (MQ=255)
                                                    gCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATcccc                                                                                         >  1:807820/1‑151 (MQ=255)
                                                    gCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATccc                                                                                          <  2:347007/150‑1 (MQ=255)
                                                         aGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTg                                                                                    <  2:496285/151‑1 (MQ=255)
                                                              gcgcAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCgg                                                                               <  1:723308/151‑1 (MQ=255)
                                                                 cAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCAt                                                                            >  1:121307/1‑151 (MQ=255)
                                                                   gACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCt                                                                          >  1:956728/1‑151 (MQ=255)
                                                                      ggCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTg                                                                        >  2:767233/1‑150 (MQ=255)
                                                                       gCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTgcg                                                                      >  2:585911/1‑151 (MQ=255)
                                                                          cACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCAGGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCGTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTgcggcg                                                                   <  2:707027/151‑1 (MQ=255)
                                                                              aaCTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCa                                                                 <  1:298294/149‑1 (MQ=255)
                                                                              aaCTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCAcc                                                               <  2:720208/151‑1 (MQ=255)
                                                                               aCTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACcg                                                              <  2:258003/151‑1 (MQ=255)
                                                                                   ttctCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCgcggc                                                          >  1:738245/1‑151 (MQ=255)
                                                                                      tcCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGtt                                                       >  2:932784/1‑151 (MQ=255)
                                                                                       ccAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCgcggcg                                                         <  1:665333/148‑1 (MQ=255)
                                                                                       ccAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGt                                                        <  1:1030317/149‑1 (MQ=255)
                                                                                        cAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTc                                                      >  2:22141/1‑150 (MQ=255)
                                                                                            gggAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACcgc                                                 <  2:741036/151‑1 (MQ=255)
                                                                                                gaCGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCAACGCGaat                                             >  1:805393/1‑151 (MQ=255)
                                                                                                  cgggcgggTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGaataa                                           <  2:635416/151‑1 (MQ=255)
                                                                                                   gggcgggTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGaataat                                          <  1:439257/151‑1 (MQ=255)
                                                                                                       gggTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGc                                      >  1:699270/1‑151 (MQ=255)
                                                                                                           gACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAAc                                  >  1:719314/1‑151 (MQ=255)
                                                                                                             cAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCg                                <  1:166935/151‑1 (MQ=255)
                                                                                                                ggATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTcc                             <  1:430189/151‑1 (MQ=255)
                                                                                                                ggATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTcc                             <  2:986885/151‑1 (MQ=255)
                                                                                                                 gATACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCg                            <  1:812388/151‑1 (MQ=255)
                                                                                                                  aTACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCgg                           >  1:662017/1‑151 (MQ=255)
                                                                                                                   tACTCATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGgc                          >  2:1114452/1‑151 (MQ=255)
                                                                                                                       cATGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGgcgcg                       <  2:719477/150‑1 (MQ=255)
                                                                                                                        aTGCCGGGGCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGcc                     <  2:55067/151‑1 (MQ=255)
                                                                                                                               ggCGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAgc              <  1:311115/151‑1 (MQ=255)
                                                                                                                                 cGTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAgcgc            <  1:410563/151‑1 (MQ=255)
                                                                                                                                  gTCTCCGTTATAGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAgcgc            >  1:597610/1‑150 (MQ=255)
                                                                                                                                           ataGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGcc    >  1:767411/1‑149 (MQ=255)
                                                                                                                                           ataGACTGCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAg  <  2:1123848/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
GGTAATTGTTGACCCGGAGAAAACTCAATCAGCGGAAAATGCCAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAG  >  NC_000913/3989707‑3989996

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