Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,996,506 G→T 8.8% T194T (ACC→ACA fliC ← flagellar filament structural protein (flagellin)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,996,5060GT 8.8% 813.1 / 36.6 216T194T (ACC→ACAfliCflagellar filament structural protein (flagellin)
Reads supporting (aligned to +/- strand):  ref base G (87/110);  new base T (13/6);  total (100/116)
Fisher's exact test for biased strand distribution p-value = 5.41e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 227 total aligned reads displayed.

TCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTT  >  CP009273/1996358‑1996655
                                                                                                                                                    |                                                                                                                                                     
tCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCACTTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTgtgg                                                                                                                                                      <  2:1621422/150‑1 (MQ=255)
tCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTgtgg                                                                                                                                                      <  1:952469/150‑1 (MQ=255)
     tCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTa                                                                                                                                                 <  1:4130172/150‑1 (MQ=255)
     tCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTa                                                                                                                                                 <  1:3158485/150‑1 (MQ=255)
      cGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCCCTTGAAGCTGGGTTAGGTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAg                                                                                                                                                <  2:3109656/150‑1 (MQ=255)
       gTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGc                                                                                                                                               <  1:4184288/150‑1 (MQ=255)
       gTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGc                                                                                                                                               <  1:3440690/150‑1 (MQ=255)
        ttATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTTGTGGTAGCa                                                                                                                                              <  2:4466531/150‑1 (MQ=255)
         tATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCAc                                                                                                                                             <  1:3342438/150‑1 (MQ=255)
          aTCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCAcc                                                                                                                                            >  1:1932303/1‑150 (MQ=255)
             accaccGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCaaa                                                                                                                                         <  1:803287/150‑1 (MQ=255)
             accaccGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCaaa                                                                                                                                         <  1:2026060/150‑1 (MQ=255)
             accaccGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCaaa                                                                                                                                         <  1:3656427/150‑1 (MQ=255)
                   ggTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGTTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGt                                                                                                                                   <  2:4514174/150‑1 (MQ=255)
                    gTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTTGTGGTAGCACCAAAAGCAGtt                                                                                                                                  >  2:3026333/1‑150 (MQ=255)
                       aTTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACt                                                                                                                               <  1:2692767/150‑1 (MQ=255)
                        ttttCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTg                                                                                                                              <  1:3245761/150‑1 (MQ=255)
                           tCGCATAGTAATCATTNCCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTNAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAg                                                                                                                           >  2:4342824/1‑150 (MQ=255)
                            cGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTTGTGGTAGCACCAAAAGCAGTTACTGGAGc                                                                                                                          >  2:3653500/1‑150 (MQ=255)
                             gCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCa                                                                                                                         >  1:2576493/1‑150 (MQ=255)
                              cATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCAc                                                                                                                        >  2:4132361/1‑150 (MQ=255)
                              cATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCAc                                                                                                                        >  2:1805068/1‑150 (MQ=255)
                               aTAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACt                                                                                                                       <  2:4020062/150‑1 (MQ=255)
                                 aGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAg                                                                                                                     >  1:2443960/1‑150 (MQ=255)
                                  gTAATCATTACCATTATCAGTATACACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGt                                                                                                                    <  2:1517413/150‑1 (MQ=255)
                                  gTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGt                                                                                                                    <  1:859481/150‑1 (MQ=255)
                                   tAATCATTACCATTCTCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTg                                                                                                                   <  2:2202142/150‑1 (MQ=255)
                                   tAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTg                                                                                                                   <  2:1094635/150‑1 (MQ=255)
                                    aaTCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTgg                                                                                                                  >  1:1730918/1‑150 (MQ=255)
                                       cATTACCATTATCAGTATAAACACCCTCACTTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTaa                                                                                                               <  2:3186726/150‑1 (MQ=255)
                                         ttACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTTGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACt                                                                                                             >  2:2916295/1‑150 (MQ=255)
                                         ttACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACt                                                                                                             <  1:1138232/150‑1 (MQ=255)
                                           aCCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGt                                                                                                           <  1:852209/150‑1 (MQ=255)
                                            ccATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTa                                                                                                          >  1:1539034/1‑150 (MQ=255)
                                             cATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTAt                                                                                                         >  2:3806594/1‑150 (MQ=255)
                                              aTTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATc                                                                                                        <  2:4221593/150‑1 (MQ=255)
                                              aTTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATc                                                                                                        <  2:2711344/150‑1 (MQ=255)
                                                     gTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTAttt                                                                                                 <  1:3681602/150‑1 (MQ=255)
                                                      tataAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTAtttt                                                                                                >  1:4256585/1‑150 (MQ=255)
                                                        taAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTa                                                                                              >  2:2350399/1‑150 (MQ=255)
                                                        taAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTa                                                                                              >  2:2349703/1‑150 (MQ=255)
                                                         aaaCACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTaa                                                                                             >  2:1325013/1‑150 (MQ=255)
                                                         aaaCACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTaa                                                                                             >  2:2633635/1‑150 (MQ=255)
                                                         aaaCACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTaa                                                                                             >  2:2826331/1‑150 (MQ=255)
                                                          aaCACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAAc                                                                                            >  2:2962353/1‑150 (MQ=255)
                                                             acCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCt                                                                                         <  1:1015176/150‑1 (MQ=255)
                                                              cccTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTa                                                                                        <  1:1993920/150‑1 (MQ=255)
                                                                cTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTTGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTaaa                                                                                      >  1:2015457/1‑150 (MQ=255)
                                                                   aaTTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAAcc                                                                                   >  1:3890479/1‑150 (MQ=255)
                                                                        aaGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCaa                                                                              >  1:1753552/1‑150 (MQ=255)
                                                                           cTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGc                                                                           >  1:1967319/1‑150 (MQ=255)
                                                                            tGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGcc                                                                          <  1:1157310/150‑1 (MQ=255)
                                                                              tgTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCaa                                                                        <  2:2698589/149‑1 (MQ=255)
                                                                              tgTTAGTTCCGCCAGTATCAGTGGCTGCTTCAGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGTTAGCACCAAAATCATTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCaa                                                                        <  2:1107077/149‑1 (MQ=255)
                                                                                   gTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGtt                                                                   <  1:1639717/150‑1 (MQ=255)
                                                                                   gTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGtt                                                                   <  1:3837529/150‑1 (MQ=255)
                                                                                   gTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGtt                                                                   <  1:4408795/150‑1 (MQ=255)
                                                                                    ttCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGttt                                                                  <  1:2997375/150‑1 (MQ=255)
                                                                                    ttCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGttt                                                                  <  2:2285162/150‑1 (MQ=255)
                                                                                        gCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGc                                                              >  2:123771/1‑150 (MQ=255)
                                                                                        gCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGc                                                              <  2:3567670/150‑1 (MQ=255)
                                                                                        gCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGc                                                              <  2:3523141/150‑1 (MQ=255)
                                                                                            gTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCa                                                          >  2:27643/1‑150 (MQ=255)
                                                                                            gTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCa                                                          <  1:3410508/150‑1 (MQ=255)
                                                                                              aTCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAAt                                                        <  2:419086/150‑1 (MQ=255)
                                                                                              aTCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAAt                                                        >  2:1952948/1‑150 (MQ=255)
                                                                                                cAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCt                                                      <  1:4519979/150‑1 (MQ=255)
                                                                                                   tGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCt                                                   <  1:327600/150‑1 (MQ=255)
                                                                                                    ggCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTTGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCtt                                                  >  1:520430/1‑150 (MQ=255)
                                                                                                    ggCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCtt                                                  >  2:844215/1‑150 (MQ=255)
                                                                                                    ggCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCtt                                                  >  1:2847880/1‑150 (MQ=255)
                                                                                                      ctgcgTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCa                                                <  2:3666890/150‑1 (MQ=255)
                                                                                                       tgctTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAg                                               <  1:2066560/150‑1 (MQ=255)
                                                                                                           tCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCg                                           <  1:24076/150‑1 (MQ=255)
                                                                                                           tCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCg                                           <  1:24086/150‑1 (MQ=255)
                                                                                                              gTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCgata                                        >  2:2793503/1‑150 (MQ=255)
                                                                                                                aGAAATGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCgatagt                                      >  1:4087229/1‑150 (MQ=255)
                                                                                                                aGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTTGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCgatagt                                      >  1:697056/1‑150 (MQ=255)
                                                                                                                aGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCgatagt                                      <  2:313961/150‑1 (MQ=255)
                                                                                                                   aaGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCgatagtgat                                   <  2:3976197/150‑1 (MQ=255)
                                                                                                                    aGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCgatagtgata                                  >  2:2049303/1‑150 (MQ=255)
                                                                                                                    aGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCgatagtgata                                  >  2:1713140/1‑150 (MQ=255)
                                                                                                                     gggTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCgatagtgatag                                 <  2:2096211/150‑1 (MQ=255)
                                                                                                                     gggTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCgatagtgatag                                 <  2:1386756/150‑1 (MQ=255)
                                                                                                                          aTTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGGTTCAGATCGATAGTGATAGTCTgg                            <  1:2623582/150‑1 (MQ=255)
                                                                                                                          aTTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTgg                            <  1:3567150/150‑1 (MQ=255)
                                                                                                                            tCCAGTAAGTTTAATATTGTTTGTTGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGtt                          <  1:2785/150‑1 (MQ=255)
                                                                                                                             ccAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTa                         >  1:2817877/1‑150 (MQ=255)
                                                                                                                                  aaGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCAtt                    <  2:545392/150‑1 (MQ=255)
                                                                                                                                  aaGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCAtt                    <  2:1931905/150‑1 (MQ=255)
                                                                                                                                   aGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCAttt                   <  1:3133090/150‑1 (MQ=255)
                                                                                                                                   aGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCAttt                   <  2:813203/150‑1 (MQ=255)
                                                                                                                                   aGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCAttt                   <  2:2493224/150‑1 (MQ=255)
                                                                                                                                     tttAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTgc                 <  1:2108253/150‑1 (MQ=255)
                                                                                                                                       tAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTgcgc               >  1:1154337/1‑150 (MQ=255)
                                                                                                                                        aaTATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGcc              >  1:3085342/1‑150 (MQ=255)
                                                                                                                                         atatTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTACCGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCa             <  2:4366428/150‑1 (MQ=255)
                                                                                                                                           atTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAAc           <  1:720811/150‑1 (MQ=255)
                                                                                                                                               tttgtTGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTgg       <  2:3898343/150‑1 (MQ=255)
                                                                                                                                                   tggtggTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGAttt   >  1:878032/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TCCCATCGTTATCACCACCGGTGATTTTCGCATAGTAATCATTACCATTATCAGTATAAACACCCTCAATTGAAGCTGGGTTAGTTCCGCCAGTATCAGTGGCTGCTTCCGTAGAAAGGGTAATTCCAGTAAGTTTAATATTGTTTGTGGTGGTAGCACCAAAAGCAGTTACTGGAGCACTAGTGGTAACTGTATCGTTATTTTTAACGCTAAAACCATCAAGGCCAAGAGTTTTAGCATCAATCTGCTTCAGATCGATAGTGATAGTCTGGTTATCATTTGCGCCAACCTGGATTTT  >  CP009273/1996358‑1996655

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: