Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,457,530 G→T 5.1% V12L (GTA→TTA) ‡ paaF → putative 2,3‑dehydroadipyl‑CoA hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,457,5300GT 5.1% 526.3 / 9.6 155V12L (GTA→TTA) ‡paaFputative 2,3‑dehydroadipyl‑CoA hydratase
Reads supporting (aligned to +/- strand):  ref base G (68/78);  new base T (2/6);  total (70/85)
Fisher's exact test for biased strand distribution p-value = 2.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

Only 100 of 156 total aligned reads displayed.

ATGGAAACCAATTACAGTCTGGAACCGGATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCTGTGTGATTACCGGTAATGCACGCTTTTTTGCCGC  >  NC_000913/1457381‑1457678
                                                                                                                                                     |                                                                                                                                                    
aTGGAAACCAATTACAGTCTGGAACCGGATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAg                                                                                                                                                      >  1:894103/1‑150 (MQ=255)
   gAAACCAAATACAGTCTGGAACCGGCTCAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGAGGTGGTGGGTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGGCCGCAACGAGTAt                                                                                                                                                   <  2:105644/150‑1 (MQ=255)
    aaaCCAATTACAGTCTGGAACCGGATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGCCCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTAtt                                                                                                                                                  >  2:1060189/1‑150 (MQ=255)
      aCCAATTACAGTCTGGAACCGGATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTAttgt                                                                                                                                                <  2:877006/150‑1 (MQ=255)
       ccAATTACAGTCTGGAACCGGATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTAttgtt                                                                                                                                               >  1:77181/1‑150 (MQ=255)
        cAATTACAGTCTGGAACCGGATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTAttgttg                                                                                                                                              <  1:158912/150‑1 (MQ=255)
        cAATTACAGTCTGGAACCGGATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGAACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTAttgttg                                                                                                                                              >  2:913316/1‑150 (MQ=255)
          aTTACAGTCTGGAACCGGATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCt                                                                                                                                            <  2:1155991/150‑1 (MQ=255)
            tACAGTCTGGAACCGGATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGa                                                                                                                                          <  1:78584/150‑1 (MQ=255)
             aCAGTCTGGAACCGGATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGAc                                                                                                                                         <  2:1176529/150‑1 (MQ=255)
             aCAGTCTGGAACCGGATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACAATACTCGTTGCTGAc                                                                                                                                         <  1:2495/150‑1 (MQ=255)
             aCAGTCTGGAACCGGATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGAGGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACAATACTCGTTGCTGAc                                                                                                                                         <  1:1081498/150‑1 (MQ=255)
              cAGTCTGGAACCGGATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTAGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACAATACTCGTTGCTGAcg                                                                                                                                        <  2:274542/150‑2 (MQ=255)
               aGTCTGGAACCGGATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGAccc                                                                                                                                       <  1:110969/150‑1 (MQ=255)
               aGTCTGGAACCGGATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACAATACTCGTTGCTGAcgg                                                                                                                                       <  1:483272/150‑3 (MQ=255)
                   tGGAACCGGATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTaa                                                                                                                                   <  1:239626/150‑1 (MQ=255)
                      aaCCGGATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCg                                                                                                                                <  2:1174560/150‑1 (MQ=255)
                       aCCGGATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGt                                                                                                                               <  2:21789/150‑1 (MQ=255)
                        ccGGATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTc                                                                                                                              <  1:1197950/150‑1 (MQ=255)
                         cGGATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTcc                                                                                                                             >  1:645468/1‑150 (MQ=255)
                          ggATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTccc                                                                                                                            >  1:172337/1‑150 (MQ=255)
                          ggATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTccc                                                                                                                            >  2:772161/1‑150 (MQ=255)
                             tGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCccgcc                                                                                                                         <  2:9398/150‑1 (MQ=255)
                             tGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCccgcc                                                                                                                         <  2:1372618/150‑1 (MQ=255)
                             ggAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCccgcc                                                                                                                         <  2:150163/149‑1 (MQ=255)
                                aCTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCa                                                                                                                      <  2:1363369/150‑1 (MQ=255)
                                aCTGGCCGCAGGTTATGTGTTGAGTTGACAGGCACTGCCGCTGACCAGCGATGTGGTGCTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACCAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCa                                                                                                                      >  2:1164527/1‑150 (MQ=255)
                                 cTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCAc                                                                                                                     >  2:1356636/1‑150 (MQ=255)
                                 cTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCAc                                                                                                                     >  1:935484/1‑150 (MQ=255)
                                   ggCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGt                                                                                                                   <  2:1271556/150‑1 (MQ=255)
                                       gCAGGTTATGTGTTGAGTTGCCAAGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATg                                                                                                               <  2:173666/150‑1 (MQ=255)
                                          ggTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGCCCCTTACCCGTCCCGCCGCACGTAATgcgc                                                                                                            >  1:183765/1‑150 (MQ=255)
                                          ggTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATgcgc                                                                                                            >  2:664753/1‑150 (MQ=255)
                                             tATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTaa                                                                                                         <  1:1027600/150‑1 (MQ=255)
                                              aTGTCTTGAGTTGCCAGGCACTGCCGCTGACCAGCGAAGTGGTCGTTGACATTGACCCGAAGGGAATGGCATCAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTTCTTACCCTTAACCGTCCCGCCGCACGTAATGCGCTaaa                                                                                                        >  2:271372/1‑150 (MQ=255)
                                                gtgtTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAaata                                                                                                      >  1:4141/1‑150 (MQ=255)
                                                 tgtTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAAc                                                                                                                                  >  1:395001/1‑121 (MQ=255)
                                                 tgtTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAaataa                                                                                                     >  1:1419707/1‑150 (MQ=255)
                                                  gtTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAaataat                                                                                                    >  1:761744/1‑150 (MQ=255)
                                                    tGAGTTGCCAGGCACTGCCGCTGCCCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGc                                                                                                  >  2:30100/1‑150 (MQ=255)
                                                    tGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGc                                                                                                  >  1:291181/1‑150 (MQ=255)
                                                    tGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGc                                                                                                  >  1:404594/1‑150 (MQ=255)
                                                     gAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGcc                                                                                                 >  1:131580/1‑150 (MQ=255)
                                                         tGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGATCTGATCGTCAGCCGTCAGCAAGGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCctgc                                                                                             <  2:1305719/150‑1 (MQ=255)
                                                         tGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCctgc                                                                                             <  1:127406/150‑1 (MQ=255)
                                                          gCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCctgct                                                                                            <  2:163187/150‑1 (MQ=255)
                                                            cAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGa                                                                                          >  2:1191633/1‑150 (MQ=255)
                                                              ggCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATg                                                                                        <  1:427168/150‑1 (MQ=255)
                                                               gCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGc                                                                                       <  2:845389/150‑1 (MQ=255)
                                                                    gCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTg                                                                                  <  1:1339988/150‑1 (MQ=255)
                                                                    gCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTg                                                                                  >  1:848885/1‑150 (MQ=255)
                                                                     ccGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTgg                                                                                 >  1:987326/1‑150 (MQ=255)
                                                                      cGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGt                                                                                <  1:669565/150‑1 (MQ=255)
                                                                      cGCTGACCAGCGATGTGGGGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGt                                                                                <  2:1221335/150‑1 (MQ=255)
                                                                      cGCTGACCAGCGAGGTGGGGGTTGACTTTGACGCGAAGGGGAGGGCATGAGCGAACTGATCGGCAGCCGTCAGCAACGAGTATTGTGGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGt                                                                                <  2:567197/150‑1 (MQ=255)
                                                                         tGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTaaa                                                                             <  2:356290/150‑1 (MQ=255)
                                                                          gACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAAt                                                                            <  1:905304/150‑1 (MQ=255)
                                                                          gACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAAt                                                                            >  1:1226482/1‑150 (MQ=255)
                                                                          gACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAAt                                                                            <  2:913647/150‑1 (MQ=255)
                                                                           aCCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATg                                                                           <  1:993464/150‑1 (MQ=255)
                                                                            ccAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGa                                                                          <  1:945633/150‑1 (MQ=255)
                                                                            ccAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGa                                                                          <  2:483190/150‑1 (MQ=255)
                                                                                  aTGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGa                                                                    <  2:916809/150‑1 (MQ=255)
                                                                                    gtggtgGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAAt                                                                                                                >  2:1043790/1‑104 (MQ=255)
                                                                                     tggtggTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGACGc                                                                 <  1:1060526/150‑1 (MQ=255)
                                                                                     tggtggTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGc                                                                 >  1:836093/1‑150 (MQ=255)
                                                                                      ggtggtTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCt                                                                >  2:1221508/1‑150 (MQ=255)
                                                                                       gtggtTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTg                                                               <  2:213793/150‑1 (MQ=255)
                                                                                         ggtTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCg                                                             <  2:225581/150‑1 (MQ=255)
                                                                                          gtTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCgg                                                            >  1:862824/1‑150 (MQ=255)
                                                                                          gtTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCgg                                                            <  2:404594/150‑1 (MQ=255)
                                                                                               cTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTAcc                                                       >  1:1035009/1‑150 (MQ=255)
                                                                                               cTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTAcc                                                       >  2:395143/1‑150 (MQ=255)
                                                                                               cTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTAcc                                                       >  1:646013/1‑150 (MQ=255)
                                                                                                tttGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCg                                                      <  1:1096900/150‑1 (MQ=255)
                                                                                                tttGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCg                                                      >  1:1303086/1‑150 (MQ=255)
                                                                                                     cgcgAAGGGGAGGGCATCACGGAACTGATCGTCAGCCGTCATCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCCCCTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATAcc                                                 <  2:482764/150‑1 (MQ=255)
                                                                                                        gAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGc                                              >  1:979590/1‑150 (MQ=255)
                                                                                                          aGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCAt                                            >  2:961196/1‑150 (MQ=255)
                                                                                                             ggATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTc                                         <  1:1167595/150‑1 (MQ=255)
                                                                                                               aTGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCgg                                       <  2:646013/150‑1 (MQ=255)
                                                                                                                tGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGt                                      >  2:466065/1‑150 (MQ=255)
                                                                                                                     tGAGCGAACTGATCGTCAGCCGGCAGCAACGAGTATTGTTGCTGACCCTAAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCtgtg                                 <  2:1408073/150‑1 (MQ=255)
                                                                                                                       aGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCtgtgtg                               >  2:522904/1‑150 (MQ=255)
                                                                                                                        gCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTAcc                                                       >  1:498306/1‑125 (MQ=255)
                                                                                                                             cTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCTGTGTGATTAcc                         <  1:961196/150‑1 (MQ=255)
                                                                                                                              tGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCTGTGTGATTACCg                        >  1:1362953/1‑150 (MQ=255)
                                                                                                                              tGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCTGTGTGATTACCg                        >  1:1182771/1‑150 (MQ=255)
                                                                                                                              tGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCTGTGTGATTACCg                        >  2:1082514/1‑150 (MQ=255)
                                                                                                                                   gtcagccgtcagcAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCTGTGTGATTACCGGTAAt                   <  2:195933/150‑1 (MQ=255)
                                                                                                                                   gtcagccgtcagcAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCTGTGTGATTACCGGTAAt                   >  1:309739/1‑150 (MQ=255)
                                                                                                                                       gccgtcagcAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCTGTGTGATTACCGGTAATGCAc               >  2:361508/1‑150 (MQ=255)
                                                                                                                                         ggtcagcAACAATACTCGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCTGTGTGATTACCGGTAATGCACGc             >  1:256420/2‑150 (MQ=255)
                                                                                                                                         ggtcagcAACAATACTCGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATATCGGTCTGTGTGATTACCGGTAATGCACGc             >  1:923435/2‑150 (MQ=255)
                                                                                                                                            cagcAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCTGTGTGATTACCGGTAATGCACGCttt          <  2:1068688/150‑1 (MQ=255)
                                                                                                                                               cAACGAGTATTGTTGCTGACCCTTAACCGTGCCGCCGCACGTAATGCGCTAAATAATGCCCTTCTTATTCAACTTGTAAATGAACTGGTAGCTGCGGCTACCGATTCCAGCATTTCGGTCTGTGTGATTACCGGTAATGCACGCtttttt       <  1:1103432/150‑1 (MQ=255)
                                                                                                                                                aaCGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCTGTGTGATTACCGGTAATGCa                <  2:1348651/140‑1 (MQ=255)
                                                                                                                                                aaCGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCTGTGTGATTACCGGTAATGCACGCTTTTTTg      <  2:309739/150‑1 (MQ=255)
                                                                                                                                                 aCGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCTGTGTGATTACCGGTAATGCACGCTTTTTTGc     <  2:952515/150‑1 (MQ=255)
                                                                                                                                                    aGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGGCTGGGTGATTACCGGTAATGCACGCTTTTTTGCCGc  <  2:979590/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATGGAAACCAATTACAGTCTGGAACCGGATGAACTGGCCGCAGGTTATGTGTTGAGTTGCCAGGCACTGCCGCTGACCAGCGATGTGGTGGTTGACTTTGACGCGAAGGGGATGGCATGAGCGAACTGATCGTCAGCCGTCAGCAACGAGTATTGTTGCTGACCCTTAACCGTCCCGCCGCACGTAATGCGCTAAATAATGCCCTGCTGATGCAACTGGTAAATGAACTGGAAGCTGCGGCTACCGATACCAGCATTTCGGTCTGTGTGATTACCGGTAATGCACGCTTTTTTGCCGC  >  NC_000913/1457381‑1457678

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