Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,681,990 C→A 5.4% K215N (AAG→AAT fumA ← fumarate hydratase (fumarase A), aerobic Class I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,681,9900CA 5.4% 789.4 / 13.1 204K215N (AAG→AATfumAfumarate hydratase (fumarase A), aerobic Class I
Reads supporting (aligned to +/- strand):  ref base C (98/95);  new base A (3/8);  total (101/103)
Fisher's exact test for biased strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 217 total aligned reads displayed.

ACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGC  >  CP009273/1681843‑1682137
                                                                                                                                                   |                                                                                                                                                   
 ccGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCAtct                                                                                                                                                    >  1:1144385/1‑148 (MQ=255)
 ccGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCAtct                                                                                                                                                    >  1:2496225/1‑148 (MQ=255)
    tGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCAtcttct                                                                                                                                                 <  1:824872/148‑1 (MQ=255)
    tGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCAtcttct                                                                                                                                                 >  1:2296069/1‑148 (MQ=255)
     gCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTtctc                                                                                                                                                >  1:773739/1‑148 (MQ=255)
      cTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCa                                                                                                                                               <  2:894767/148‑1 (MQ=255)
          ttCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCa                                                                                                                                           >  2:78880/1‑148 (MQ=255)
            cccTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAgg                                                                                                                                         <  2:3409050/148‑1 (MQ=255)
            cccTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAgg                                                                                                                                         <  2:3046783/148‑1 (MQ=255)
            cccTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGAACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAgg                                                                                                                                         >  2:2740160/1‑148 (MQ=255)
             ccTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGt                                                                                                                                        >  1:3412153/1‑148 (MQ=255)
               ttCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTaa                                                                                                                                      <  2:283931/148‑1 (MQ=255)
                tCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAAt                                                                                                                                     >  1:2669350/1‑148 (MQ=255)
                  cGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAAttt                                                                                                                                   >  1:699265/1‑148 (MQ=255)
                     tGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTc                                                                                                                                >  2:1606665/1‑148 (MQ=255)
                       gCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAg                                                                                                                              <  2:1031575/148‑1 (MQ=255)
                       gCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGACGCCGTACCCAGCGAGCGCATCTTATCAACCAGGTAATTTTTCAg                                                                                                                              >  2:3275637/1‑148 (MQ=255)
                          gTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGttt                                                                                                                           >  1:1694674/1‑148 (MQ=255)
                          gTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGttt                                                                                                                           >  1:1695960/1‑148 (MQ=255)
                                   aGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGt                                                                                                                  >  2:1429569/1‑148 (MQ=255)
                                        ttCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTa                                                                                                             >  1:3254417/1‑148 (MQ=255)
                                        ttCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTa                                                                                                             <  1:1748650/148‑1 (MQ=255)
                                            cgGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGc                                                                                                         >  2:1656849/1‑148 (MQ=255)
                                             ggAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCt                                                                                                        <  1:1428777/148‑1 (MQ=255)
                                               aaGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTACCGCttt                                                                                                      >  2:2771635/1‑148 (MQ=255)
                                                 gCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTgg                                                                                                    >  1:2170210/1‑148 (MQ=255)
                                                 gCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTgg                                                                                                    <  1:1464401/148‑1 (MQ=255)
                                                 gCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTgg                                                                                                    <  1:1464373/148‑1 (MQ=255)
                                                  ccAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGt                                                                                                   >  2:1345525/1‑148 (MQ=255)
                                                   cAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGtt                                                                                                  >  2:640692/1‑148 (MQ=255)
                                                     gTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTc                                                                                                <  2:1484708/148‑1 (MQ=255)
                                                     gTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTc                                                                                                <  1:135850/148‑1 (MQ=255)
                                                      tttCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTcc                                                                                               <  2:2495215/148‑1 (MQ=255)
                                                      tttCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTcc                                                                                               >  1:1029289/1‑148 (MQ=255)
                                                      tttCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTcc                                                                                               >  1:2735471/1‑148 (MQ=255)
                                                        tCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCTGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTg                                                                                             <  2:1715704/148‑1 (MQ=255)
                                                          aCCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGAt                                                                                           >  2:2885563/1‑148 (MQ=255)
                                                            cGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATag                                                                                         <  1:1966589/148‑1 (MQ=255)
                                                              ttttAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATagag                                                                                       <  1:1365922/148‑1 (MQ=255)
                                                              ttttAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATagag                                                                                       >  1:1101631/1‑148 (MQ=255)
                                                               tttAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAgaga                                                                                      <  2:1043388/148‑1 (MQ=255)
                                                                 tAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATa                                                                                    <  2:1472809/148‑1 (MQ=255)
                                                                   aGGTTCGTGTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACg                                                                                  <  2:3269009/148‑1 (MQ=255)
                                                                    ggTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGt                                                                                 <  2:2511797/148‑1 (MQ=255)
                                                                     gTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTc                                                                                >  2:375568/1‑148 (MQ=255)
                                                                     gTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTGTGGTTTCCTGATAGAGATACGTc                                                                                >  2:497471/1‑148 (MQ=255)
                                                                      ttCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCGTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCt                                                                               <  2:2970849/148‑1 (MQ=255)
                                                                      ttCGTTTCAGCAGAAGTTCCACCAATAACGAACGCAAAATTATACGGAGGACAGGCCGCCGTACCCAGCTTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTACTGATAGAGATACGTCt                                                                               >  1:892230/1‑148 (MQ=255)
                                                                        cGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCACCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTg                                                                             >  2:3555654/1‑148 (MQ=255)
                                                                         gTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGt                                                                            <  1:2983639/148‑1 (MQ=255)
                                                                         gTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGt                                                                            <  1:2457473/148‑1 (MQ=255)
                                                                            tCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTg                                                                         <  1:2138037/148‑1 (MQ=255)
                                                                            tCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTg                                                                         <  1:1318845/148‑1 (MQ=255)
                                                                               gCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCg                                                                      >  1:3463234/1‑148 (MQ=255)
                                                                               gCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCg                                                                      >  1:285525/1‑148 (MQ=255)
                                                                                  gAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAac                                                                   <  2:2460276/148‑1 (MQ=255)
                                                                                   aaGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAacc                                                                  >  1:665898/1‑148 (MQ=255)
                                                                                    aGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGGGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAacca                                                                 <  1:3046903/148‑1 (MQ=255)
                                                                                     gTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATATTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAaccac                                                                >  2:1295058/1‑148 (MQ=255)
                                                                                       tCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAaccacca                                                              >  2:3393052/1‑148 (MQ=255)
                                                                                         caccaATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAaccaccacc                                                            >  2:561454/1‑148 (MQ=255)
                                                                                         caccaATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAaccaccacc                                                            >  2:1868848/1‑148 (MQ=255)
                                                                                            caATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCttt                                                         >  2:2452592/1‑148 (MQ=255)
                                                                                                aacgaacgCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATATTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCg                                                     <  1:1599547/148‑1 (MQ=255)
                                                                                                   gaacgCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATa                                                  <  2:1144385/148‑1 (MQ=255)
                                                                                                   gaacgCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATa                                                  <  2:11595/148‑1 (MQ=255)
                                                                                                     acgCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACa                                                >  1:426974/1‑148 (MQ=255)
                                                                                                       gCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACaga                                              <  2:892230/148‑1 (MQ=255)
                                                                                                         aaTATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAgg                                            <  1:3165448/148‑1 (MQ=255)
                                                                                                          atatGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATATTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGa                                           <  2:3024641/148‑1 (MQ=255)
                                                                                                            atGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAAt                                         >  2:1247433/1‑148 (MQ=255)
                                                                                                            atGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAAt                                         >  1:1061325/1‑148 (MQ=255)
                                                                                                             tGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAAtt                                        <  1:2197339/148‑1 (MQ=255)
                                                                                                             tGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAAtt                                        <  2:421313/148‑1 (MQ=255)
                                                                                                                  cGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTAt                                   <  2:2296069/148‑1 (MQ=255)
                                                                                                                  cGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTAt                                   <  1:9302/148‑1 (MQ=255)
                                                                                                                       gaCAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATtcgtc                              >  1:3541405/1‑148 (MQ=255)
                                                                                                                       gaCAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATtcgtc                              >  1:1219237/1‑148 (MQ=255)
                                                                                                                           ggCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCa                          >  1:2307749/1‑148 (MQ=255)
                                                                                                                           ggCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCa                          <  1:587994/148‑1 (MQ=255)
                                                                                                                            gccgccGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCAt                         >  2:3107814/1‑148 (MQ=255)
                                                                                                                             ccgccgTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATc                        >  1:1586463/1‑148 (MQ=255)
                                                                                                                                ccgTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAAc                     >  1:1287657/1‑148 (MQ=255)
                                                                                                                                ccgTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAAc                     <  1:2152362/148‑1 (MQ=255)
                                                                                                                                ccgTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAAc                     >  2:1599814/1‑148 (MQ=255)
                                                                                                                                  gTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACgg                   <  1:2943112/148‑1 (MQ=255)
                                                                                                                                    aCCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCa                 >  1:2434407/1‑148 (MQ=255)
                                                                                                                                         gCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAg            <  2:2735471/148‑1 (MQ=255)
                                                                                                                                         gCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAg            <  1:2744778/148‑1 (MQ=255)
                                                                                                                                         gCGTGCGCATATTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAg            <  2:23112/148‑1 (MQ=255)
                                                                                                                                          cGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAga           <  1:1877431/148‑1 (MQ=255)
                                                                                                                                          cGTGCGCATATTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAga           <  1:954126/148‑1 (MQ=255)
                                                                                                                                           gTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAgat          <  1:2554271/148‑1 (MQ=255)
                                                                                                                                           gTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAgat          <  1:514705/148‑1 (MQ=255)
                                                                                                                                            tGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAgatc         >  1:3238692/1‑148 (MQ=255)
                                                                                                                                             gcgcATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAgatcg        <  2:1992618/148‑1 (MQ=255)
                                                                                                                                              cgcATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAgatcga       <  2:266370/148‑1 (MQ=255)
                                                                                                                                               gcATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAgatcgat      <  1:1692852/148‑1 (MQ=255)
                                                                                                                                                 aTCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCt    >  1:2756758/1‑148 (MQ=255)
                                                                                                                                                  tcttctCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACAACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTg   >  1:286326/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGC  >  CP009273/1681843‑1682137

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