Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,462,029 T→C 6.2% G379G (GGT→GGC) ‡ paaJ → 3‑oxoadipyl‑CoA/3‑oxo‑5,6‑dehydrosuberyl‑CoA thiolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,462,0290TC 6.2% 375.6 / 9.5 113G379G (GGT→GGC) ‡paaJ3‑oxoadipyl‑CoA/3‑oxo‑5,6‑dehydrosuberyl‑CoA thiolase
Reads supporting (aligned to +/- strand):  ref base T (55/51);  new base C (3/4);  total (58/55)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 113 total aligned reads displayed.

TTGGGTGTACTACGCGAATTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCGATTGAAACCGCGTCTGTTGATGAGT  >  NC_000913/1461883‑1462176
                                                                                                                                                  |                                                                                                                                                   
      gtACTACGCGAATTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGACCGTTACGc                                                                                                                                             <  3:276341/149‑1 (MQ=255)
         ctccgcgAATTGGGACTGCCTGATGATGCCCCACATGTTAACCCCGACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTACCCGCCTGGCACTGGCTGCCAGCCGTGAGCTGCAACGGCGTAACGACCGTTACGccga                                                                                                                                          <  4:543781/146‑4 (MQ=255)
         ctacGCGAATTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGACCGTTACGccga                                                                                                                                          <  2:587594/149‑4 (MQ=255)
             gcgAATTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGc                                                                                                                                      >  4:215428/1‑149 (MQ=255)
             gcgAATTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGc                                                                                                                                      >  2:513685/1‑149 (MQ=255)
             gcgAATTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGACTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGc                                                                                                                                      >  2:742130/1‑149 (MQ=255)
             gcgAATTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCACGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGc                                                                                                                                      >  1:638648/1‑149 (MQ=255)
               gAATTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCAc                                                                                                                                    >  1:682986/1‑149 (MQ=255)
                  ttGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTTTCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCAt                                                                                                                                 >  4:202630/1‑149 (MQ=255)
                    ggggCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCAtgt                                                                                                                               <  3:225817/149‑1 (MQ=255)
                    ggggCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCAtgt                                                                                                                               <  4:409197/149‑1 (MQ=255)
                    ggggCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCAtgt                                                                                                                               <  2:412582/149‑1 (MQ=255)
                     gggCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCAtgtg                                                                                                                              >  4:770848/1‑149 (MQ=255)
                       gCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCa                                                                                                                            >  4:74164/1‑149 (MQ=255)
                         tGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATc                                                                                                                          <  1:615994/149‑1 (MQ=255)
                            cTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGgt                                                                                                                       >  1:146628/1‑149 (MQ=255)
                             tgatgaTGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGgtg                                                                                                                      >  1:505285/1‑149 (MQ=255)
                              gatgatGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGgtgt                                                                                                                     >  3:244951/1‑149 (MQ=255)
                               atgatgCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTc                                                                                                                                                    <  2:815367/117‑1 (MQ=255)
                               atgatgCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTc                                                                                                                                                    >  1:815367/1‑117 (MQ=255)
                                      cccACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAgg                                                                                                             <  1:305378/149‑1 (MQ=255)
                                      cccACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCAGTGTGCACCATGTGCATCGGTGTCGGTCatg                                                                                                             >  3:527485/1‑147 (MQ=255)
                                      cccACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGCTGCTCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAgg                                                                                                             <  2:638648/149‑1 (MQ=255)
                                       ccACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAggg                                                                                                            >  2:101707/1‑149 (MQ=255)
                                       ccACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAggg                                                                                                            >  1:856756/1‑149 (MQ=255)
                                             gTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGc                                                                                                      >  3:394553/1‑149 (MQ=255)
                                              ttAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGcc                                                                                                     >  3:131364/1‑149 (MQ=255)
                                              ttAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGcc                                                                                                     >  1:420219/1‑149 (MQ=255)
                                                aaCCCCAACGTAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGTCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCAt                                                                                                   >  2:826373/1‑149 (MQ=255)
                                                aaCCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCAt                                                                                                   >  2:518691/1‑149 (MQ=255)
                                                  ccccCACGGAGGCGCTATCGCCTTAAGCCATCCGTTGGGAATGAGTGGTGCCCCCCTGGCACTGGCTTCGCGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCACTGTGCACCGTGAGCATCGATGTCGGTCAGGGCACCGACATGa                                                                                                 <  2:569299/149‑1 (MQ=255)
                                                  ccccAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGa                                                                                                 >  3:755286/1‑149 (MQ=255)
                                                   cccAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGAt                                                                                                >  4:226240/1‑149 (MQ=255)
                                                   cccAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGAt                                                                                                >  2:554584/1‑149 (MQ=255)
                                                    ccAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGAtt                                                                                               >  1:681011/1‑149 (MQ=255)
                                                    ccAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGAtt                                                                                               <  3:2384/149‑1 (MQ=255)
                                                    ccAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGAtt                                                                                               >  3:315389/1‑149 (MQ=255)
                                                    ccAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGAtt                                                                                               >  4:300691/1‑149 (MQ=255)
                                                     cAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTc                                                                                              >  1:672454/1‑149 (MQ=255)
                                                      aaCGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCt                                                                                             >  3:300293/1‑149 (MQ=255)
                                                       aCGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTg                                                                                            <  4:363217/149‑1 (MQ=255)
                                                             gcgcTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCgt                                                                                      >  1:251840/1‑149 (MQ=255)
                                                              cgcTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCgtg                                                                                     <  4:251564/149‑1 (MQ=255)
                                                              cgcTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCgtg                                                                                     >  3:113880/1‑149 (MQ=255)
                                                                 tATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGttt                                                                                  >  1:137686/1‑149 (MQ=255)
                                                                   tCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGa                                                                                <  4:345294/149‑1 (MQ=255)
                                                                   tCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGa                                                                                <  2:572200/149‑1 (MQ=255)
                                                                      ccTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCa                                                                             >  3:160482/1‑149 (MQ=255)
                                                                              cATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAAcc                                                                     <  4:113880/149‑1 (MQ=255)
                                                                               aTCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCt                                                                    >  1:455826/1‑149 (MQ=255)
                                                                               aTCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCt                                                                    <  2:455826/149‑1 (MQ=255)
                                                                                 ccGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGc                                                                  <  1:518691/149‑1 (MQ=255)
                                                                                 ccGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGc                                                                  <  2:109462/149‑1 (MQ=255)
                                                                                 ccGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGc                                                                  <  1:826373/149‑1 (MQ=255)
                                                                                   gTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGa                                                                <  4:837674/149‑1 (MQ=255)
                                                                                    ttGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAg                                                               >  2:778945/1‑149 (MQ=255)
                                                                                    ttGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAg                                                               <  2:672454/149‑1 (MQ=255)
                                                                                    ttGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAg                                                               >  3:686387/1‑149 (MQ=255)
                                                                                      gggAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTa                                                             >  4:414272/1‑149 (MQ=255)
                                                                                        gAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTAcc                                                           >  4:104497/1‑149 (MQ=255)
                                                                                         aaTGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTAccc                                                          <  3:226240/149‑1 (MQ=255)
                                                                                         aaTGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTAccc                                                          <  1:15079/149‑1 (MQ=255)
                                                                                           tGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAgg                                                                                                             <  2:531630/96‑1 (MQ=255)
                                                                                           tGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAgg                                                                                                             >  1:531630/1‑96 (MQ=255)
                                                                                           tGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAgg                                                                                                             >  3:331996/1‑96 (MQ=255)
                                                                                           tGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAgg                                                                                                             <  4:331996/96‑1 (MQ=255)
                                                                                            gAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTAc                                                       >  3:431300/1‑149 (MQ=255)
                                                                                             aGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACa                                                      <  3:224931/149‑1 (MQ=255)
                                                                                               tggtgCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAAt                                                    <  3:840563/149‑1 (MQ=255)
                                                                                                    cccGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATaa                                               <  2:206327/149‑1 (MQ=255)
                                                                                                     ccGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAAc                                              <  3:74164/149‑1 (MQ=255)
                                                                                                        ccTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGa                                                                                >  2:215531/1‑112 (MQ=255)
                                                                                                        ccTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGa                                                                                <  1:215531/112‑1 (MQ=255)
                                                                                                        ccTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCaa                                           <  3:770848/149‑1 (MQ=255)
                                                                                                         cTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAAt                                          >  4:802523/1‑149 (MQ=255)
                                                                                                               cTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAg                                    <  3:802523/149‑1 (MQ=255)
                                                                                                                    tGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGa                               <  1:101707/149‑1 (MQ=255)
                                                                                                                    tGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGa                               >  4:381409/1‑149 (MQ=255)
                                                                                                                      ccagccaTGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGAcc                             >  4:219856/1‑149 (MQ=255)
                                                                                                                         gccaTGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCGa                          <  2:363366/149‑1 (MQ=255)
                                                                                                                           caTGAGCTGCAACGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTAGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCGAtt                        >  4:606308/1‑149 (MQ=255)
                                                                                                                            aTGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCGATTg                       <  1:327799/149‑1 (MQ=255)
                                                                                                                               aGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGGCCCGATTGaaa                    <  2:816261/149‑1 (MQ=255)
                                                                                                                               aGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGGCCCGATTGaaa                    >  1:816261/1‑149 (MQ=255)
                                                                                                                               aGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCGATTGaaa                    <  1:742130/149‑1 (MQ=255)
                                                                                                                                   gCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCGATTGAAACcgc                <  1:59499/149‑1 (MQ=255)
                                                                                                                                   gCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCGATTGAAACcgc                <  4:321085/149‑1 (MQ=255)
                                                                                                                                      aatgCGTAACGACCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCGAtt                        >  2:119747/4‑138 (MQ=255)
                                                                                                                                      aatgCGTAACGACCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCGAtt                        <  1:119747/135‑1 (MQ=255)
                                                                                                                                      aatgCGTAACGACCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCGATTGAAACCGCGTc             >  3:598655/4‑149 (MQ=255)
                                                                                                                                       cGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCGATTGAAACCGCGTCt            <  1:164352/149‑1 (MQ=255)
                                                                                                                                       atgCGTAACGACCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCGATTGAAACCGCGTCt            >  4:185499/3‑149 (MQ=255)
                                                                                                                                         gCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTAc                                                            <  2:582052/99‑1 (MQ=255)
                                                                                                                                         gCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTAc                                                            >  1:582052/1‑99 (MQ=255)
                                                                                                                                         gCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCGATTGAAACCGCGTCTGt          >  1:132743/1‑149 (MQ=255)
                                                                                                                                           gTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCGATTGAAACCGCGTCTGTtg        >  2:605514/1‑149 (MQ=255)
                                                                                                                                               cGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCTTATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCGATTGAAACCGCGTCTGTtgatga    <  4:755286/149‑1 (MQ=255)
                                                                                                                                               cGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCGATTGAAACCGCGTCTGTtgatga    <  3:390604/149‑1 (MQ=255)
                                                                                                                                                ggTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCGATTGAAACCGCGTCTGTTGATGAg   >  3:711724/1‑149 (MQ=255)
                                                                                                                                                 gTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCGATTGAAACCGCGTCTGTTGATGAGt  <  3:219856/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TTGGGTGTACTACGCGAATTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCGATTGAAACCGCGTCTGTTGATGAGT  >  NC_000913/1461883‑1462176

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