Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,576,325 T→A 100% N153Y (AAT→TAT)  ydeN ← putative Ser‑type periplasmic non‑aryl sulfatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,576,3250TA99.0% 650.9 / ‑2.4 197N153Y (AAT→TAT) ydeNputative Ser‑type periplasmic non‑aryl sulfatase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base A (134/61);  total (135/62)
Fisher's exact test for biased strand distribution p-value = 5.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 203 total aligned reads displayed.

GCGTTTGTTTATCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGC  >  CP009273/1576181‑1576469
                                                                                                                                                |                                                                                                                                                
gcgTTTGTTTATCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGa                                                                                                                                               <  2:2503787/148‑1 (MQ=255)
  gtttgtttATCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAAt                                                                                                                                             <  1:507217/148‑1 (MQ=255)
    ttgtttATCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAg                                                                                                                                           >  1:2551462/1‑148 (MQ=255)
           tCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCaa                                                                                                                                    >  2:3368916/1‑148 (MQ=255)
           tCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCaa                                                                                                                                    >  2:684047/1‑148 (MQ=255)
              tCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAgc                                                                                                                                 >  1:3650631/1‑148 (MQ=255)
               ccGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATTGGAATAGACACCAAAgcg                                                                                                                                <  1:2825699/148‑1 (MQ=255)
               ccGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAgcg                                                                                                                                <  1:1737056/148‑1 (MQ=255)
                   tACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAgcgggcg                                                                                                                            >  2:3026033/1‑148 (MQ=255)
                    aCCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAgcgggcgg                                                                                                                           <  2:3207443/148‑1 (MQ=255)
                     ccGGCACATTACTGATTTTTGACACGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGcgggcggg                                                                                                                          <  1:1609404/148‑1 (MQ=255)
                     ccGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGcgggcggg                                                                                                                          <  1:1609353/148‑1 (MQ=255)
                        gCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGc                                                                                                                       <  1:2960012/148‑1 (MQ=255)
                        gCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGc                                                                                                                       >  1:1197628/1‑148 (MQ=255)
                         cacaTTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCt                                                                                                                      >  1:1342246/1‑148 (MQ=255)
                            aTTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGa                                                                                                                   >  1:3183742/1‑148 (MQ=255)
                            aTTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGa                                                                                                                   >  2:3247965/1‑148 (MQ=255)
                              tACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATGATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGAcc                                                                                                                 >  2:49871/1‑148 (MQ=255)
                              tACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGAcc                                                                                                                 >  1:1461245/1‑148 (MQ=255)
                              tACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGAcc                                                                                                                 >  1:1461157/1‑148 (MQ=255)
                               aCTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCg                                                                                                                >  2:2141803/1‑148 (MQ=255)
                                cTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCgg                                                                                                               <  1:994515/148‑1 (MQ=255)
                                  gATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTc                                                                                                             >  1:2080592/1‑148 (MQ=255)
                                    tttttGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCat                                                                                                           >  2:631112/1‑148 (MQ=255)
                                       ttGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCattat                                                                                                        <  2:579275/148‑1 (MQ=255)
                                          aCAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTgc                                                                                                     >  2:1411013/1‑148 (MQ=255)
                                           cAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTgcg                                                                                                    >  2:392903/1‑148 (MQ=255)
                                            aaGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACAACAAAGCGGGCGGGAGCTCGACCGGTCATTATTgcgg                                                                                                   <  1:313459/148‑1 (MQ=255)
                                              gTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTgcggcg                                                                                                 <  1:3549099/148‑1 (MQ=255)
                                              gTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTgcggcg                                                                                                 <  1:3694360/148‑1 (MQ=255)
                                              gTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTgcggcg                                                                                                 >  2:2468428/1‑148 (MQ=255)
                                                 ccATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCgg                                                                                              <  2:2149693/148‑1 (MQ=255)
                                                    tttACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAg                                                                                           >  2:1286031/1‑148 (MQ=255)
                                                    tttACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAg                                                                                           >  2:2199875/1‑148 (MQ=255)
                                                    tttACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAg                                                                                           >  1:3632546/1‑148 (MQ=255)
                                                     ttACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAgg                                                                                          >  2:1340885/1‑148 (MQ=255)
                                                     ttACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAgg                                                                                          <  2:314767/148‑1 (MQ=255)
                                                     ttACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAgg                                                                                          >  2:1788603/1‑148 (MQ=255)
                                                     ttACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAgg                                                                                          <  1:2344429/148‑1 (MQ=255)
                                                     ttACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAgg                                                                                          >  1:2345752/1‑148 (MQ=255)
                                                         cTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGcc                                                                                      >  1:305946/1‑148 (MQ=255)
                                                          tACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCg                                                                                     >  2:3570558/1‑148 (MQ=255)
                                                          tACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCg                                                                                     >  1:2501667/1‑148 (MQ=255)
                                                            ctgctgCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGa                                                                                   >  1:198132/1‑148 (MQ=255)
                                                             tgctgcAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGaa                                                                                  <  1:3601631/148‑1 (MQ=255)
                                                             tgctgcAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGaa                                                                                  <  1:2449659/148‑1 (MQ=255)
                                                               ctgcAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAac                                                                                <  1:2466355/148‑1 (MQ=255)
                                                                tgcAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAaca                                                                               >  2:1126588/1‑148 (MQ=255)
                                                                tgcAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAaca                                                                               >  2:3383558/1‑148 (MQ=255)
                                                                  cAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACAcc                                                                             >  2:2973677/1‑148 (MQ=255)
                                                                    gtgtAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCgt                                                                           >  1:803662/1‑148 (MQ=255)
                                                                       taataaCCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGtgtg                                                                        >  2:915051/1‑148 (MQ=255)
                                                                         ataaCCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGcc                                                                      >  1:79917/1‑148 (MQ=255)
                                                                           aaCCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCcac                                                                    >  1:502615/1‑148 (MQ=255)
                                                                               aTGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAg                                                                >  1:533710/1‑148 (MQ=255)
                                                                               aTGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAg                                                                >  2:860650/1‑148 (MQ=255)
                                                                                 gATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGcc                                                              >  1:1546849/1‑148 (MQ=255)
                                                                                  aTTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCg                                                             <  1:1725762/148‑1 (MQ=255)
                                                                                    tCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGtt                                                           <  1:1365763/148‑1 (MQ=255)
                                                                                     cTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTa                                                          <  1:3303176/148‑1 (MQ=255)
                                                                                         aataatTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTaa                                                      <  1:1052586/148‑1 (MQ=255)
                                                                                          ataatTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGATCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTaaa                                                     <  2:2987138/148‑1 (MQ=255)
                                                                                          ataatTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTaaa                                                     >  2:1450605/1‑148 (MQ=255)
                                                                                           taatTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGTTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTaaaa                                                    <  1:3522100/148‑1 (MQ=255)
                                                                                           taatTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTaaaa                                                    >  2:1901343/1‑148 (MQ=255)
                                                                                            aatTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAAc                                                   >  1:2652608/1‑148 (MQ=255)
                                                                                               tCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTa                                                >  1:1112522/1‑148 (MQ=255)
                                                                                               tCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTa                                                >  2:3602856/1‑148 (MQ=255)
                                                                                               tCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTa                                                >  1:2811675/1‑148 (MQ=255)
                                                                                                  ggCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGc                                             <  1:3252084/148‑1 (MQ=255)
                                                                                                   gCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGcc                                            >  1:599640/1‑148 (MQ=255)
                                                                                                    cAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCt                                           >  2:3708499/1‑148 (MQ=255)
                                                                                                     aaGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCtt                                          <  1:2751860/148‑1 (MQ=255)
                                                                                                        aaaGTTTCTGTTAGAGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCAt                                       >  2:3533577/1‑148 (MQ=255)
                                                                                                           gTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCa                                    >  2:2078263/1‑148 (MQ=255)
                                                                                                            tttCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCAt                                   <  2:3220694/148‑1 (MQ=255)
                                                                                                            tttCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCAt                                   <  2:3694234/148‑1 (MQ=255)
                                                                                                            tttCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCAt                                   >  2:3259459/1‑148 (MQ=255)
                                                                                                             ttCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCAtt                                  >  2:3161268/1‑148 (MQ=255)
                                                                                                              tCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTa                                 <  2:2551462/148‑1 (MQ=255)
                                                                                                                 gTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATg                              >  2:1561230/1‑148 (MQ=255)
                                                                                                                   tAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGaa                            >  2:2873463/1‑148 (MQ=255)
                                                                                                                   tAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGaa                            >  2:962137/1‑148 (MQ=255)
                                                                                                                    aGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGaaa                           >  1:1355017/1‑148 (MQ=255)
                                                                                                                     gCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAg                          >  2:2264126/1‑148 (MQ=255)
                                                                                                                        gAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAg                       <  2:3527102/148‑1 (MQ=255)
                                                                                                                         aaTACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGc                      >  2:2646949/1‑148 (MQ=255)
                                                                                                                          aTACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCg                     >  2:1480315/1‑148 (MQ=255)
                                                                                                                           tACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGt                    <  1:992031/148‑1 (MQ=255)
                                                                                                                           tACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGt                    >  2:7641/1‑148 (MQ=255)
                                                                                                                            aCCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTc                   >  1:2902814/1‑148 (MQ=255)
                                                                                                                            aCCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTc                   >  1:2902636/1‑148 (MQ=255)
                                                                                                                            aCCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACAGAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTc                   >  2:239985/1‑148 (MQ=255)
                                                                                                                                tCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCg               <  2:1676578/148‑1 (MQ=255)
                                                                                                                                         tCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGAttttt      >  1:1779840/1‑148 (MQ=255)
                                                                                                                                         tCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGAttttt      >  1:1336898/1‑148 (MQ=255)
                                                                                                                                             tataGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGc  >  1:687639/1‑148 (MQ=255)
                                                                                                                                             tataGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGc  <  2:3629842/148‑1 (MQ=255)
                                                                                                                                             tataGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGc  <  2:1197628/148‑1 (MQ=255)
                                                                                                                                             tataGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTGGATTTTTGTGc  >  1:3519325/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
GCGTTTGTTTATCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGC  >  CP009273/1576181‑1576469

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