Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,529,866 C→T 6.7% intergenic (‑548/‑908) nanC ← / → fimB N‑acetylnuraminic acid outer membrane channel protein/tyrosine recombinase/inversion of on/off regulator of fimA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,529,8660CT 6.7% 1100.0 / 44.2 317intergenic (‑548/‑908)nanC/fimBN‑acetylnuraminic acid outer membrane channel protein/tyrosine recombinase/inversion of on/off regulator of fimA
Reads supporting (aligned to +/- strand):  ref base C (154/141);  new base T (14/7);  total (168/149)
Fisher's exact test for biased strand distribution p-value = 2.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 324 total aligned reads displayed.

TGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAA  >  CP009273/4529717‑4530015
                                                                                                                                                     |                                                                                                                                                     
 gTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAaca                                                                                                                                                      <  2:1581508/150‑1 (MQ=255)
      aTTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGc                                                                                                                                                 >  1:2505067/1‑150 (MQ=255)
      aTTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAAAAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCAATGTGAAGCCAGACGAACACTTga                                                                                                                                                 >  2:2803487/1‑149 (MQ=255)
           cTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAATACTTGCggtgg                                                                                                                                            >  1:1413740/1‑150 (MQ=255)
           cTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCggtgg                                                                                                                                            <  1:3898108/150‑1 (MQ=255)
            tATCTTGGTGCACCTTTATACCTGTTCTACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAAGTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCggtggt                                                                                                                                           <  2:4006255/150‑1 (MQ=255)
             aTCTTGGTGCACCTTTATCCCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTc                                                                                                                                          <  2:1685243/150‑1 (MQ=255)
             aTCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTc                                                                                                                                          <  2:579689/150‑1 (MQ=255)
              tCTTTGTGCACCTTTATACCTGTTATACCAGATCAAAAATAACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCt                                                                                                                                         <  2:4529749/150‑1 (MQ=255)
                ttGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTc                                                                                                                                       <  2:4025669/150‑1 (MQ=255)
                    tGCACCTTTCTACCTGTTATAACAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCaaaa                                                                                                                                   <  2:1275269/150‑1 (MQ=255)
                     gCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAATACTTGCGGTGGTCTTCaaaaa                                                                                                                                  >  1:1372567/1‑150 (MQ=255)
                      cACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAAc                                                                                                                                 <  1:2438315/150‑1 (MQ=255)
                          tttATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTaaa                                                                                                                             <  1:4103360/150‑1 (MQ=255)
                            tataCCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGa                                                                                                                           >  1:2203328/1‑150 (MQ=255)
                             ataCCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCATATTTGCATTTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGAt                                                                                                                          <  1:2842388/150‑1 (MQ=255)
                             ataCCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGAt                                                                                                                          <  1:4123839/150‑1 (MQ=255)
                              taCCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGAta                                                                                                                         >  2:4839202/1‑149 (MQ=255)
                              taCCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATc                                                                                                                         >  1:2069987/1‑150 (MQ=255)
                                     tataCCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGttt                                                                                                                  >  1:8532/1‑150 (MQ=255)
                                       taCCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAATACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTaa                                                                                                                >  2:3225192/1‑150 (MQ=255)
                                       taCCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTaa                                                                                                                <  2:2477641/150‑1 (MQ=255)
                                          cAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTa                                                                                                             >  1:4970399/1‑150 (MQ=255)
                                           aGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTAt                                                                                                            <  1:2923642/150‑1 (MQ=255)
                                             ataaaaaaTCACGCAATCCATACAACAAAACCTGTTTTGCAATTCGTGTCACAAAATTTGTCGTTCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTAttt                                                                                                          <  1:3504398/147‑1 (MQ=255)
                                              tCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTg                                                                                                         <  2:5039484/150‑1 (MQ=255)
                                               cAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAATACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGt                                                                                                        <  1:4556932/150‑1 (MQ=255)
                                                 aaaaTCACGCAATCCATACAAAAAAACCAGATTTTCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAATATGCCATGTGTAGCCAGACGAACACTTTCGGTGGTCTTAAAAAACTAAATATCTTATTTTAACTATTTGttt                                                                                                      <  2:4286158/150‑1 (MQ=255)
                                                   aaTCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAATACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTta                                                                                                    <  2:822036/150‑1 (MQ=255)
                                                    aTCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTtat                                                                                                   <  1:3488125/150‑1 (MQ=255)
                                                    aTCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTtat                                                                                                   <  1:4835476/150‑1 (MQ=255)
                                                     tCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTtata                                                                                                  >  1:1855520/1‑150 (MQ=255)
                                                     tCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTtata                                                                                                  <  1:3206338/150‑1 (MQ=255)
                                                       aCGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTataaa                                                                                                <  1:2983859/150‑1 (MQ=255)
                                                           aaTCCATACCACACACCCAGATCTGCACCTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAaataat                                                                                            <  1:2661449/150‑1 (MQ=255)
                                                           aaTCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAaataat                                                                                            >  1:536525/1‑150 (MQ=255)
                                                            aTCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAAtt                                                                                           >  1:755718/1‑150 (MQ=255)
                                                            aTCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAAtt                                                                                           >  1:3820989/1‑150 (MQ=255)
                                                               cATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATttat                                                                                        >  1:3617929/1‑150 (MQ=255)
                                                                  acaacaAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTaa                                                                                     <  1:1088747/150‑1 (MQ=255)
                                                                   caacaaAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAag                                                                                    >  1:1731160/1‑150 (MQ=255)
                                                                   caacaaAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAag                                                                                    <  1:1762209/150‑1 (MQ=255)
                                                                   caacaaAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAag                                                                                    <  2:3362659/150‑1 (MQ=255)
                                                                    aacaaAACCAGATTTGCAATTCGTGTCACAAAATATTTCTATCTTTTTCTAAGTGGTAGATGCCTTGTGAAGCCAGACTAACACTTGCGGTTGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTTaga                                                                                   <  1:1553185/150‑1 (MQ=255)
                                                                    aacaaAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAaga                                                                                   >  1:957813/1‑150 (MQ=255)
                                                                      caaAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGt                                                                                 >  2:2319336/1‑150 (MQ=255)
                                                                       aaaaCCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTc                                                                                <  2:4977108/150‑1 (MQ=255)
                                                                        aaaCCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTACTTGTTTTATAAATAATTTATTAAGAGtat                                                                               >  2:1451081/1‑148 (MQ=255)
                                                                         aaCCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAATACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTCAGAGTCTa                                                                              >  2:92678/1‑150 (MQ=255)
                                                                         aaCCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAATACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTa                                                                              >  2:347608/1‑150 (MQ=255)
                                                                         aaCCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTa                                                                              <  2:2963593/150‑1 (MQ=255)
                                                                         aaCCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTa                                                                              <  2:1625054/150‑1 (MQ=255)
                                                                           ccAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTaaa                                                                            <  2:160677/150‑1 (MQ=255)
                                                                           ccAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTaaa                                                                            <  2:4427146/150‑1 (MQ=255)
                                                                           ccAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTaaa                                                                            >  2:2049424/1‑150 (MQ=255)
                                                                           ccAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTGTTATAAATAATTTCTTCAGAGTCTaga                                                                            >  2:1132716/1‑148 (MQ=255)
                                                                             aGATTTGCAGTTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACa                                                                          >  2:3893161/1‑150 (MQ=255)
                                                                             aGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACa                                                                          >  1:3361855/1‑150 (MQ=255)
                                                                                 ttGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAggg                                                                      <  1:4330272/150‑1 (MQ=255)
                                                                                    cAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAg                                                                   >  1:500163/1‑150 (MQ=255)
                                                                                       ttCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCtt                                                                >  1:79630/1‑150 (MQ=255)
                                                                                        tCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCttt                                                               <  2:2195891/150‑1 (MQ=255)
                                                                                          gtgaCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGc                                                             <  1:3527107/150‑1 (MQ=255)
                                                                                           tgtCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCa                                                            >  1:472929/1‑150 (MQ=255)
                                                                                            gtCACAAAATATGTCGATCTTGTTATAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAACACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGCGTCTAAACCAGGGGAGCTTTGCaa                                                           <  2:47195/150‑1 (MQ=255)
                                                                                             tCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTATAAACAAGGGGAGCTTTGCAAg                                                          >  2:733086/1‑150 (MQ=255)
                                                                                             tCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCATCAAAAAATAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAg                                                          >  2:4648854/1‑150 (MQ=255)
                                                                                                       tGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGt                                                <  2:2863936/150‑1 (MQ=255)
                                                                                                            aTCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCt                                           >  2:4615121/1‑150 (MQ=255)
                                                                                                             tCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCtt                                          <  1:4214045/150‑1 (MQ=255)
                                                                                                              cTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTg                                         <  1:3711553/150‑1 (MQ=255)
                                                                                                                   tCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTACGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGTTGaa                                    >  2:355455/1‑150 (MQ=255)
                                                                                                                     tAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGaaaa                                  >  2:264/1‑150 (MQ=255)
                                                                                                                     tAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGaaaa                                  <  1:160937/150‑1 (MQ=255)
                                                                                                                     tAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGaaaa                                  <  1:2269292/150‑1 (MQ=255)
                                                                                                                       agagGAAGATGCCATGTGAAGCCAGACGAATACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATc                                <  2:1372567/150‑1 (MQ=255)
                                                                                                                         cggAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAg                              <  1:4854640/149‑1 (MQ=255)
                                                                                                                          ggAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAgt                             >  1:4156618/1‑150 (MQ=255)
                                                                                                                               aTGCCATGTTAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTATTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTa                        <  2:2084119/150‑1 (MQ=255)
                                                                                                                               aTGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATAAGTGTTTa                        >  2:2070478/1‑150 (MQ=255)
                                                                                                                                   cATGAGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATATGTTTTATAAATAAATTAAAAAGAGTCTAAACAAGGGGAGCATATCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTAccag                    >  2:384701/1‑148 (MQ=255)
                                                                                                                                    aTGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGc                   >  1:3054792/1‑150 (MQ=255)
                                                                                                                                         aaGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCa              >  1:663617/1‑150 (MQ=255)
                                                                                                                                         aaGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCa              >  1:4303463/1‑150 (MQ=255)
                                                                                                                                           tccAGAAGACGACTTGCTGTTGTCTTAAAATACTCCAGAACTGACTTTAACTATTTGTTTTATAAATAATTTATTTAGAGTCTAAACCAGGGGAGCTTTGCAAGCTAACTCAGTGCGCTTGGTGAAAATCAGTGTTTTCGCGCCATCAgg            <  2:3240105/149‑1 (MQ=255)
                                                                                                                                            ccAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAATAGTCTAAACAAGGGGAGCTTTGCAAGCTAACACAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGc           >  1:516135/1‑150 (MQ=255)
                                                                                                                                             cAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCt          <  1:2764373/150‑1 (MQ=255)
                                                                                                                                              aGACGAACACTTGCGGAGGTCTTCAAAAAATAAAGATCTTAGTTTAACTATTTGTTTTATAAAAAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGATAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTg         >  1:2295855/1‑150 (MQ=255)
                                                                                                                                                aCGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAg       <  1:3860374/150‑1 (MQ=255)
                                                                                                                                                aCGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAg       <  1:1670926/150‑1 (MQ=255)
                                                                                                                                                 cGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGTGc      >  1:4927713/1‑150 (MQ=255)
                                                                                                                                                 cGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGc      <  2:2019749/150‑1 (MQ=255)
                                                                                                                                                 cGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGc      <  2:2505067/150‑1 (MQ=255)
                                                                                                                                                  gAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCa     >  1:2389299/1‑150 (MQ=255)
                                                                                                                                                   aaCACTTGCGGTTGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCAt    <  2:1892133/150‑1 (MQ=255)
                                                                                                                                                   aaCACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCAt    <  2:3745222/150‑1 (MQ=255)
                                                                                                                                                   aaCACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCAt    <  1:3950624/150‑1 (MQ=255)
                                                                                                                                                   aaCACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCAt    >  1:2908755/1‑150 (MQ=255)
                                                                                                                                                    acacTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATa   >  2:3986799/1‑150 (MQ=255)
                                                                                                                                                     cacTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATaa  <  2:273416/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAA  >  CP009273/4529717‑4530015

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: