Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,823,777 A→T 5.5% E88V (GAA→GTA)  cho → endonuclease of nucleotide excision repair

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,823,7770AT 5.5% 236.9 / 3.4 73E88V (GAA→GTA) choendonuclease of nucleotide excision repair
Reads supporting (aligned to +/- strand):  ref base A (35/34);  new base T (2/2);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.64e-01

TCATGGTGAAAGTGACACCATGCCGCTCTATATCGGCAAAAGCGTTAACATCCGCAGCCGCGTCCTTTCTCATTTACGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGTTTAATAAACGGTTGCGCCGCAATCGCCAGCTCTGTGCCCTGCAATTAAATGAAAAGCGCGTCGATGTGGTGTATGCCAAAGAGGTGGATTTTTC  >  NC_000913/1823634‑1823924
                                                                                                                                               |                                                                                                                                                   
tCATGGTGAAAGTGACACCATGCCGCTCTATATCGGCAAAAGCGTTAACATCCGCAGCCGCGTCCTTTCTCATTTACGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATcg                                                                                                                                                <  3:70821‑M2/149‑3 (MQ=255)
  aTGGTGAAAGTGACACCATGCCGCTCTATATCGGCAAAAGCGTTAACATCCGCCGCCGCGTCCTTTCTCATTTACGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATcggc                                                                                                                                              <  5:236879‑M2/149‑5 (MQ=255)
  aTGGTGAAAGTGACACCATGCCGCTCTATATCGGCAAAAGCGTTAACATCCGCAGCCGCGTCCTTTCTCATTTACGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATcggc                                                                                                                                              >  6:149057‑M2/1‑145 (MQ=255)
  aTGGTGAAAGTGACACCATGCCGCTCTATATCGGCAAAAGCGTTAACATCCGCAGCCGCGTCCTTTCTCATTTACGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATcgg                                                                                                                                               >  8:52969‑M2/1‑145 (MQ=255)
          aGTGACACCATGCCGCTCTATATCGGCAAAAGCGTTAACATCCGCAGCCGCGTCCTTTCTCATTTACGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCt                                                                                                                                      >  8:243345/1‑149 (MQ=255)
           gTGACACCATGCCGCTCTATATCGGCAAAAGCGTTAACATCCGCAGCCGCGTCCTTTCTCATTTACGTACCCCGGATGAAGCCGCCATGATACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGCTGAAATCGGCGCTCAGCta                                                                                                                                     >  3:97204/1‑148 (MQ=255)
            tGACACCATGCCGCTCTATATCGGCAAAAGCGTTAACATCCGCAGCCGCGTCCTTTCTCATTTACGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTcc                                                                                                                                    >  2:333868/1‑149 (MQ=255)
                 ccATGCCGCTCTATATCGGCAAAAGCGTTAACATCCGCAGCCGCGGCCTTTCTCATATCCGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGaa                                                                                                                               <  6:161159/149‑1 (MQ=255)
                  cATGCCGCTCTATATCGGCAAAAGCGTTAACATCCGCAGCCGCGTCCTTTCTCATTTACGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTg                                                                                                                                 <  7:230839/146‑1 (MQ=255)
                  cATGCCGCTCTATATCGGCAAAAGCGTTAACATCCGCAGCCGCGTCCTTTCTCATTTACGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTg                                                                                                                                 >  8:230839/1‑146 (MQ=255)
                       cGCTCTATATCGGCAAAAGCGTTAACATCCGCAGCCGCGTCCTTTCTCATTTACGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGa                                                                                                                         <  7:77741/149‑1 (MQ=255)
                         ctctATATCGGCAAAAGCGTTAACATCCGCAGCCGCGTCCGTTCTCATTTACGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGAt                                                                                                                        >  2:47256/1‑148 (MQ=255)
                          tctATATCGGCAAAAGCGTTAACATCCGCAGCCGCGTCCTTTCTCATTTACGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGAGCAGCTGGAGATGTACAGCCGGTGAAATCGGCGCTCTGCTCCGTGACGCGCGAtg                                                                                                                       >  4:170357/1‑147 (MQ=255)
                           ctATATCGGCAAAAGCGTTAACATCCGCAGCCGCGTCCTTTCTCATTTACGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTaa                                                                                                                     >  2:262766/1‑149 (MQ=255)
                            tataTCGGCAAAAGCGTTAACATCCGCAGCCGCGTCCTTTCTCATTTACGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTc                                                                                                                                     <  2:43226/132‑1 (MQ=255)
                            tataTCGGCAAAAGCGTTAACATCCGCAGCCGCGTCCTTTCTCATTTACGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTc                                                                                                                                     >  1:43226/1‑132 (MQ=255)
                            tataTCGGCAAAAGCGTTAACATCCGCAGCCGCGTCCTTTCTCATTTACGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAAt                                                                                                                    >  1:137704/1‑149 (MQ=255)
                               atCGGCAAAAGCGTTAACATCCGCAGCCGCGTCCTTTCTCATTTACGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAAt                                                                                                                    <  7:117138/146‑1 (MQ=255)
                               atCGGCAAAAGCGTTAACATCCGCAGCCGCGTCCTTTCTCATTTACGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAAt                                                                                                                    >  8:117138/1‑146 (MQ=255)
                                  ggCAAAAGCGTTAACATCCGCAGCCGCGTCCTTTCTCATTTACGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGAt                                                                                                                        >  6:168461/1‑139 (MQ=255)
                                  ggCAAAAGCGTTAACATCCGCAGCCGCGTCCTTTCTCATTTACGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGAt                                                                                                                        <  5:168461/139‑1 (MQ=255)
                                   gCAAAAGCGTTAACATCCGCAGCCGCGTCCTTTCTCATTTACGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGaa                                                                                                             >  7:2725/1‑149 (MQ=255)
                                     aaaaGCGTTAACATCCGCAGCCGCGTCCTTTCTCATTTACGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGaac                                                                                                            >  1:57117/1‑148 (MQ=255)
                                       aaGCGTTAACATCCGCAGCCGCGTCCTTTCTCATTTACGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGaa                                                                                                                               >  7:233870/1‑127 (MQ=255)
                                       aaGCGTTAACATCCGCAGCCGCGTCCTTTCTCATTTACGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGaa                                                                                                                               <  8:233870/127‑1 (MQ=255)
                                              aaCATCCGCAGCCGCGTCCTTTCTCATTTACGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGATTTCACcggcggtacatatccagctgatccgtcgggattgccgtag                                                                                                          <  1:49981‑M2/141‑41 (MQ=255)
                                              aaCATCCGCAGCCGCGTCCTTTCTCATTTACGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGATTTCACcggcggtacatatccagctgatccgtcgggattgccgtag                                                                                                          >  2:49981‑M2/1‑101 (MQ=255)
                                               aCATCCGCAGCCGCGTCCTTTCTCATTTACGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCt                                                                                                  <  4:70821/148‑1 (MQ=255)
                                                cATCCGCAGCCGCGTCCTTTCTCATTTACGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGt                                                                                                <  7:243345/149‑1 (MQ=255)
                                                 aTCCGCAGCCGCGTCCTTTCTCATTTACGTACCCCGGATGCAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGATTTCACcggcggtacatatccagctga                                                                                                                             >  5:134752‑M2/1‑98 (MQ=255)
                                                 aTCCGCAGCCGCGTCCTTTCTCATTTACGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGATTTCACcggcggtacatatccagctga                                                                                                                             <  6:134752‑M2/119‑22 (MQ=255)
                                                    cGCAGCCGCGTCCTTTCTCATTTACGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTTCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGTTtaa                                                                                            <  5:109397/149‑1 (MQ=255)
                                                    cGCAGCCGCGTCCTTTCTCATTTACGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGTTtaa                                                                                            <  2:96002/149‑1 (MQ=255)
                                                     gCAGCCGCGTCCTTTCTCATTTTCGTACCCCGGATGCAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTCGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGTTtaat                                                                                           <  6:90528/149‑1 (MQ=255)
                                                     gCAGCCGCGTCCTTTCTCATTTACGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGTTtaat                                                                                           <  6:52852/149‑1 (MQ=255)
                                                         ccGCGTCCTTTCTCATTTACGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGTTTAATAAAc                                                                                       >  1:55680/1‑149 (MQ=255)
                                                           gcgTCCTTTCTCATTTACGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGTTTAATAAACgg                                                                                     <  1:333868/149‑1 (MQ=255)
                                                             gTCCTTTCTCATTTACGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATTAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACACCAGCCGCTGTTTAATAAACGGtt                                                                                   <  7:181193/149‑1 (MQ=255)
                                                                cTTTCTCATTTACGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGTTTAATAAACGGTTgcg                                                                                <  8:2725/149‑1 (MQ=255)
                                                                   tctcATTTACGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGTTTAATAAACGGTTGcgcc                                                                              <  1:47256/148‑1 (MQ=255)
                                                                           aCGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACATCAGGTGAAATCGGCGCTCTGCTCCTTGAAGAGCGATTAATCAAAGAACAACAGCCGCTGTTTAATAAACTGTTGCGCCGCAATCGc                                                                      >  4:110596/1‑148 (MQ=255)
                                                                           aCGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACAGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCCCTGTTTAATAAACGGTTGCGCCTCAATCGc                                                                      >  7:240324/1‑148 (MQ=255)
                                                                             gTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGTTTAATAAACGGTTGCGCCGCAATCGCCAg                                                                   >  5:125436/1‑149 (MQ=255)
                                                                                 cccGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGTTTAATAAACGGTTGCGCCGCAATCGCCAGctct                                                               <  3:159345/149‑1 (MQ=255)
                                                                                     gATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGTTTAATAAACGGTTGCGCCGCAATCGCCAGCTCTGTGc                                                           >  5:162742/1‑149 (MQ=255)
                                                                                           gccgccATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGTTTAATAAACGGTTGCGCCGCAATCGCCAGCTCTGTGCCCTGc                                                      >  5:194054/1‑148 (MQ=255)
                                                                                                    cTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGTTTAATAAACGGTTGCGCCGCAATCGCCAGCTCTGTGCCCTGCAATTAAATGa                                            <  2:137704/149‑1 (MQ=255)
                                                                                                       cGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGTTTAATAAACGGTTGCGCCGCAATCGCCAGCTCTGTGCCCTGCAATTAAATGaaa                                          <  1:262766/148‑1 (MQ=255)
                                                                                                        ggCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGTTTAATAAACGGTTGCGCCGCAATCGCCAGCTCTGTGCCCTGCAATTAAATGAAAAg                                        >  3:316067/1‑149 (MQ=255)
                                                                                                             tCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGTTTAATAAACGGTTGCGCCGCAATCGCCAGCTCTGTGCCCTGCAATTAAATGAAAAGCGCGt                                   >  3:65888/1‑149 (MQ=255)
                                                                                                                cGACGGATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGTTTAATAAACGGTTGCGCCGCAATCGCCAGCTCTGTGCCCTGCAATTAAATGAAAAGCGCGTCGa                                >  7:19018/1‑149 (MQ=255)
                                                                                                                 gACGGATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGTTTAATAAACGGTTGCGCCGCAATCGACCGCTCTGTGCCCTGCAATTAAATGAAAAGCGCGTCGAt                               >  1:31061/1‑149 (MQ=255)
                                                                                                                    ggATCAGCTGGATATGTACCGCCGTGTAAATCGGCTCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCAGCTGTTTAATAAACGGTTGCGCCGCAATCGCCAGCTCTGTGCCCTGCAATTAAATGAAAAGCGCGTCGAtgtg                            <  4:97204/149‑1 (MQ=255)
                                                                                                                     gATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGTTTAATAAACGGTTGCGCCGCAATCGCCAGCTCTGTGCCCTGCAATTAAATGAAAAGCGCGTCGATgtgg                           <  6:194054/149‑1 (MQ=255)
                                                                                                                      aTCAGCTTGTTATGTACCGCCGTTGAAATCGGCGCTCTGCTGCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGTTTAATAAACGGTTGCGCCGCAATCGCCAGCTCTGTGCCCTGCATTTAAATGAAAAGCGCGTCGATgtggt                          <  8:240324/149‑1 (MQ=255)
                                                                                                                         aGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTTTTTAATAAACGGTTGCGCCGCAATCGCCAGCTCTGTGCCCTGCAATTAACTGAAAAGCGCGTCGATGTGgtgt                        <  6:236879/148‑1 (MQ=255)
                                                                                                                         aGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGTTTAATAAACGGTTGCGCCGCAATCGCCAGCTCTGTGCCCTGCAATTAAATGAAAAGCGCGTCGATGTGGTGTa                       >  4:257343/1‑149 (MQ=255)
                                                                                                                           cTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGTTTAATATACGGATGCGCCGCAATCGCCAGCTCGGTGCCCTGCTATTAAATGAAAAGCTCGTCTATGTGGTGTATg                     >  2:170538/1‑149 (MQ=255)
                                                                                                                              gATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGTTTAATAAACGGTTGCGCCGCAATCGCCAGCTCTGTGCCCTGCAATTAAATGAAAAGCGCGTCGATGTGGTGTATGcc                   <  8:184974/148‑1 (MQ=255)
                                                                                                                              gATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGTTTAATAAACGGTTGCGCCGCAATCGCCAGCTCTGTGCCCTGCAATTAAATGAAAAGCGCGTCGATGTGGTGTATGcc                   >  8:121086/1‑148 (MQ=255)
                                                                                                                              gATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGTTTAATAAACGGTTGCGCCGCAATCGCCAGCTCTGTGCCCTGCAATTAAATGAAAAGCGCGTCGATGTGGTGTATGc                    <  2:55680/147‑1 (MQ=255)
                                                                                                                                 atGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGAt                                                                                                                        >  4:126769/1‑44 (MQ=255)
                                                                                                                                 atGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGAt                                                                                                                        <  3:126769/44‑1 (MQ=255)
                                                                                                                                 atGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGTTTAATAAACGGTTGCGCCGCAATCGCCAGCTCTGTGCCCTGCAATTAAATGAAAAGCGCGTCGATGTGGTGTATGCCa                  >  6:165947/1‑146 (MQ=255)
                                                                                                                                  tGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGTTTAATAAACGGTTGCGCCGCAATCGCCAGCTCTGTGCCCTGCAATTAAATGAAAAGCGCGTCGAtgtg                            <  5:18611/135‑1 (MQ=255)
                                                                                                                                  tGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGTTTAATAAACGGTTGCGCCGCAATCGCCAGCTCTGTGCCCTGCAATTAAATGAAAAGCGCGTCGAtgtg                            >  6:18611/1‑135 (MQ=255)
                                                                                                                                  tGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGTTTAATAAACGGTTGCGCCGCAATCGCCAGCTCTGTGCCCTGCAATTAAATGAAAAGCGCGTCGATGTGGTGTATGCCAAaga              <  3:257343/149‑1 (MQ=255)
                                                                                                                                    tACCGCCGGTGAAATCGGCGCTCTGCTCCTTGCAGCTCGATTAAGCAAAGAACAACAGCCGCTGTTTAATAAACGGTTGCTCCGCAATCGCCAGCTCTGTGCCCTGCAATTAAATGAAAAGCGCGTCGATGTGGTGTATGCCAAagag             <  3:170357/148‑1 (MQ=255)
                                                                                                                                     aCCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGTTTAATAAACGGTTGCGCCGCAATCGCCAGCTCTGTGCCCTGCAATTAAATGAAAAGCGCGTCGATGTGGTGTATGCCAAAGAGGt           >  8:107471/1‑149 (MQ=255)
                                                                                                                                     aCCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGTTTAATAAACGGTTGCGCCGCAATCGCCAGCTCTGTGCCCTGCAATTAAATGAAAAGCGCGTCGATGTGGTGTATGCCAAAGAGGt           >  6:117202/1‑149 (MQ=255)
                                                                                                                                     aCCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGTTTAATAAACGGTTGCGCCGCAATCGCCAGCTCTGTGCCCTGCAATTAAATGAAAAGCGCGTCGATGTGGTGTATGCCAAAGAGGt           >  4:21249/1‑149 (MQ=255)
                                                                                                                                             tGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGTTTAATAAACGGTTGCGCCGCAATCGCCAGCTCTGTGCCCTGCAATTAAATGAAAAGCGCGTCGATGTGGTGTATGCCAAAGAGGTGGAttttt   <  3:21249/149‑1 (MQ=255)
                                                                                                                                             tGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGTTTAATAAACGGTTGCGCCGCAATCGCCAGCTCTGTGCCCTGCAATTAAATGAAAAGCGCGTCGATGTGGTGTATGCCAAAGAGGTGGAttttt   <  5:117202/149‑1 (MQ=255)
                                                                                                                                               aaaTCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGTTTAATAAACGGTTGCGCCGCAATCGCCAGCTCTGTGCCCTGCAATTAAATGAAAAGCGCGTCAATGTGGTGTATGCCAAAGAGGTGGATTTTTc  <  2:57117/148‑1 (MQ=255)
                                                                                                                                               aaaTCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGTTTAATAAACGGTTGCGCCGCAATCGCCAGCTCTGTGCCCTGCAATTAAATGAAAAGCGCGGCGATGTGGTGTATGCCAAAGAGGTGGATTTTTc  <  5:165947/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
TCATGGTGAAAGTGACACCATGCCGCTCTATATCGGCAAAAGCGTTAACATCCGCAGCCGCGTCCTTTCTCATTTACGTACCCCGGATGAAGCCGCCATGCTACGGCAATCCCGACGGATCAGCTGGATATGTACCGCCGGTGAAATCGGCGCTCTGCTCCTTGAAGCGCGATTAATCAAAGAACAACAGCCGCTGTTTAATAAACGGTTGCGCCGCAATCGCCAGCTCTGTGCCCTGCAATTAAATGAAAAGCGCGTCGATGTGGTGTATGCCAAAGAGGTGGATTTTTC  >  NC_000913/1823634‑1823924

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