Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,067,638 C→T 6.3% T873T (ACC→ACT flu → CP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,067,6380CT 6.3% 779.3 / 27.3 239T873T (ACC→ACTfluCP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter
Reads supporting (aligned to +/- strand):  ref base C (114/110);  new base T (9/6);  total (123/116)
Fisher's exact test for biased strand distribution p-value = 5.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 242 total aligned reads displayed.

AGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGCACAACATGTGCGTGCCGGTTTCC  >  CP009273/2067489‑2067786
                                                                                                                                                     |                                                                                                                                                    
 gCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCg                                                                                                                                                     <  1:2321889/150‑1 (MQ=255)
     gggCGGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCt                                                                                                                                                 >  2:930356/1‑150 (MQ=255)
     gggCGGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCt                                                                                                                                                 <  2:2517853/150‑1 (MQ=255)
              ccTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAg                                                                                                                                        >  2:3035532/1‑150 (MQ=255)
                  aaTCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATc                                                                                                                                    >  2:2800047/1‑150 (MQ=255)
                   aTCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCa                                                                                                                                   <  1:4349426/150‑1 (MQ=255)
                   aTCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCa                                                                                                                                   >  2:4709367/1‑150 (MQ=255)
                    tCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCAc                                                                                                                                  <  1:3558700/150‑1 (MQ=255)
                      tGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTg                                                                                                                                >  2:2107654/1‑150 (MQ=255)
                      tGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTg                                                                                                                                >  2:591902/1‑150 (MQ=255)
                        gTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGAc                                                                                                                              >  1:2002586/1‑150 (MQ=255)
                         tACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACa                                                                                                                             >  2:4389729/1‑150 (MQ=255)
                                gTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACTGGTCTGCCCTTCAGTATCACTGACAACCTGAt                                                                                                                      >  1:4754577/1‑150 (MQ=255)
                                gTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACTGGTCTGCCCTTCAGTATCACTGACAACCTGAt                                                                                                                      >  1:2169512/1‑150 (MQ=255)
                                 tcctccGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATg                                                                                                                     <  1:1880894/150‑1 (MQ=255)
                                   ctccGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCt                                                                                                                   <  1:1836208/150‑1 (MQ=255)
                                      cGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGa                                                                                                                <  2:1004522/150‑1 (MQ=255)
                                         ccTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGcc                                                                                                             >  2:2884603/1‑150 (MQ=255)
                                               ggCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACt                                                                                                       <  1:4530416/150‑1 (MQ=255)
                                                gCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTg                                                                                                      >  1:4934105/1‑150 (MQ=255)
                                                gCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTg                                                                                                      >  1:3302271/1‑150 (MQ=255)
                                                  tGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCa                                                                                                    >  1:3791636/1‑150 (MQ=255)
                                                     cATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGta                                                                                                 >  1:773423/1‑150 (MQ=255)
                                                     cATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGta                                                                                                 <  1:4276686/150‑1 (MQ=255)
                                                       ttGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGtata                                                                                               <  1:580593/150‑1 (MQ=255)
                                                          tgGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCt                                                                                            <  2:1076879/150‑1 (MQ=255)
                                                               caGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAg                                                                                       <  1:5172380/150‑1 (MQ=255)
                                                                aGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAgg                                                                                      >  1:414938/1‑150 (MQ=255)
                                                                 gggAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAggg                                                                                     >  2:3945893/1‑150 (MQ=255)
                                                                   gAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGAc                                                                                   >  2:4604870/1‑150 (MQ=255)
                                                                    aaCCCGCCACAGCATGAAAGCGTCATCGGACAATACCGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACt                                                                                  <  2:3832176/150‑1 (MQ=255)
                                                                    aaCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACt                                                                                  >  2:4487360/1‑150 (MQ=255)
                                                                     aCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACtt                                                                                 >  2:129605/1‑150 (MQ=255)
                                                                      cccGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACttt                                                                                >  2:4226317/1‑150 (MQ=255)
                                                                          ccACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACTGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCt                                                                            <  2:1237843/150‑1 (MQ=255)
                                                                          ccACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCt                                                                            >  2:4103459/1‑150 (MQ=255)
                                                                           cacaGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTg                                                                           <  2:2181585/150‑1 (MQ=255)
                                                                            acaGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACTGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTgg                                                                          >  1:487018/1‑150 (MQ=255)
                                                                            acaGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTgg                                                                          <  1:1628017/150‑1 (MQ=255)
                                                                                 aTGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGAc                                                                     >  1:4008543/1‑150 (MQ=255)
                                                                                   gAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACTGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACgg                                                                   >  2:492282/1‑150 (MQ=255)
                                                                                   gAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACgg                                                                   >  1:260113/1‑150 (MQ=255)
                                                                                    aaaGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGt                                                                  >  1:1681371/1‑150 (MQ=255)
                                                                                    aaaGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGt                                                                  <  1:1236430/150‑1 (MQ=255)
                                                                                     aaGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCGTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTa                                                                 <  1:375312/150‑1 (MQ=255)
                                                                                         gTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGa                                                             <  1:327086/150‑1 (MQ=255)
                                                                                         gTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGa                                                             >  2:3649703/1‑150 (MQ=255)
                                                                                          tCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGAc                                                            <  2:2810802/150‑1 (MQ=255)
                                                                                            aTCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACaa                                                          >  2:2010070/1‑150 (MQ=255)
                                                                                               ggACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGc                                                       >  1:1073287/1‑150 (MQ=255)
                                                                                                 aCAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCg                                                     >  1:1093116/1‑150 (MQ=255)
                                                                                                   aaTAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGt                                                   >  2:4768197/1‑150 (MQ=255)
                                                                                                    aTAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGtt                                                  <  2:5064740/150‑1 (MQ=255)
                                                                                                     tAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTa                                                 >  1:3668524/1‑150 (MQ=255)
                                                                                                     tAACGACTTCCGCGCCCGGGGCTGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTAt                                                <  2:969917/150‑1 (MQ=255)
                                                                                                      aaCGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTAt                                                <  2:1507211/150‑1 (MQ=255)
                                                                                                       aCGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTAtg                                               <  2:4816928/150‑1 (MQ=255)
                                                                                                       aCGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTAtg                                               <  2:2926827/150‑1 (MQ=255)
                                                                                                       aCGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTAtg                                               >  2:2347241/1‑150 (MQ=255)
                                                                                                          aCTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGa                                            <  1:2808213/150‑1 (MQ=255)
                                                                                                          aCTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGa                                            >  2:3243629/1‑150 (MQ=255)
                                                                                                           cTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGaa                                           >  2:3090795/1‑150 (MQ=255)
                                                                                                           cTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGaa                                           >  2:708021/1‑150 (MQ=255)
                                                                                                           cTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGaa                                           >  2:3381070/1‑150 (MQ=255)
                                                                                                              ccGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGtt                                        >  2:4993106/1‑150 (MQ=255)
                                                                                                              ccGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGGATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGtt                                        <  1:1533943/150‑1 (MQ=255)
                                                                                                               cgcgCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTc                                       <  2:3519007/150‑1 (MQ=255)
                                                                                                               cgcgCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTc                                       >  1:4375677/1‑150 (MQ=255)
                                                                                                                gcgcCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCg                                      >  1:4976946/1‑150 (MQ=255)
                                                                                                                gcgcCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCg                                      >  1:5120895/1‑150 (MQ=255)
                                                                                                                gcgcCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCg                                      >  2:2451766/1‑150 (MQ=255)
                                                                                                                 cgcCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCgg                                     >  2:706626/1‑150 (MQ=255)
                                                                                                                 cgcCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCgg                                     >  1:2236616/1‑150 (MQ=255)
                                                                                                                 cgcCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCgg                                     <  1:4900245/150‑1 (MQ=255)
                                                                                                                 cgcCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCgg                                     >  2:2742239/1‑150 (MQ=255)
                                                                                                                  gcCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCggg                                    <  2:828007/150‑1 (MQ=255)
                                                                                                                      ggggctggggctGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATg                                >  1:1278206/1‑150 (MQ=255)
                                                                                                                       gggctggggctgGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATgg                               >  1:686608/1‑150 (MQ=255)
                                                                                                                        ggctggggctggCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGc                              >  2:3298542/1‑150 (MQ=255)
                                                                                                                         gctggggctggCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCa                             <  1:1918300/150‑1 (MQ=255)
                                                                                                                          ctggggctggCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAg                            <  1:2969767/150‑1 (MQ=255)
                                                                                                                           tggggctggCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGt                           <  1:692582/150‑1 (MQ=255)
                                                                                                                           tggggctggCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGt                           <  2:4297682/150‑1 (MQ=255)
                                                                                                                              ggctggctGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGca                        <  1:644387/150‑1 (MQ=255)
                                                                                                                               gctggctgGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGcac                       <  1:1432415/150‑1 (MQ=255)
                                                                                                                                  ggctggGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGCacaac                    <  2:1691783/150‑1 (MQ=255)
                                                                                                                                   gctggGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGCacaaca                   >  2:3572999/1‑150 (MQ=255)
                                                                                                                                    ctggGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGCACAACAt                  <  1:2108130/150‑1 (MQ=255)
                                                                                                                                      gggCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGCACAACAtgt                <  2:5134463/150‑1 (MQ=255)
                                                                                                                                       ggCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGCACAACAtgtg               >  2:4039740/1‑150 (MQ=255)
                                                                                                                                           cACTGGAAACTGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGCACAACATgtgcgtg           >  1:4454628/1‑150 (MQ=255)
                                                                                                                                             cTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTCAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGCACAACATGTGCGTGcc         >  1:5081012/1‑150 (MQ=255)
                                                                                                                                             cTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGCACAACATGTGCGTGcc         <  1:2507907/150‑1 (MQ=255)
                                                                                                                                              tGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGCACAACATGTGCGTGCCg        <  2:2981703/150‑1 (MQ=255)
                                                                                                                                               ggAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGCACAACATGTGCGTGCCgg       >  1:3509332/1‑150 (MQ=255)
                                                                                                                                               ggAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGCACAACATGTGCGTGCCgg       >  2:147381/1‑150 (MQ=255)
                                                                                                                                                gAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGCACAACATGTGCGTGCCGGt      >  2:1551164/1‑150 (MQ=255)
                                                                                                                                                  aaCTGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGCACAACATGTGCGTGCCGGttt    <  2:885723/150‑1 (MQ=255)
                                                                                                                                                  aaCCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGCACAACATGTGCGTGCCGGttt    >  1:2067387/1‑150 (MQ=255)
                                                                                                                                                  aaCCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGCACAACATGTGCGTGCCGGttt    >  1:295863/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGCACAACATGTGCGTGCCGGTTTCC  >  CP009273/2067489‑2067786

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: