Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 210,017 C→T 5.6% I113I (ATC→ATT ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913210,0170CT 5.6% 282.7 / 4.5 89I113I (ATC→ATTldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base C (48/36);  new base T (3/2);  total (51/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

TCGATTTATGTAGCGATATCAATCAGCTTAATGAATATCTCCCGCTTTATGCCTTCATCAACACCCACTCGACGATGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTGTACGCCGGGGCATATGGGCGGCACCGCATATCAAAAAAGCCC  >  NC_000913/209869‑210163
                                                                                                                                                    |                                                                                                                                                  
tcGATTTATGTAGCGATATCAATCAGCTTAATGAATATCTCCCGCTTTATGCCTTCATCAACACCCACTCGACGATGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATc                                                                                                                                                    >  1:307616/1‑149 (MQ=255)
   aTTTATGTAGCGATATCAATCAGCTTAATGAATATCTCCCGCTTTATGCCTTCATCAACACCCACTCGACGATGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGcc                                                                                                                                                 >  2:777224/1‑149 (MQ=255)
       aTGTAGCGATATCAATCAGCTTAATGAATATCTCCCGCTTTATGCCTTCATCAACACCCACTCGACGATGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTc                                                                                                                                             <  3:357031/149‑1 (MQ=255)
               atatCAATCAGCTTAATGAATATCTCCCGCTTTATGCCTTCATCAACACCCACTCGACGATGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGt                                                                                                                                     <  1:128896/149‑1 (MQ=255)
                tatCAATCAGCTTAATGAATATCTCCCGCTTTATGCCTTCATCAACACCCACTCGACGATGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTc                                                                                                                                    <  1:295373/149‑1 (MQ=255)
                  tcaatcaGCTTAATGAATATCTCCCGCTTTATGCCTTCATCAACACCCACTCGACGATGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAg                                                                                                                                  >  4:810583/1‑149 (MQ=255)
                    aatcaGCTTAATGAATATCTCCCGCTTTATGCCTTCATCAACACCCACTCGACGATGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATTGCCATTCGTATGCGTCAGTa                                                                                                                                >  3:846548/1‑149 (MQ=255)
                    aatcaGCTTAATGAATATCTCCCGCTTTATGCCTTCATCAACACCCACTCGACGATGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATTGCCATTCGTATGCGTCAGTa                                                                                                                                >  1:333599/1‑149 (MQ=255)
                    aatcaGCTTAATGAATATCTCCCGCTTTATGCCTTCATCAACACCCACTCGACGATGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTa                                                                                                                                >  2:263808/1‑149 (MQ=255)
                     atcaGCTTAATGAATATCTCCCGCTTTATGCCTTCATCAACACCCACTCGACGATGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTac                                                                                                                               >  1:569308/1‑149 (MQ=255)
                     atcaGCTTAATGAATATCTCCCGCTTTATGCCTTCATCAACACCCACTCGACGATGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTac                                                                                                                               >  1:3863/1‑149 (MQ=255)
                           ttAATGAATATCTCCCGCTTTATGCCTTCATCAACACCCACTCGACGATGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACcgac                                                                                                                         <  2:795047/149‑1 (MQ=255)
                              aTGAATATCTCCCGCTTTATGCCTTCATCAACACCCACTCGACGATGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGaa                                                                                                                      >  1:796149/1‑149 (MQ=255)
                              aTGAATATCTCCCGCTTTATGCCTTCATCAACACCCACTCGACGATGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGaa                                                                                                                      <  4:656855/149‑1 (MQ=255)
                               tGAATATCTCCCGCGTTATGCCTTCGTCAACACCCGCTCGACGATGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAat                                                                                                                     <  3:663083/149‑1 (MQ=255)
                                gAATATCTCCCGCTTTATGCCTTCATCAACACCCACTCGACGATGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAata                                                                                                                    >  1:576877/1‑149 (MQ=255)
                                     tctcCCGCTTTATGCCTTCATCAACACCCACTCGACGATGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTg                                                                                                               >  2:423280/1‑149 (MQ=255)
                                      ctcCCGCTTTATGCCTTCATCAACACCCACTCGACGATGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGa                                                                                                              >  4:827076/1‑149 (MQ=255)
                                       tcCCGCTTTATGCCTTCATCAACACCCACTCGACGATGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGAt                                                                                                             >  4:444858/1‑149 (MQ=255)
                                       tcCCGCTTTATGCCTTCATCAACACCCACTCGACGATGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGAt                                                                                                             >  1:479044/1‑149 (MQ=255)
                                         ccGCTTTATGCCTTCATCAACACCCACTCGACGATGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATaa                                                                                                           >  4:571028/1‑149 (MQ=255)
                                         ccGCTTTATGCCTTCATCAACACCCACTCGACGATGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATaa                                                                                                           >  2:226595/1‑149 (MQ=255)
                                           gCTTTATGCCTTCATCAACACCCACTCGACGATGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACa                                                                                                         <  2:447272/149‑1 (MQ=255)
                                                  gCCTTCATCAACACCCACTCGACGATGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACAcc                                                                                                  >  4:271290/1‑149 (MQ=255)
                                                    cTTCATCAACACCCACTCGACGATGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACAccgc                                                                                                >  3:296086/1‑149 (MQ=255)
                                                    cTTCATCAACACCCACTCGACGATGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACAccgc                                                                                                >  3:182564/1‑149 (MQ=255)
                                                     ttCATCAACACCCACTCGACGATGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACAccgcc                                                                                               >  2:877259/1‑149 (MQ=255)
                                                     ttCATCAACACCCACTCGACGATGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCAATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACAccgcc                                                                                               >  3:772573/1‑149 (MQ=255)
                                                         tcaACACCCACTCGACGATGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTc                                                                                           >  4:641854/1‑149 (MQ=255)
                                                         tcaACACCCACTCGACGATGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTc                                                                                           >  3:440480/1‑149 (MQ=255)
                                                         tcaACACCCACTCGACGATGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTc                                                                                           >  2:647717/1‑149 (MQ=255)
                                                          caACACCCACTCGACGATGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCa                                                                                          >  4:268608/1‑149 (MQ=255)
                                                            acacCCACTCGACGATGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGTGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACg                                                                                        >  3:295681/1‑149 (MQ=255)
                                                            acacCCACTCGACGATGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACg                                                                                        >  3:747379/1‑149 (MQ=255)
                                                               cccACTCGACGATGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGaaa                                                                                     <  1:397317/149‑1 (MQ=255)
                                                                    tCGACGATGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCtt                                                                                <  1:178268/149‑1 (MQ=255)
                                                                    tCGACGATGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCtt                                                                                <  4:303728/149‑1 (MQ=255)
                                                                    tCGACGATGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCtt                                                                                <  3:641854/149‑1 (MQ=255)
                                                                      gacgaTGGATGTCAGCGTGCAGGATATGAGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGt                                                                              >  3:58987/1‑149 (MQ=255)
                                                                        cgaTGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGttt                                                                            <  2:281344/149‑1 (MQ=255)
                                                                         ggtggatgTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGaaa                                                                                     >  1:743968/3‑139 (MQ=255)
                                                                         gatggatgTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGaat                                                                                     <  2:743968/139‑2 (MQ=255)
                                                                         gatggatgTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTa                                                                           >  1:870911/1‑149 (MQ=255)
                                                                           tggatgTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATTGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTAcc                                                                         <  2:333599/149‑1 (MQ=255)
                                                                           tggatgTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATTGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTAcc                                                                         <  4:846548/149‑1 (MQ=255)
                                                                            ggatgTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAGCATTACACCGCCGTTCACGAAAGCCTTGTTTACCt                                                                        >  3:592103/1‑149 (MQ=255)
                                                                            ggatgTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCt                                                                        <  3:268608/149‑1 (MQ=255)
                                                                              atgTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTAc                                                                      <  1:647717/149‑1 (MQ=255)
                                                                              atgTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTAc                                                                      <  3:271290/149‑1 (MQ=255)
                                                                                gTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGt                                                                    >  2:747722/1‑149 (MQ=255)
                                                                                  cAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCa                                                                  >  2:256870/1‑149 (MQ=255)
                                                                                    gCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCaaa                                                                <  3:246382/149‑1 (MQ=255)
                                                                                        gCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTAAAAGAGc                                                            <  2:479044/149‑1 (MQ=255)
                                                                                         cAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCg                                                           >  4:87346/1‑149 (MQ=255)
                                                                                            gATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGaa                                                        <  1:777224/149‑1 (MQ=255)
                                                                                            gATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGaa                                                        <  4:131442/149‑1 (MQ=255)
                                                                                            gATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGaa                                                        <  2:3863/149‑1 (MQ=255)
                                                                                                  cGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTacac                                                  <  4:766696/149‑1 (MQ=255)
                                                                                                  cGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTacac                                                  <  4:223380/149‑1 (MQ=255)
                                                                                                       ggCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCtttt                                             <  2:796149/149‑1 (MQ=255)
                                                                                                       ggCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCtttt                                             <  1:216955/149‑1 (MQ=255)
                                                                                                         cgcTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTGt                                           <  1:423280/149‑1 (MQ=255)
                                                                                                         cgcTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTGt                                           <  1:263808/149‑1 (MQ=255)
                                                                                                          gcTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTGTa                                          >  2:192450/1‑149 (MQ=255)
                                                                                                          gcTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGGTCACGAAAGCCTTGCTTACCTACGTCAAAGAGCGGAAGTACACCTTTTGTa                                          >  2:633894/1‑149 (MQ=255)
                                                                                                           ctctGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCt                                                                        <  1:638848/118‑1 (MQ=255)
                                                                                                           ctctGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCt                                                                        >  2:638848/1‑118 (MQ=255)
                                                                                                            tctGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTacac                                                                                                   >  2:156102/1‑90 (MQ=255)
                                                                                                            tctGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCg                                                           <  1:215332/130‑1 (MQ=255)
                                                                                                            tctGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCg                                                           >  2:215332/1‑130 (MQ=255)
                                                                                                            tctGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTGTACg                                        >  3:182508/1‑149 (MQ=255)
                                                                                                            tctGGTTTTTTGAATATGCGCTGGAGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTacac                                                                                                   <  1:156102/90‑1 (MQ=255)
                                                                                                                gTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATTGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTGTACGCCgg                                    >  3:822605/1‑149 (MQ=255)
                                                                                                                 ttttttGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTGTACGCCggg                                   >  4:273502/1‑149 (MQ=255)
                                                                                                                   ttttGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTGTTCGCCGGGGc                                 >  2:116083/1‑149 (MQ=255)
                                                                                                                   ttttGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTGTACGCCGGGGc                                 >  3:365070/1‑149 (MQ=255)
                                                                                                                   ttttGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTGTACGCCGGGGc                                 >  3:159851/1‑149 (MQ=255)
                                                                                                                        aaTATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTGTACGCCGGGGCATATg                            >  1:123400/1‑149 (MQ=255)
                                                                                                                           atGCGCTGGGGCAGGCGTAAGATATCGCCATTCATCTGCGGCAATACACTCACGAATATCTTGATAACATTACACCGCCGTTCACGATAGGCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTCGTACGCCGGATCATATGggc                         <  4:564390/149‑1 (MQ=255)
                                                                                                                            tGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAGCATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTGTTCGCCGGGGCATATGggcg                        <  1:256870/149‑1 (MQ=255)
                                                                                                                             gcgcTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTGTACGCCGGGGCATATGggcgg                       >  2:218448/1‑149 (MQ=255)
                                                                                                                               gcTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTGTACGCCGGGGCATATGGGCGGCa                     >  1:782377/1‑149 (MQ=255)
                                                                                                                                  ggggCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTGTACGCCGGGGCATATGGGCGGCACCg                  <  4:182508/149‑1 (MQ=255)
                                                                                                                                     gCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTTATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGATGTACACCTTTTGTACGCCGGGGCATATGGGCGGCACCGCat               >  3:671281/1‑149 (MQ=255)
                                                                                                                                      cAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTGTACGCCGGGGCATATGGGCGGCACCGCata              <  3:248728/149‑1 (MQ=255)
                                                                                                                                         gCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTGTACGCCGGGGCATATGGGCGGCACCGCATATCa           >  1:55537/1‑149 (MQ=255)
                                                                                                                                           ggAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTGTACGCCGGGGCATATGGGCGGCACCGCATATCaaa         <  3:420643/149‑1 (MQ=255)
                                                                                                                                            gAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTGTACGCCGGGGCATATGGGCGGCACCGCATATCaaaa        >  4:358185/1‑149 (MQ=255)
                                                                                                                                                 tatCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTGTACGCCGGGGCATATGGGCGGCACCGCATATCAAAAAAGcc   <  3:621959/149‑1 (MQ=255)
                                                                                                                                                  atCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTGTACGCCGGGGCATATGGGCGGCACCGCATATCAAAAAAGccc  >  3:516928/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TCGATTTATGTAGCGATATCAATCAGCTTAATGAATATCTCCCGCTTTATGCCTTCATCAACACCCACTCGACGATGGATGTCAGCGTGCAGGATATGCGGATGGCGCTCTGGTTTTTTGAATATGCGCTGGGGCAGGCGGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTGTACGCCGGGGCATATGGGCGGCACCGCATATCAAAAAAGCCC  >  NC_000913/209869‑210163

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