Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,471,930 T→A 88.2% P575P (CCA→CCT valS ← valyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,471,9300TA88.2% 542.7 / 51.5 185P575P (CCA→CCTvalSvalyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (13/9);  new base A (74/89);  total (87/98)
Fisher's exact test for biased strand distribution p-value = 2.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 185 total aligned reads displayed.

TAGACGCCAGCGCCGCCAGGGTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCTGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCACCTGCGGTTTGCCATTTTCATCTTTGA  >  CP009273/4471785‑4472075
                                                                                                                                                 |                                                                                                                                                 
tAGACGCCAGCGCCGCCAGGGTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAgg                                                                                                                                                 <  2:694099/148‑1 (MQ=255)
tAGACGCCAGCGCCGCCAGGGTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAgg                                                                                                                                                 <  2:2608053/148‑1 (MQ=255)
     gCCAGCGCCGCCAGGGTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCg                                                                                                                                            <  2:2693908/148‑1 (MQ=255)
       cAGCGCCGCCAGGGTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGaa                                                                                                                                          >  1:1380686/1‑148 (MQ=255)
       cAGCGCCGCCAGGGTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGaa                                                                                                                                          <  1:1182101/148‑1 (MQ=255)
        aGCGCCGCCAGGGTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGaaa                                                                                                                                         >  1:1026777/1‑148 (MQ=255)
        aGCGCCGCCAGGGTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGaaa                                                                                                                                         <  1:812847/148‑1 (MQ=255)
         gcgcCGCCAGGGTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAAt                                                                                                                                        >  1:412087/1‑148 (MQ=255)
          cgccgcCAGGGTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCTGGCAGCGAAATa                                                                                                                                       >  1:897520/1‑148 (MQ=255)
           gccgccAGGGTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATAc                                                                                                                                      <  1:972005/148‑1 (MQ=255)
           gccgccAGGGTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATAc                                                                                                                                      <  1:482917/148‑1 (MQ=255)
            ccgccAGGGTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATAcc                                                                                                                                     <  1:1712093/148‑1 (MQ=255)
                cAGGGTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTCCGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCa                                                                                                                                 >  1:2846363/1‑148 (MQ=255)
                       aaGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCatat                                                                                                                          >  1:2119827/1‑148 (MQ=255)
                          cgcagcgcGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCa                                                                                                                       >  1:1981430/1‑148 (MQ=255)
                            cagcgcGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGt                                                                                                                     <  2:1241777/148‑1 (MQ=255)
                            cagcgcGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGt                                                                                                                     <  2:1241804/148‑1 (MQ=255)
                               cgcgTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTggg                                                                                                                  <  1:109817/148‑1 (MQ=255)
                               cgcgTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTggg                                                                                                                  <  1:25712/148‑1 (MQ=255)
                                     aGTACCTTGCGGCTAAATAACGTTCGGGAACAGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGAATCGAAATACAGTCAACCATATCCAGTGGGTCGATa                                                                                                            >  2:2916581/1‑148 (MQ=255)
                                      gTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATaa                                                                                                           <  1:120644/148‑1 (MQ=255)
                                       tACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAAc                                                                                                          >  2:16324/1‑148 (MQ=255)
                                         ccGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGt                                                                                                        >  2:1330702/1‑148 (MQ=255)
                                         ccGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGt                                                                                                        >  2:1033246/1‑148 (MQ=255)
                                            tGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTAc                                                                                                     >  1:1762327/1‑148 (MQ=255)
                                            tGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAAAGTTAc                                                                                                     >  1:1227265/1‑148 (MQ=255)
                                             gCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTAcc                                                                                                    <  1:692990/148‑1 (MQ=255)
                                              cGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTAccc                                                                                                   <  2:3241209/148‑1 (MQ=255)
                                              cGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTAccc                                                                                                   <  1:3083814/148‑1 (MQ=255)
                                                   cAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCTGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGa                                                                                              >  2:910137/1‑148 (MQ=255)
                                                   cAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCTGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGa                                                                                              >  2:910109/1‑148 (MQ=255)
                                                     aTACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCTGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGAtt                                                                                            <  1:1834828/148‑1 (MQ=255)
                                                     aTACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGAtt                                                                                            >  2:1665537/1‑148 (MQ=255)
                                                      tACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGAttt                                                                                           >  1:1404059/1‑148 (MQ=255)
                                                      tACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGAttt                                                                                           >  1:3313255/1‑148 (MQ=255)
                                                      tACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGAttt                                                                                           >  1:326018/1‑148 (MQ=255)
                                                      tACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGAttt                                                                                           >  2:573159/1‑148 (MQ=255)
                                                      tACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGAttt                                                                                           >  2:22220/1‑148 (MQ=255)
                                                       aCCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTg                                                                                          <  1:3215375/148‑1 (MQ=255)
                                                        ccGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTgg                                                                                         <  2:149213/148‑1 (MQ=255)
                                                         cGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGa                                                                                        <  2:1521242/148‑1 (MQ=255)
                                                             cGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACAtc                                                                                    <  2:897994/148‑1 (MQ=255)
                                                                gAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACAtcttc                                                                                 >  2:170298/1‑148 (MQ=255)
                                                                gAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACAtcttc                                                                                 >  2:1106131/1‑148 (MQ=255)
                                                                 aaCTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACAtcttct                                                                                <  2:973009/148‑1 (MQ=255)
                                                                  aCTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTg                                                                               >  1:3158398/1‑148 (MQ=255)
                                                                  aCTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTg                                                                               >  2:2930458/1‑148 (MQ=255)
                                                                     gCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGcc                                                                            <  1:2666100/148‑1 (MQ=255)
                                                                      cTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCt                                                                           <  1:3133653/148‑1 (MQ=255)
                                                                       ttCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCtt                                                                          >  1:3008042/1‑148 (MQ=255)
                                                                       ttCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCtt                                                                          >  2:2126906/1‑148 (MQ=255)
                                                                           cGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGtc                                                                      >  2:826846/1‑148 (MQ=255)
                                                                           cGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGtc                                                                      >  2:2485057/1‑148 (MQ=255)
                                                                              tGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGtcatc                                                                   >  1:429036/1‑148 (MQ=255)
                                                                              tGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGtcatc                                                                   >  1:952000/1‑148 (MQ=255)
                                                                               gcgcTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGtcatca                                                                  >  1:2620776/1‑148 (MQ=255)
                                                                                cgcTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCAc                                                                 >  1:1705427/1‑148 (MQ=255)
                                                                                  cTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGa                                                               <  1:977310/148‑1 (MQ=255)
                                                                                    tACGGATTTTGTCCTCCAGCTGCGGCTGAATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTAACCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAAt                                                             <  1:1774633/148‑1 (MQ=255)
                                                                                    tACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAAt                                                             <  1:1256598/148‑1 (MQ=255)
                                                                                       ggATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAg                                                          <  2:2356046/148‑1 (MQ=255)
                                                                                        gATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAgg                                                         <  1:258894/148‑1 (MQ=255)
                                                                                         aTTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGc                                                        <  2:274605/148‑1 (MQ=255)
                                                                                              gTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTc                                                   <  1:1033246/148‑1 (MQ=255)
                                                                                                   ccAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTa                                              <  2:2822271/148‑1 (MQ=255)
                                                                                                      gCTGCGGCTGCATGATAGTGGCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCGAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGTTCATGTAAAc                                           <  2:2875621/148‑1 (MQ=255)
                                                                                                       cTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACg                                          <  1:746753/148‑1 (MQ=255)
                                                                                                            gCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGtgtg                                     >  2:1474857/1‑148 (MQ=255)
                                                                                                             cTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTgg                                    >  1:888642/1‑148 (MQ=255)
                                                                                                               gCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGaa                                  <  1:1595680/148‑1 (MQ=255)
                                                                                                               gCATCATATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGaa                                  <  1:2162388/148‑1 (MQ=255)
                                                                                                                catcatATTGCCGGTACGTTTTTCCAGCAGTTCTGGCAGCGAAATACCGTCAACCATTTCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAAc                                 >  2:1233800/1‑148 (MQ=255)
                                                                                                                catcatATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAAc                                 >  1:1344734/1‑148 (MQ=255)
                                                                                                                catcatATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAAc                                 >  2:2840497/1‑148 (MQ=255)
                                                                                                                 atcatATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACg                                >  2:820924/1‑148 (MQ=255)
                                                                                                                   catATTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGc                              <  1:2235284/148‑1 (MQ=255)
                                                                                                                      atTGCCGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCAcc                           <  1:1028713/148‑1 (MQ=255)
                                                                                                                          ccGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCACCTGCg                       >  1:2358103/1‑148 (MQ=255)
                                                                                                                           cGGTACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCACCTGCgg                      >  1:733836/1‑148 (MQ=255)
                                                                                                                              tACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCACCTGCGGttt                   >  2:1587005/1‑148 (MQ=255)
                                                                                                                              tACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCACCTGCGGttt                   >  2:817198/1‑148 (MQ=255)
                                                                                                                              tACGTTTTTCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCACCTGCGGttt                   >  2:1587319/1‑148 (MQ=255)
                                                                                                                                cGTTTTTCCAGCAGTTCTGGCAGCGAAATACCGTCAACCATATAAAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGACGGTCATGTAAACGGTTTGGAACGGCACCTGCGGTTTGc                 >  2:351544/1‑148 (MQ=255)
                                                                                                                                  tttttCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCACCTGCGGTTTGCCa               <  2:3218061/148‑1 (MQ=255)
                                                                                                                                   ttttCCAGCAGTTCTGGCAGCGAAATACCGTCAACCATTTCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCACCTGCGGTTTGCCAt              >  1:2146198/1‑148 (MQ=255)
                                                                                                                                   ttttCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCACCTGCGGTTTGCCAt              >  1:2902785/1‑148 (MQ=255)
                                                                                                                                   ttttCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCACCTGCGGTTTGCCAt              >  1:18174/1‑148 (MQ=255)
                                                                                                                                    tttCCAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCACCTGCGGTTTGCCAtt             >  1:1227388/1‑148 (MQ=255)
                                                                                                                                       ccAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCACCTGCGGTTTGCCATTTTc          <  2:412074/148‑1 (MQ=255)
                                                                                                                                       ccAGCAGTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCACCTGCGGTTTGCCATTTTc          <  2:3174771/148‑1 (MQ=255)
                                                                                                                                          gcagTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCACCTGCGGTTTGCCATTTTCATc       <  1:62187/148‑1 (MQ=255)
                                                                                                                                           cagTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCTGGCCGGTCATGTAAACGGTGTGGAACGGCACCTGCGGTTTGCCATTTTCATCt      <  1:337901/148‑1 (MQ=255)
                                                                                                                                           cagTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCACCTGCGGTTTGCCATTTTCATCt      <  1:266152/148‑1 (MQ=255)
                                                                                                                                           cagTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCACCTGCGGTTTGCCATTTTCATCt      <  1:2961119/148‑1 (MQ=255)
                                                                                                                                             gTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCACCTGCGGTTTGCCATTTTCATCttt    >  2:431519/1‑148 (MQ=255)
                                                                                                                                             gTTCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCACCTGCGGTTTGCCATTTTCATCttt    <  2:2501645/148‑1 (MQ=255)
                                                                                                                                              ttCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCACCTGCGGTTTGCCCTTTTCATCTTtg   >  2:1764402/1‑148 (MQ=255)
                                                                                                                                              ttCAGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCACCTGCGGTTTGCCATTTTCATCTTtg   <  1:2625977/148‑1 (MQ=255)
                                                                                                                                               tCTGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCACCTGCGGTTTGCCATTTTCATCTTtga  <  1:910109/148‑1 (MQ=255)
                                                                                                                                               tCTGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCACCTGCGGTTTGCCATTTTCATCTTtga  <  1:910137/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TAGACGCCAGCGCCGCCAGGGTGAAGCGCAGCGCGTCAGTACCGTGCGGCTCAATACCGTTCGGGAACTGCTTCTCGGTGCGCTTACGGATTTTGTCCGCCAGCTGCGGCTGCATCATATTGCCGGTACGTTTTTCCAGCAGTTCTGGCAGCGAAATACCGTCAACCATATCCAGTGGGTCGATAACGTTACCCTTGGATTTGGACATCTTCTGGCCTTCGTCATCACGAATCAGGCCGGTCATGTAAACGGTGTGGAACGGCACCTGCGGTTTGCCATTTTCATCTTTGA  >  CP009273/4471785‑4472075

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: