Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 250,909 A→T 100% L432Q (CTG→CAG)  pepD ← aminoacyl‑histidine dipeptidase (peptidase D)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273250,9090AT96.6% 592.6 / 5.4 192L432Q (CTG→CAG) pepDaminoacyl‑histidine dipeptidase (peptidase D)
Reads supporting (aligned to +/- strand):  ref base A (2/5);  new base T (94/88);  total (96/96)
Fisher's exact test for biased strand distribution p-value = 2.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 192 total aligned reads displayed.

TTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTT  >  CP009273/250764‑251053
                                                                                                                                                 |                                                                                                                                                
tttcAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTgg                                                                                                                                                >  2:1935800/1‑148 (MQ=255)
 ttcAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGTATACGGTTTTTTGAACAGACCACATTCCGGGc                                                                                                                                               <  1:499107/148‑1 (MQ=255)
 ttcAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGc                                                                                                                                               >  2:745869/1‑148 (MQ=255)
 ttcAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGc                                                                                                                                               <  1:1635012/148‑1 (MQ=255)
  tcAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGcc                                                                                                                                              <  2:1000556/148‑1 (MQ=255)
  tcAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGcc                                                                                                                                              <  1:1923917/148‑1 (MQ=255)
    agcagTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCcg                                                                                                                                            >  1:1582794/1‑148 (MQ=255)
      cagTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCcgcg                                                                                                                                          >  2:716137/1‑148 (MQ=255)
      cagTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCcgcg                                                                                                                                          >  1:1724202/1‑148 (MQ=255)
         ttCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTgg                                                                                                                                       >  2:1221824/1‑148 (MQ=255)
         ttCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTgg                                                                                                                                       >  2:3308693/1‑148 (MQ=255)
         ttCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTgg                                                                                                                                       <  1:3172546/148‑1 (MQ=255)
         ttCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTgg                                                                                                                                       >  1:2736358/1‑148 (MQ=255)
            agtcagCTGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATa                                                                                                                                    <  1:1159253/148‑1 (MQ=255)
             gtcagCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCTCGCGTGGATaa                                                                                                                                   <  2:1407969/148‑1 (MQ=255)
              tcagCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCcgcg                                                                                                                                          >  1:1221921/1‑140 (MQ=255)
              tcagCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAAt                                                                                                                                  <  2:1198841/148‑1 (MQ=255)
               cagcagTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATc                                                                                                                                 >  2:642885/1‑148 (MQ=255)
                  cagTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTgg                                                                                                                              <  1:3333380/148‑1 (MQ=255)
                  cagTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTgg                                                                                                                              >  1:3490982/1‑148 (MQ=255)
                     tgtCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATg                                                                                                                           >  2:2381029/1‑148 (MQ=255)
                      gtCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGt                                                                                                                          <  2:1961061/148‑1 (MQ=255)
                      gtCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCAGTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGt                                                                                                                          <  2:2281710/148‑1 (MQ=255)
                       tCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGtt                                                                                                                         <  1:13705/148‑1 (MQ=255)
                           gTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGc                                                                                                                     <  2:467709/148‑1 (MQ=255)
                              aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTc                                                                                                                  >  1:1443735/1‑148 (MQ=255)
                              aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTc                                                                                                                  >  2:3080174/1‑148 (MQ=255)
                               tGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCt                                                                                                                 <  1:125388/148‑1 (MQ=255)
                               tGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCt                                                                                                                 <  1:2343031/148‑1 (MQ=255)
                                gACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCtt                                                                                                                >  1:2261768/1‑148 (MQ=255)
                                gACCTACGCTTTCGATGTGAACTCGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCtt                                                                                                                >  2:697545/1‑148 (MQ=255)
                                  ccTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCttgt                                                                                                              <  2:2161291/148‑1 (MQ=255)
                                     aCGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGa                                                                                                           >  2:1128060/1‑148 (MQ=255)
                                      cGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGaa                                                                                                          <  1:3079648/148‑1 (MQ=255)
                                         tttCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAg                                                                                                       <  1:2121700/148‑1 (MQ=255)
                                          ttCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAgg                                                                                                      >  2:2994835/1‑148 (MQ=255)
                                              aTGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCt                                                                                                  >  1:1245036/1‑148 (MQ=255)
                                               tgtgAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTg                                                                                                 <  2:2429817/148‑1 (MQ=255)
                                                 tgAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGAt                                                                                               <  1:1082285/148‑1 (MQ=255)
                                                    aCTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAgg                                                                                            <  1:2441717/148‑1 (MQ=255)
                                                     cTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGt                                                                                           <  1:124881/148‑1 (MQ=255)
                                                        gCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTc                                                                                        <  1:1515035/148‑1 (MQ=255)
                                                        gCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTc                                                                                        <  2:3451094/148‑1 (MQ=255)
                                                         cTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCa                                                                                       <  1:3536051/148‑1 (MQ=255)
                                                           cATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACg                                                                                     >  2:2875300/1‑148 (MQ=255)
                                                              ccGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTAc                                                                                  <  2:1930107/148‑1 (MQ=255)
                                                               cGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTAcc                                                                                 <  1:2908234/148‑1 (MQ=255)
                                                                 gagaGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAg                                                                               >  1:1206529/1‑148 (MQ=255)
                                                                 gagaGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAg                                                                               >  1:3064213/1‑148 (MQ=255)
                                                                       gtgGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCAt                                                                         >  1:700740/1‑148 (MQ=255)
                                                                        tgGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATc                                                                        >  2:1570098/1‑148 (MQ=255)
                                                                         ggACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCa                                                                       <  1:3293830/148‑1 (MQ=255)
                                                                                tGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGgaga                                                                <  1:396982/148‑1 (MQ=255)
                                                                                       tGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCg                                                         <  1:2584920/148‑1 (MQ=255)
                                                                                       tGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCg                                                         >  1:2091815/1‑148 (MQ=255)
                                                                                         gCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTc                                                       >  1:3337538/1‑148 (MQ=255)
                                                                                         gCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTc                                                       >  2:3222435/1‑148 (MQ=255)
                                                                                            cGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCgg                                                    >  1:2706233/1‑148 (MQ=255)
                                                                                             gATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGc                                                   <  2:2028471/148‑1 (MQ=255)
                                                                                                aGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTgc                                                >  1:3471371/1‑148 (MQ=255)
                                                                                                aGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTgc                                                >  1:1164973/1‑148 (MQ=255)
                                                                                                 gAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTgcc                                               >  1:1983120/1‑148 (MQ=255)
                                                                                                   aaCCATGTCCATTTCCGGATACGTTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGTATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTgccag                                             <  1:2418698/148‑1 (MQ=255)
                                                                                                    aCCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTgccagc                                            >  1:575173/1‑148 (MQ=255)
                                                                                                     ccATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTgccagcc                                           >  1:2953840/1‑148 (MQ=255)
                                                                                                       aTGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGccagccag                                         >  1:432631/1‑148 (MQ=255)
                                                                                                        tGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAgg                                        <  1:3095397/148‑1 (MQ=255)
                                                                                                         gTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACCGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGa                                       <  1:2985643/148‑1 (MQ=255)
                                                                                                          tCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGat                                      <  1:2383750/148‑1 (MQ=255)
                                                                                                           ccATTTCCGTATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGata                                     <  2:2449970/148‑1 (MQ=255)
                                                                                                           ccATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGata                                     <  1:268215/148‑1 (MQ=255)
                                                                                                           ccATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGata                                     <  1:3185881/148‑1 (MQ=255)
                                                                                                             aTTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGACTTAGCGTCCGGCTGCCAGCCAGGATATg                                   >  1:3218552/1‑148 (MQ=255)
                                                                                                             aTTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATg                                   >  1:824911/1‑148 (MQ=255)
                                                                                                              tttCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATgc                                  >  1:2490974/1‑148 (MQ=255)
                                                                                                               ttCCGTATACTGTTTTTTGAACAGACAACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATTGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATgcg                                 <  2:1592410/148‑1 (MQ=255)
                                                                                                               ttCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTCCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATgcg                                 >  1:319971/1‑148 (MQ=255)
                                                                                                                   ggATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCtt                             >  1:3273856/1‑148 (MQ=255)
                                                                                                                   ggATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCtt                             <  2:1724202/148‑1 (MQ=255)
                                                                                                                     aTACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTc                           >  2:263168/1‑148 (MQ=255)
                                                                                                                     aTACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTc                           >  1:1735702/1‑148 (MQ=255)
                                                                                                                      tACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCg                          >  2:2394403/1‑148 (MQ=255)
                                                                                                                      tACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCg                          >  2:2258461/1‑148 (MQ=255)
                                                                                                                          gTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTc                      >  2:1765828/1‑148 (MQ=255)
                                                                                                                          gTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTGTCGCTTc                      >  1:2639911/1‑148 (MQ=255)
                                                                                                                            tttttGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCgg                    >  1:1743895/1‑148 (MQ=255)
                                                                                                                            tttttGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCgg                    <  2:1703869/148‑1 (MQ=255)
                                                                                                                               ttGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTTATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGttt                 <  1:1721230/148‑1 (MQ=255)
                                                                                                                                  aaCAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTcg              >  1:863892/1‑148 (MQ=255)
                                                                                                                                   aCAGACCACATTCCAGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTcgc             >  2:1319545/1‑148 (MQ=255)
                                                                                                                                      gACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCgcc          >  1:1915406/1‑148 (MQ=255)
                                                                                                                                         cacaTTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCgccagc       <  2:3042966/148‑1 (MQ=255)
                                                                                                                                           cagtCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCgccagcca     <  2:3114896/145‑1 (MQ=255)
                                                                                                                                           caTTCCAGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCgccagcca     <  1:2790072/148‑1 (MQ=255)
                                                                                                                                             ttCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGt   <  2:1245036/148‑1 (MQ=255)
                                                                                                                                              tCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGtt  <  2:904240/148‑1 (MQ=255)
                                                                                                                                              tCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGtt  >  2:2424905/1‑148 (MQ=255)
                                                                                                                                              tCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGtt  <  2:2091815/148‑1 (MQ=255)
                                                                                                                                              tCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGtt  >  1:1793732/1‑148 (MQ=255)
                                                                                                                                              tCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGtt  <  1:1094698/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTT  >  CP009273/250764‑251053

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: ‑