Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,730,144 A→T 5.1% A170A (GCA→GCT)  sodB → superoxide dismutase, Fe

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,730,1440AT 5.1% 617.1 / 18.0 196A170A (GCA→GCT) sodBsuperoxide dismutase, Fe
Reads supporting (aligned to +/- strand):  ref base A (98/88);  new base T (6/4);  total (104/92)
Fisher's exact test for biased strand distribution p-value = 7.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 196 total aligned reads displayed.

GCTGGACCTGGCTGGTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGG  >  CP009273/1729999‑1730290
                                                                                                                                                 |                                                                                                                                                  
gctggACCTGGCTGGTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACg                                                                                                                                                  >  2:1272665/1‑148 (MQ=255)
  tggACCTGGCTGGTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTc                                                                                                                                                <  1:2808023/148‑1 (MQ=255)
  tggACCTGGCTGGTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTc                                                                                                                                                <  1:3416552/148‑1 (MQ=255)
   ggACCTGGCTGGTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTcc                                                                                                                                               <  1:2646744/148‑1 (MQ=255)
         ggctggTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTa                                                                                                                                         >  2:540485/1‑148 (MQ=255)
          gctggTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTAt                                                                                                                                        >  1:2499424/1‑148 (MQ=255)
          gctggTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTAt                                                                                                                                        >  1:781509/1‑148 (MQ=255)
                gAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAg                                                                                                                                  >  2:2225974/1‑148 (MQ=255)
                 aaaaaCAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGc                                                                                                                                 >  2:1166034/1‑148 (MQ=255)
                  aaaaCAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCa                                                                                                                                >  2:304353/1‑148 (MQ=255)
                   aaaCAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCACTGCTCGTCCTGGCTATCTGGAGCAc                                                                                                                               >  2:2605575/1‑148 (MQ=255)
                   aaaCAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCAc                                                                                                                               >  2:2128263/1‑148 (MQ=255)
                    aaCAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACt                                                                                                                              >  2:1416945/1‑148 (MQ=255)
                     aCAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACtt                                                                                                                             >  2:2282811/1‑148 (MQ=255)
                     aCAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACtt                                                                                                                             >  2:2786387/1‑148 (MQ=255)
                     aCAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACtt                                                                                                                             <  1:1062016/148‑1 (MQ=255)
                        gCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTCTCTGGAGCACTTCTg                                                                                                                          >  2:1668697/1‑148 (MQ=255)
                        gCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTg                                                                                                                          >  2:1668726/1‑148 (MQ=255)
                         cGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTgg                                                                                                                         >  2:1891290/1‑148 (MQ=255)
                         cGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTgg                                                                                                                         >  1:2994650/1‑148 (MQ=255)
                         cGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTgg                                                                                                                         >  1:9487/1‑148 (MQ=255)
                         cGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTgg                                                                                                                         >  1:2969532/1‑148 (MQ=255)
                         cGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTgg                                                                                                                         >  1:10114/1‑148 (MQ=255)
                           aTGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGATATCTGGAGCACTTCTGGgc                                                                                                                       >  2:1068418/1‑148 (MQ=255)
                            tGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGgcg                                                                                                                      <  2:2575654/148‑1 (MQ=255)
                            tGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGgcg                                                                                                                      <  1:1999119/148‑1 (MQ=255)
                              gCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCt                                                                                                                    <  1:3050226/148‑1 (MQ=255)
                              gCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCt                                                                                                                    >  2:2815819/1‑148 (MQ=255)
                                 aaCTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGt                                                                                                                 <  1:1079299/148‑1 (MQ=255)
                                 aaCTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGt                                                                                                                 <  2:3341115/148‑1 (MQ=255)
                                 aaCTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGt                                                                                                                 >  1:2366029/1‑148 (MQ=255)
                                 aaCTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGt                                                                                                                 >  1:485141/1‑148 (MQ=255)
                                  aCTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTg                                                                                                                >  2:2939045/1‑148 (MQ=255)
                                    tGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGaa                                                                                                              <  1:942905/148‑1 (MQ=255)
                                     ggCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCTCGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAAc                                                                                                             >  1:219990/1‑148 (MQ=255)
                                         aTCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTggg                                                                                                         >  2:2038613/1‑148 (MQ=255)
                                         aTCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTggg                                                                                                         >  1:2028367/1‑148 (MQ=255)
                                          tCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGa                                                                                                        <  1:1272665/148‑1 (MQ=255)
                                          tCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGa                                                                                                        >  1:2271832/1‑148 (MQ=255)
                                            gTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAAt                                                                                                      <  2:1990907/148‑1 (MQ=255)
                                              ttCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTc                                                                                                    >  1:2299561/1‑148 (MQ=255)
                                                  aCCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAg                                                                                                >  2:598805/1‑148 (MQ=255)
                                                      ctAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGaa                                                                                            >  2:2468537/1‑148 (MQ=255)
                                                            cgGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAAtct                                                                                      <  2:2489294/148‑1 (MQ=255)
                                                            cgGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAAtct                                                                                      >  2:523911/1‑148 (MQ=255)
                                                               gTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGc                                                                                   <  2:2340693/148‑1 (MQ=255)
                                                                  cTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTgg                                                                                                                         <  1:911438/107‑1 (MQ=255)
                                                                  cTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGc                                                                                >  1:1884440/1‑148 (MQ=255)
                                                                  cTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGc                                                                                <  2:1524832/148‑1 (MQ=255)
                                                                  cTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGc                                                                                >  2:1024477/1‑148 (MQ=255)
                                                                     cGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCata                                                                             >  1:2730733/1‑148 (MQ=255)
                                                                        tGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCataata                                                                          <  2:418157/148‑1 (MQ=255)
                                                                         gACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCTCGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCAtaataa                                                                         >  1:1006982/1‑148 (MQ=255)
                                                                         gACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCAtaataa                                                                         <  2:1267416/148‑1 (MQ=255)
                                                                          accaccGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAAc                                                                        <  2:2717685/148‑1 (MQ=255)
                                                                          accaccGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAAc                                                                        >  2:3059639/1‑148 (MQ=255)
                                                                          accaccGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAAc                                                                        <  2:3468154/148‑1 (MQ=255)
                                                                           ccaccGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGACCTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACt                                                                       >  1:3548116/1‑148 (MQ=255)
                                                                           ccaccGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACt                                                                       >  1:1334145/1‑148 (MQ=255)
                                                                             accGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGa                                                                     <  1:874016/148‑1 (MQ=255)
                                                                             accGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGa                                                                     <  2:709562/148‑1 (MQ=255)
                                                                              ccGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCTCGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGAt                                                                    <  2:3394409/148‑1 (MQ=255)
                                                                                 aTGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGc                                                                 >  1:43374/1‑148 (MQ=255)
                                                                                    cGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCCATGCACGTCCTGGCTATCTGGAGCACGTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCaat                                                              >  1:543617/1‑147 (MQ=255)
                                                                                    cGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCaaa                                                              <  1:2825857/148‑1 (MQ=255)
                                                                                    cGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCaaa                                                              >  1:1487974/1‑148 (MQ=255)
                                                                                      aCTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATg                                                            <  1:1329394/148‑1 (MQ=255)
                                                                                        tCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATgca                                                          <  1:1222031/148‑1 (MQ=255)
                                                                                         ccGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATgcag                                                         <  2:3524081/148‑1 (MQ=255)
                                                                                         ccGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATgcag                                                         <  2:255126/148‑1 (MQ=255)
                                                                                                tgACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGcc                                                  <  2:566468/148‑1 (MQ=255)
                                                                                                  aCCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTg                                                >  1:712548/1‑148 (MQ=255)
                                                                                                   ccGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGa                                               <  2:1942000/148‑1 (MQ=255)
                                                                                                    cGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGaa                                              <  1:1977645/148‑1 (MQ=255)
                                                                                                         aTGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCTCGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCt                                         <  2:1308706/148‑1 (MQ=255)
                                                                                                         aTGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCt                                         >  2:1467909/1‑148 (MQ=255)
                                                                                                            tCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACg                                      >  2:1483953/1‑148 (MQ=255)
                                                                                                            tCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACg                                      <  2:1323703/148‑1 (MQ=255)
                                                                                                             cTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCACATGCAGCATTGCCTGAAGCGCTACGc                                     >  2:1774654/1‑148 (MQ=255)
                                                                                                                   acacGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCa                               >  1:157913/1‑148 (MQ=255)
                                                                                                                    cacGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAg                              >  1:539362/1‑148 (MQ=255)
                                                                                                                      cGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGc                            <  1:3277397/148‑1 (MQ=255)
                                                                                                                        cTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCt                          <  2:2188212/148‑1 (MQ=255)
                                                                                                                        cTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCt                          <  2:852136/148‑1 (MQ=255)
                                                                                                                         ttattaCATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTa                         >  1:2825773/1‑148 (MQ=255)
                                                                                                                         ttattaCATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTa                         <  1:1829013/148‑1 (MQ=255)
                                                                                                                         ttattaCATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTa                         <  1:1829333/148‑1 (MQ=255)
                                                                                                                         ttattaCATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTa                         <  1:3136579/148‑1 (MQ=255)
                                                                                                                             taCATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTAcgcg                     >  2:309181/1‑148 (MQ=255)
                                                                                                                             taCATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTAcgcg                     >  2:820005/1‑148 (MQ=255)
                                                                                                                                aTCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCCAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGat                  >  2:3178161/1‑148 (MQ=255)
                                                                                                                                aTCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGat                  >  1:381389/1‑148 (MQ=255)
                                                                                                                                 tCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGatc                 <  1:2068802/148‑1 (MQ=255)
                                                                                                                                 tCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGatc                 >  1:659074/1‑148 (MQ=255)
                                                                                                                                 tCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGatc                 >  2:948116/1‑148 (MQ=255)
                                                                                                                                     cTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCg             <  2:1884440/148‑1 (MQ=255)
                                                                                                                                          gCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTa        <  2:781509/148‑1 (MQ=255)
                                                                                                                                             aTGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGt     <  2:1717583/148‑1 (MQ=255)
                                                                                                                                             aTGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGt     <  2:1717329/148‑1 (MQ=255)
                                                                                                                                                cACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCgg  <  1:3005307/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GCTGGACCTGGCTGGTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGG  >  CP009273/1729999‑1730290

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