Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 557,592 G→T 100% D814Y (GAC→TAC)  sfmD → putative outer membrane export usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273557,5920GT97.3% 805.1 / 10.3 259D814Y (GAC→TAC) sfmDputative outer membrane export usher protein
Reads supporting (aligned to +/- strand):  ref base G (3/4);  new base T (129/123);  total (132/127)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 261 total aligned reads displayed.

CCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTG  >  CP009273/557445‑557737
                                                                                                                                                   |                                                                                                                                                 
ccGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTAc                                                                                                                                                 <  2:2153315/150‑1 (MQ=255)
ccGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTAc                                                                                                                                                 <  1:5090085/150‑1 (MQ=255)
 cGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACg                                                                                                                                                >  1:3303595/1‑150 (MQ=255)
 cGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACg                                                                                                                                                >  1:504964/1‑150 (MQ=255)
  gACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACgg                                                                                                                                               >  1:2585522/1‑150 (MQ=255)
      cAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTCTGGTGGGTGATTACGGGCaa                                                                                                                                           >  2:1401529/1‑150 (MQ=255)
      cAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCaa                                                                                                                                           >  1:1368935/1‑150 (MQ=255)
       aaGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAg                                                                                                                                          >  1:3952108/1‑150 (MQ=255)
       aaGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAg                                                                                                                                          <  2:1559862/150‑1 (MQ=255)
        aGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCTGGGAGGAAAACCGGTGCCGTTTGGATCACTGGTCCGGGACAACAGTACCGGAATAACCAGTCTGGTGGGTGATTACGGGCACGt                                                                                                                                         >  2:3743252/1‑150 (MQ=255)
         ggCGCGTTGGTTCGTGCCAATTTTGTTACACGCATAGGCGTGCGGGCGCTCATTACCGTTAACCAGGGAGGAAAACCGGTGCAGTTTGGATAACTGGTACGGGAAAACAGTACCGGAATAACAAGTATGGTGTGTGATTAAGGGCAAGtt                                                                                                                                        >  1:3144973/1‑150 (MQ=255)
         ggCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGtt                                                                                                                                        >  1:2163308/1‑150 (MQ=255)
           cgcgTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGgtta                                                                                                                                      >  2:4799810/1‑150 (MQ=255)
            gcgTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGtttat                                                                                                                                     >  2:1123056/1‑150 (MQ=255)
              gTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGtttattt                                                                                                                                   >  1:4144388/1‑150 (MQ=255)
               ttGGTTCGAGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGtttattta                                                                                                                                  <  1:4328906/150‑1 (MQ=255)
                 ggTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAg                                                                                                                                >  1:2417209/1‑150 (MQ=255)
                  gTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCCAGTTTATTTAAgt                                                                                                                               >  2:4616414/1‑150 (MQ=255)
                  gTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAgt                                                                                                                               <  2:838090/150‑1 (MQ=255)
                   ttCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAAACAGTATTGTGGGTGATTACGGGCAAGTTAATTTAAgtg                                                                                                                              >  2:2625688/1‑150 (MQ=255)
                     cGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAgtggt                                                                                                                            >  1:2701696/1‑150 (MQ=255)
                      gTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTGCGGGCAAGTTTATTTAAgtggtg                                                                                                                           >  1:1633686/1‑150 (MQ=255)
                      gTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAgtggtg                                                                                                                           >  1:1562823/1‑150 (MQ=255)
                            aTTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCAt                                                                                                                     >  1:453394/1‑150 (MQ=255)
                             ttttGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCAtt                                                                                                                    <  2:101566/150‑1 (MQ=255)
                              tttGAGACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTg                                                                                                                   <  1:962321/150‑1 (MQ=255)
                                tGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTc                                                                                                                 <  2:4523822/150‑1 (MQ=255)
                                     cccGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTCCCGGAATAACCAGTATGGTGGGTGAGTACGGGCAAGGTGATTTAAGTGGTGCGCCATGGTCGGGTg                                                                                                            >  2:2634109/1‑150 (MQ=255)
                                     cccGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTg                                                                                                            <  1:3196441/150‑1 (MQ=255)
                                       cGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGaa                                                                                                          <  2:2409017/150‑1 (MQ=255)
                                        gCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAAt                                                                                                         <  1:4292137/150‑1 (MQ=255)
                                          aTAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTa                                                                                                       <  1:2265154/150‑1 (MQ=255)
                                          aTAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTa                                                                                                       <  2:2067542/150‑1 (MQ=255)
                                                 tGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTc                                                                                                >  1:4824154/1‑150 (MQ=255)
                                                 tGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTc                                                                                                <  2:3203993/150‑1 (MQ=255)
                                                     ggCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTg                                                                                            >  2:1844944/1‑150 (MQ=255)
                                                       cgcTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTggg                                                                                          <  2:361715/150‑1 (MQ=255)
                                                          tCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGgag                                                                                       <  1:4680560/150‑1 (MQ=255)
                                                           cATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGgaga                                                                                      >  1:1049588/1‑150 (MQ=255)
                                                            aTTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGAc                                                                                     <  2:5015279/150‑1 (MQ=255)
                                                            aTTACCGTTACCCAGGGCGGAAAACCGGTGCAGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGAc                                                                                     >  2:475032/1‑150 (MQ=255)
                                                             ttACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACg                                                                                    <  1:4765690/150‑1 (MQ=255)
                                                             ttACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACg                                                                                    >  1:1472259/1‑150 (MQ=255)
                                                              tACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACgg                                                                                   >  1:2600477/1‑150 (MQ=255)
                                                                ccGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGgcg                                                                                 <  2:3699578/150‑1 (MQ=255)
                                                                ccGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGgcg                                                                                 <  2:2118379/150‑1 (MQ=255)
                                                                ccGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGgcg                                                                                 >  2:687869/1‑150 (MQ=255)
                                                                 cGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGgcgc                                                                                >  1:5157492/1‑150 (MQ=255)
                                                                 cGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGgcgc                                                                                >  1:1286029/1‑150 (MQ=255)
                                                                 cGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGgcgc                                                                                >  1:2080913/1‑150 (MQ=255)
                                                                 cGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGgcgc                                                                                >  1:4379991/1‑150 (MQ=255)
                                                                   ttACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGa                                                                              <  1:3910172/150‑1 (MQ=255)
                                                                       ccAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTc                                                                          <  2:1690009/150‑1 (MQ=255)
                                                                         cgggCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCAc                                                                        >  2:3813777/2‑150 (MQ=255)
                                                                         aGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCAc                                                                        <  1:3368590/150‑1 (MQ=255)
                                                                         aGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCAc                                                                        <  1:4521889/150‑1 (MQ=255)
                                                                         aGGGCGGAAAACCGGTGCCGTGTGGATCACTGGTAAGGGAAAACAGTACCGGAATAACCAGGATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCGGGTGAATTACAGGTTCAGTGGGGAGACTGAGCGAACTCAc                                                                        >  1:1785775/1‑150 (MQ=255)
                                                                           ggCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTGCCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGc                                                                      <  2:4072784/150‑1 (MQ=255)
                                                                            gCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCt                                                                     >  1:3787234/1‑150 (MQ=255)
                                                                                  aaCCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTg                                                               <  2:1618149/150‑1 (MQ=255)
                                                                                    ccGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTgcg                                                             >  2:2042534/1‑150 (MQ=255)
                                                                                      ggTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCa                                                           >  2:610556/1‑150 (MQ=255)
                                                                                      ggTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCa                                                           >  2:436957/1‑150 (MQ=255)
                                                                                       gTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCAc                                                          >  1:2512135/1‑150 (MQ=255)
                                                                                          ccGTTTGGATCACTGGTACGGGAAAACCGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTAt                                                       <  1:3306247/150‑1 (MQ=255)
                                                                                           cGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATg                                                      <  2:1995367/150‑1 (MQ=255)
                                                                                           cGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATg                                                      <  2:3514932/150‑1 (MQ=255)
                                                                                           cGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATg                                                      <  1:1103664/150‑1 (MQ=255)
                                                                                                ggATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTg                                                 >  1:4939414/1‑150 (MQ=255)
                                                                                                    cACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGACGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGa                                             >  2:1113987/1‑150 (MQ=255)
                                                                                                      cTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAg                                           <  2:3303603/150‑1 (MQ=255)
                                                                                                      cTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAg                                           <  1:2748312/150‑1 (MQ=255)
                                                                                                       tGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGc                                          >  2:2770072/1‑150 (MQ=255)
                                                                                                         gTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCaa                                        <  1:1860249/150‑1 (MQ=255)
                                                                                                         gTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCaa                                        <  1:1860191/150‑1 (MQ=255)
                                                                                                           aCGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAg                                      <  2:3100090/150‑1 (MQ=255)
                                                                                                                aaaaCAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTAc                                 <  2:2364801/150‑1 (MQ=255)
                                                                                                                 aaaCAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACa                                <  1:668928/150‑1 (MQ=255)
                                                                                                                 aaaCAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACa                                <  1:3045924/150‑1 (MQ=255)
                                                                                                                   aCAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGc                              >  2:1777986/1‑150 (MQ=255)
                                                                                                                    cAGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCa                             >  1:1088980/1‑150 (MQ=255)
                                                                                                                     aGTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCaa                            >  1:4986012/1‑150 (MQ=255)
                                                                                                                      gTACCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAg                           >  2:3417810/1‑150 (MQ=255)
                                                                                                                        aCCGGAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGcc                         >  2:4242097/1‑150 (MQ=255)
                                                                                                                           ggAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTc                      >  2:4917629/1‑150 (MQ=255)
                                                                                                                            gAATAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCa                     <  2:3240209/150‑1 (MQ=255)
                                                                                                                               tAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTg                  >  2:2329791/1‑150 (MQ=255)
                                                                                                                               tAACCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTg                  >  2:3357286/1‑150 (MQ=255)
                                                                                                                                aaCCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGt                 >  1:4852475/1‑150 (MQ=255)
                                                                                                                                 aCCAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGtt                <  1:1123056/150‑1 (MQ=255)
                                                                                                                                   cAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTAt              >  2:3423733/1‑150 (MQ=255)
                                                                                                                                   cAGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTAt              <  2:1333359/150‑1 (MQ=255)
                                                                                                                                    aGTATGGTGGGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTAtt             <  1:3582813/150‑1 (MQ=255)
                                                                                                                                        tggtggGTGATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCgg         >  1:407230/1‑150 (MQ=255)
                                                                                                                                         ggtgggtgATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGc        <  2:535214/150‑1 (MQ=255)
                                                                                                                                          gtgggtgATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCa       <  2:289195/150‑1 (MQ=255)
                                                                                                                                           tgggtgATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAg      >  2:3325554/1‑150 (MQ=255)
                                                                                                                                            cggtgATTACGGGCGGGGTTTTTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGt     >  2:1310573/2‑150 (MQ=255)
                                                                                                                                             ggtgATTACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGtt    <  2:3952108/150‑1 (MQ=255)
                                                                                                                                               tgattaCGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTg  <  1:2844574/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CCGACGCAAGGCGCGTTGGTTCGTGCCAATTTTGATACCCGCATAGGCGTGCGGGCGCTCATTACCGTTACCCAGGGCGGAAAACCGGTGCCGTTTGGATCACTGGTACGGGAAAACAGTACCGGAATAACCAGTATGGTGGGTGATGACGGGCAAGTTTATTTAAGTGGTGCGCCATTGTCTGGTGAATTACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTG  >  CP009273/557445‑557737

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