Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,024,827 G→T 100% D103Y (GAT→TAT)  pepQ → proline dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,024,8270GT97.3% 681.5 / 8.7 224D103Y (GAT→TAT) pepQproline dipeptidase
Reads supporting (aligned to +/- strand):  ref base G (3/3);  new base T (102/116);  total (105/119)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 225 total aligned reads displayed.

GCAATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAACCCGAAAGGGG  >  CP009273/4024682‑4024971
                                                                                                                                                 |                                                                                                                                                
gCAATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAatat                                                                                                                                                <  2:2906267/148‑1 (MQ=255)
 ctaTTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATAtg                                                                                                                                               <  2:3278840/146‑1 (MQ=255)
  aaTTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATAtgt                                                                                                                                              <  1:3005407/148‑1 (MQ=255)
   aTTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATAtgtg                                                                                                                                             >  1:3561656/1‑148 (MQ=255)
   aTTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATAtgtg                                                                                                                                             >  1:1844814/1‑148 (MQ=255)
     tCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGa                                                                                                                                           >  2:1484428/1‑148 (MQ=255)
     tCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGa                                                                                                                                           <  2:940508/148‑1 (MQ=255)
     tCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGa                                                                                                                                           >  2:3546020/1‑148 (MQ=255)
     tCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGa                                                                                                                                           >  2:413380/1‑148 (MQ=255)
     tCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGa                                                                                                                                           >  1:531925/1‑148 (MQ=255)
     tCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGa                                                                                                                                           >  2:2210324/1‑148 (MQ=255)
         aGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTg                                                                                                                                       >  1:106192/1‑148 (MQ=255)
          gCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGa                                                                                                                                      >  1:802300/1‑148 (MQ=255)
           cGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGAt                                                                                                                                     <  1:917290/148‑1 (MQ=255)
           cGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGAt                                                                                                                                     >  2:7970/1‑148 (MQ=255)
               ggTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATcgcg                                                                                                                                 <  2:3349353/148‑1 (MQ=255)
               ggTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATcgcg                                                                                                                                 <  1:2068848/148‑1 (MQ=255)
                 tGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCt                                                                                                                               <  2:179618/148‑1 (MQ=255)
                         aCTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAg                                                                                                                       <  1:3222909/148‑1 (MQ=255)
                           tCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGcc                                                                                                                     >  2:1338789/1‑148 (MQ=255)
                               gTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATg                                                                                                                 >  2:1518110/1‑148 (MQ=255)
                                  ccAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCa                                                                                                              <  1:2775586/148‑1 (MQ=255)
                                   cAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCAt                                                                                                             >  2:1706365/1‑148 (MQ=255)
                                   cAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCAt                                                                                                             <  1:797900/148‑1 (MQ=255)
                                    aaaCTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCAtt                                                                                                            <  1:1640228/148‑1 (MQ=255)
                                       ctgctgGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGt                                                                                                         <  2:166306/148‑1 (MQ=255)
                                       ctgctgGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGt                                                                                                         >  1:2574885/1‑148 (MQ=255)
                                        tgctggttgctggtGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTa                                                                                                        >  2:2465981/1‑148 (MQ=255)
                                        tgctggttgctggtGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTa                                                                                                        >  2:2223829/1‑148 (MQ=255)
                                             gttgctggtGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAgcctg                                                                                                   <  2:115764/148‑1 (MQ=255)
                                             gttgctggtGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAgcctg                                                                                                   <  2:115749/148‑1 (MQ=255)
                                             gttgctggtGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAgcctg                                                                                                   <  1:3357512/148‑1 (MQ=255)
                                              ttgctggtGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAgcctgc                                                                                                  <  2:1384728/148‑1 (MQ=255)
                                              ttgctggtGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAgcctgc                                                                                                  <  2:2588800/148‑1 (MQ=255)
                                                   ggtggATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAgcctgctgcct                                                                                             >  1:714375/1‑148 (MQ=255)
                                                    gtggATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAgcctgctgcctg                                                                                            <  1:2841755/148‑1 (MQ=255)
                                                     tggATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAgcctgctgcctgc                                                                                           <  2:1215786/148‑1 (MQ=255)
                                                     tggATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAgcctgctgcctgc                                                                                           <  2:2427966/148‑1 (MQ=255)
                                                       gATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGcctgctgcctgctg                                                                                         >  1:315976/1‑148 (MQ=255)
                                                           gCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTgcgcg                                                                                     >  2:558156/1‑148 (MQ=255)
                                                              tGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCgg                                                                                  >  2:748767/1‑148 (MQ=255)
                                                                aaCAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCa                                                                                >  1:3268116/1‑148 (MQ=255)
                                                                aaCAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCa                                                                                >  1:2522444/1‑148 (MQ=255)
                                                                   aaGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAata                                                                             <  1:2649498/148‑1 (MQ=255)
                                                                   aaGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAata                                                                             >  2:887276/1‑148 (MQ=255)
                                                                    aGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAatat                                                                            >  2:1030203/1‑148 (MQ=255)
                                                                    aGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAatat                                                                            >  1:661538/1‑148 (MQ=255)
                                                                        gAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGt                                                                        <  1:1328761/148‑1 (MQ=255)
                                                                        gAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGt                                                                        <  1:1003099/148‑1 (MQ=255)
                                                                        gAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATAAGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGt                                                                        <  2:2263426/148‑1 (MQ=255)
                                                                         aaaCTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGtt                                                                       <  2:3507836/148‑1 (MQ=255)
                                                                           aCTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTtat                                                                     <  1:146804/148‑1 (MQ=255)
                                                                           aCTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTtat                                                                     <  2:714375/148‑1 (MQ=255)
                                                                            cTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTtata                                                                    >  1:1586348/1‑148 (MQ=255)
                                                                                 gTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGt                                                               >  1:1495294/1‑148 (MQ=255)
                                                                                 gTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGt                                                               <  1:2310998/148‑1 (MQ=255)
                                                                                  ttCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTc                                                              <  2:3050091/148‑1 (MQ=255)
                                                                                   tCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTcc                                                             >  2:664952/1‑148 (MQ=255)
                                                                                     tATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCgg                                                           >  1:1743481/1‑148 (MQ=255)
                                                                                          gCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCg                                                      >  1:2345053/1‑148 (MQ=255)
                                                                                            cGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGa                                                    <  1:1818421/148‑1 (MQ=255)
                                                                                             ggTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGaa                                                   <  1:2845846/148‑1 (MQ=255)
                                                                                              gTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGTGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAAc                                                  >  2:3237454/1‑148 (MQ=255)
                                                                                               ttGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGTAACg                                                 >  1:111955/1‑148 (MQ=255)
                                                                                               ttGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACg                                                 >  1:2260949/1‑148 (MQ=255)
                                                                                               ttGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACg                                                 >  2:2598216/1‑148 (MQ=255)
                                                                                                tGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGt                                                >  1:3487281/1‑148 (MQ=255)
                                                                                                 gATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCATGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTg                                               <  1:2651198/148‑1 (MQ=255)
                                                                                                   ttACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTgcg                                             >  1:2995015/1‑148 (MQ=255)
                                                                                                   ttACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTgcg                                             >  1:1213653/1‑148 (MQ=255)
                                                                                                     aCTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCt                                           >  2:1431153/1‑148 (MQ=255)
                                                                                                       tGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGc                                         >  1:2382925/1‑148 (MQ=255)
                                                                                                          cacaACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAAc                                      <  2:2380127/148‑1 (MQ=255)
                                                                                                           acaACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACt                                     >  1:172534/1‑148 (MQ=255)
                                                                                                            caACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTg                                    >  2:593337/1‑148 (MQ=255)
                                                                                                             aaCGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTgg                                   >  1:1781585/1‑148 (MQ=255)
                                                                                                             aaCGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTgg                                   >  1:1781472/1‑148 (MQ=255)
                                                                                                             aaCGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTgg                                   >  1:3265716/1‑148 (MQ=255)
                                                                                                              aCGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTggg                                  >  2:911138/1‑148 (MQ=255)
                                                                                                               cGTCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGt                                 >  2:1933760/1‑148 (MQ=255)
                                                                                                                 tCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTAt                               <  2:546290/148‑1 (MQ=255)
                                                                                                                 tCGAACCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTAt                               <  2:1984500/148‑1 (MQ=255)
                                                                                                                    aaCCGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGa                            >  1:3426768/1‑148 (MQ=255)
                                                                                                                       cGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGTCAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGc                         <  1:413411/148‑1 (MQ=255)
                                                                                                                       cGCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGc                         >  1:1035748/1‑148 (MQ=255)
                                                                                                                        gCTGCCGACCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGcc                        <  2:3505328/148‑1 (MQ=255)
                                                                                                                             cGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCaa                   <  1:887276/148‑1 (MQ=255)
                                                                                                                               aCCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACa                 >  2:320591/1‑148 (MQ=255)
                                                                                                                               aCCTCTTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACa                 <  1:1518110/148‑1 (MQ=255)
                                                                                                                                  tctTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCa              >  2:2683336/1‑148 (MQ=255)
                                                                                                                                  tctTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCa              >  1:2372080/1‑148 (MQ=255)
                                                                                                                                   ctgtCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCaa             <  2:2528721/145‑1 (MQ=255)
                                                                                                                                   ctTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCaa             <  1:2223829/148‑1 (MQ=255)
                                                                                                                                   ctTTCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCaa             <  1:1338789/148‑1 (MQ=255)
                                                                                                                                    tttCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAAc            >  1:3466527/1‑148 (MQ=255)
                                                                                                                                     ttCTGGACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAAcc           <  1:2467702/148‑1 (MQ=255)
                                                                                                                                          gACTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAACCCGaaa      <  2:71312/148‑1 (MQ=255)
                                                                                                                                            cTGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAACCCGAAAgg    <  2:3054855/148‑1 (MQ=255)
                                                                                                                                             tGAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAACCCGAAAggg   <  2:2720344/148‑1 (MQ=255)
                                                                                                                                              gAATATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAACCCGAAAgggg  <  2:833728/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GCAATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAACGTCGAACCGCTGCCGACCTCTTTCTGGACTGAAGATGTGGAAGTGATCGCGCTGCCGAAAGCCGATGGCATTGGTAGCCTGCTGCCTGCTGCGCGCGGCAATATCGGTTATATCGGTCCGGTGCCGGAACGTGCGCTGCAACTGGGTATTGAGGCCAGCAACATCAACCCGAAAGGGG  >  CP009273/4024682‑4024971

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: