Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,933,379 G→T 5.9% I198I (ATC→ATA hsrA ← putative multidrug or homocysteine efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,933,3790GT 5.9% 865.9 / 17.8 228I198I (ATC→ATAhsrAputative multidrug or homocysteine efflux system
Reads supporting (aligned to +/- strand):  ref base G (99/115);  new base T (7/7);  total (106/122)
Fisher's exact test for biased strand distribution p-value = 7.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 231 total aligned reads displayed.

TCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTGATTAAAAATATCCAGTGCCAGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACGGGGCCTAA  >  CP009273/3933232‑3933525
                                                                                                                                                   |                                                                                                                                                  
tcgtGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTg                                                                                                                                                    >  1:359237/1‑148 (MQ=255)
tcgtGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTg                                                                                                                                                    >  2:1600899/1‑148 (MQ=255)
 cgtGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGa                                                                                                                                                   >  1:1920940/1‑148 (MQ=255)
 cgtGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGa                                                                                                                                                   >  2:2980632/1‑148 (MQ=255)
       tGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGa                                                                                                                                             <  1:3049553/148‑1 (MQ=255)
        ggAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTTATATCGaa                                                                                                                                            >  1:1307287/1‑148 (MQ=255)
        ggAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGaa                                                                                                                                            >  2:933521/1‑148 (MQ=255)
            aaTATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTg                                                                                                                                        <  1:65684/148‑1 (MQ=255)
            aaTATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTg                                                                                                                                        >  1:927783/1‑148 (MQ=255)
            aaTATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTg                                                                                                                                        <  2:762591/148‑1 (MQ=255)
             atatAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGc                                                                                                                                       >  1:1007487/1‑148 (MQ=255)
             atatAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGc                                                                                                                                       >  1:1252182/1‑148 (MQ=255)
              tataGAGAAGCAGTAACCCGATGCTGGTGTCAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGcg                                                                                                                                      <  2:3303101/148‑1 (MQ=255)
              tataGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGcg                                                                                                                                      <  1:816005/148‑1 (MQ=255)
                 agagAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGcgtcg                                                                                                                                   <  2:3011705/148‑1 (MQ=255)
                      aGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTTATATCGAATCTGCGTCGTGCgg                                                                                                                              <  2:1578466/148‑1 (MQ=255)
                      aGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCgg                                                                                                                              <  2:2844459/148‑1 (MQ=255)
                        cAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCggtg                                                                                                                            >  2:3234488/1‑148 (MQ=255)
                         aGTAACCCGTTGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGTAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCggtgg                                                                                                                           <  2:349108/148‑1 (MQ=255)
                           tAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGgtggtg                                                                                                                         <  1:1207990/148‑1 (MQ=255)
                              cccGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGaaa                                                                                                                      >  1:1601501/1‑148 (MQ=255)
                              cccGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGaaa                                                                                                                      <  2:1602067/148‑1 (MQ=255)
                                   tGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTAGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTggg                                                                                                                 >  2:3146984/1‑148 (MQ=255)
                                     cTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCa                                                                                                               <  2:2614277/148‑1 (MQ=255)
                                      tGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCat                                                                                                              <  1:1031920/148‑1 (MQ=255)
                                       ggTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCata                                                                                                             >  1:1081412/1‑148 (MQ=255)
                                        gTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCatat                                                                                                            <  1:732518/148‑1 (MQ=255)
                                         tGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATg                                                                                                           >  2:1590514/1‑148 (MQ=255)
                                           aCAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGtt                                                                                                         <  2:1894651/148‑1 (MQ=255)
                                             aaTTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTg                                                                                                       <  2:1533447/148‑1 (MQ=255)
                                              aTTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTgc                                                                                                      >  1:3148596/1‑148 (MQ=255)
                                               ttACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTgcg                                                                                                     >  1:203942/1‑148 (MQ=255)
                                                   cGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGcgcgcg                                                                                                 <  2:648487/148‑1 (MQ=255)
                                                       aaGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTaaa                                                                                             <  1:1450897/148‑1 (MQ=255)
                                                       aaGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTaaa                                                                                             <  1:3090155/148‑1 (MQ=255)
                                                         ggCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAaaga                                                                                           <  1:1107242/148‑1 (MQ=255)
                                                            aaTCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAgg                                                                                        <  1:1659905/148‑1 (MQ=255)
                                                              tCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGcc                                                                                      >  2:87337/1‑148 (MQ=255)
                                                              tCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGcc                                                                                      >  2:1861952/1‑148 (MQ=255)
                                                               ccAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGccc                                                                                     >  1:1824464/1‑148 (MQ=255)
                                                                cAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCg                                                                                    <  2:203942/148‑1 (MQ=255)
                                                                 aGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGc                                                                                   >  1:722618/1‑148 (MQ=255)
                                                                  gCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCa                                                                                  <  1:2911426/148‑1 (MQ=255)
                                                                  gCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCa                                                                                  >  1:144678/1‑148 (MQ=255)
                                                                        gACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCt                                                                            <  2:2289924/148‑1 (MQ=255)
                                                                         aCAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTa                                                                           >  1:1065485/1‑148 (MQ=255)
                                                                         aCAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTa                                                                           >  1:1887046/1‑148 (MQ=255)
                                                                         aCAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTa                                                                           <  1:1942845/148‑1 (MQ=255)
                                                                          cAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTAt                                                                          <  1:3050232/148‑1 (MQ=255)
                                                                          aaaTGTTTTCCCCGAATAGCTATATTCCTCTTGAGAAGAGAACAAGGGTGAGGCCAAACAGCAAACAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGACCGCAATACCTAt                                                                          <  1:1406190/147‑1 (MQ=255)
                                                                           aaTCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATg                                                                         <  2:2741493/148‑1 (MQ=255)
                                                                           aaTCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATg                                                                         <  2:2183904/148‑1 (MQ=255)
                                                                           aaTCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATg                                                                         >  2:1434647/1‑148 (MQ=255)
                                                                            aTCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATgg                                                                        >  2:194421/1‑148 (MQ=255)
                                                                             tCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATggg                                                                       >  1:561369/1‑148 (MQ=255)
                                                                                tttCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGat                                                                    >  2:741692/1‑148 (MQ=255)
                                                                                tttCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGat                                                                    <  1:265938/148‑1 (MQ=255)
                                                                                    cccGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTg                                                                >  1:698570/1‑148 (MQ=255)
                                                                                    cccGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTg                                                                <  1:284256/148‑1 (MQ=255)
                                                                                    cccGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTg                                                                <  2:2265769/148‑1 (MQ=255)
                                                                                    cccGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTg                                                                <  2:1335022/148‑1 (MQ=255)
                                                                                        aaTAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTGATTa                                                            >  1:3163180/1‑148 (MQ=255)
                                                                                         aTAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTGATTaa                                                           <  1:257858/148‑1 (MQ=255)
                                                                                         aTAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTGATTaa                                                           >  2:2072401/1‑148 (MQ=255)
                                                                                          tAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTGATTaaa                                                          >  2:1473829/1‑148 (MQ=255)
                                                                                          tAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTGATTaaa                                                          >  1:1916425/1‑148 (MQ=255)
                                                                                            gCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTGATTaaaaa                                                        <  2:1010983/148‑1 (MQ=255)
                                                                                               ctATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTGATTAAAAatat                                                     <  2:1252182/148‑1 (MQ=255)
                                                                                               ctATTCCGCTTGAGAAGAGAACAACGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTGATTAAAAatat                                                     <  1:247550/148‑1 (MQ=255)
                                                                                                 aTTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTGATTAAAAATATcc                                                   >  1:169990/1‑148 (MQ=255)
                                                                                                 aTTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTGATTAAAAATATcc                                                   <  2:2585100/148‑1 (MQ=255)
                                                                                                   tCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTGATTAAAAATATCCAg                                                 <  1:2430051/148‑1 (MQ=255)
                                                                                                   tCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTGATTAAAAATATCCAg                                                 >  1:1418690/1‑148 (MQ=255)
                                                                                                   tCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTGATTAAAAATATCCAg                                                 <  1:860888/148‑1 (MQ=255)
                                                                                                   tCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTGATTAAAAATATCCAg                                                 >  2:2508634/1‑148 (MQ=255)
                                                                                                     cGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTGATTAAAAATATCCAGTg                                               <  1:1813849/148‑1 (MQ=255)
                                                                                                     cGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTGATTAAAAATATCCAGTg                                               <  1:476589/148‑1 (MQ=255)
                                                                                                      gCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTGATTAAAAATATCCAGTGc                                              <  1:1600899/148‑1 (MQ=255)
                                                                                                       cTTGAGAAGAGAACACGGCTGAGGCCAAACAGCAAAAAGCCAGTTATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTGATTAAAAATATCCAGTGcc                                             <  2:2671732/148‑1 (MQ=255)
                                                                                                        ttGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTGATTAAAAATATCCAGTGCCa                                            <  2:719495/148‑1 (MQ=255)
                                                                                                         tGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTGATTAAAAATATCCAGTGCCAg                                           >  1:2819815/1‑148 (MQ=255)
                                                                                                         tGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTGATTAAAAATATCCAGTGCCAg                                           <  2:2352190/148‑1 (MQ=255)
                                                                                                         tGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTGATTAAAAATATCCAGTGCCAg                                           >  1:1983682/1‑148 (MQ=255)
                                                                                                           agaagagaaCAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTGATTAAAAATATCAAGTGCCAGGt                                         >  1:575375/1‑148 (MQ=255)
                                                                                                                  aaCAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTGATTAAAAATATCCAGTGCCAGGTTGCCCAg                                  <  1:1862105/148‑1 (MQ=255)
                                                                                                                   aCAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTGATTAAAAATATCCAGTGCCAGGTTGCCCAgg                                 <  1:2247175/148‑1 (MQ=255)
                                                                                                                    cAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTGATTAAAAATATCCAGTGCCAGGTTGCCCAGGt                                >  2:1006192/1‑148 (MQ=255)
                                                                                                                                  aaCAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTGATTAAAAATATCCAGTGCCAGGTTGCCCAGGTGACCAGCAcgccgc                  >  2:2215575/1‑148 (MQ=255)
                                                                                                                                   aCAGCAAAAAGCCAGTTATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTGATTAAAAATATCCAGTGCCAGGTTGCCCAGGTGACCAGCACgccgcc                 >  2:917652/1‑148 (MQ=255)
                                                                                                                                   aCAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGGGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGCCCCGCAATACCTATGGGGATATTGATTAAAAATATCCAGTGCCAGGTTGCCCAGGTGACCAGCACgccgcc                 <  2:410932/148‑1 (MQ=255)
                                                                                                                                      gCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTGATTAAAAATATCCAGTGCCAGGTTGCCCAGGTGACCAGCACGCCGCCAAg              >  1:2209159/1‑148 (MQ=255)
                                                                                                                                       cAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTGATTAAAAATATCCAGTGCCAGGTTGCCCAGGTGACCAGCACGCCGCCAAGa             <  2:2101103/148‑1 (MQ=255)
                                                                                                                                       cAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTGATTAAAAATATCCAGTGCCAGGTTGCCCAGGTGACCAGCACGCCGCCAAGa             <  2:3102313/148‑1 (MQ=255)
                                                                                                                                         aaaaGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTGATTAAAAATATCCAGTGCCAGGTTGCCCAGGTGACCAGCACGCCGCCAAGAAc           >  1:3061021/1‑148 (MQ=255)
                                                                                                                                          aaaGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTGATTAAAAATATCCAGTGCCAGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACg          <  2:1823313/148‑1 (MQ=255)
                                                                                                                                             gCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTGATTAAAAATATCCAGTGCCAGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACgggg       >  1:1210015/1‑148 (MQ=255)
                                                                                                                                             gCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTGATTAAAAATATCCAGTGCCAGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACgggg       >  1:2341660/1‑148 (MQ=255)
                                                                                                                                             gCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTGATTAAAAATATCCAGTGCCAGATTGCCCAGGTGACCAGCACGCCGCCAAGAACgggg       >  1:438848/1‑148 (MQ=255)
                                                                                                                                                aGTTATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTGATTAAAAATATCCAGTGCCAGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACGGGGCCt    >  2:625858/1‑148 (MQ=255)
                                                                                                                                                  tGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTGATTAAAAATATCCAGTGCCAGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACGGGGCCTaa  >  2:2026475/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGGCTGAGGCCAAACAGCAAAAAGCCAGTGATATCGAATCTGCGTCGTGCGGTGGTGAAATTGGGCATATGTTTGCGCGCGTAAAGAAGGCCCGCAATACCTATGGGGATATTGATTAAAAATATCCAGTGCCAGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACGGGGCCTAA  >  CP009273/3933232‑3933525

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: