Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,545,608 C→A 100% R177I (AGA→ATA)  yjiC ← uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,545,6080CA97.8% 859.8 / 8.2 272R177I (AGA→ATA) yjiCuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (3/3);  new base A (128/138);  total (131/141)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 272 total aligned reads displayed.

TTCTATACCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCATAAATAAATTTA  >  CP009273/4545463‑4545752
                                                                                                                                                 |                                                                                                                                                
ttCTATACCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGACTCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATAtt                                                                                                                                                >  1:2369577/1‑148 (MQ=255)
ttCTATACCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATAtt                                                                                                                                                <  2:1729904/148‑1 (MQ=255)
ttCTATACCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATAtt                                                                                                                                                >  1:1751530/1‑148 (MQ=255)
 tCTATACCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATAttt                                                                                                                                               <  1:2753907/148‑1 (MQ=255)
  cTATACCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATAtttt                                                                                                                                              <  2:3108608/148‑1 (MQ=255)
      aCCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGcc                                                                                                                                          <  1:3593318/148‑1 (MQ=255)
           ttattGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTgaga                                                                                                                                     <  2:2820048/148‑1 (MQ=255)
           ttattGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTgaga                                                                                                                                     <  2:2235393/148‑1 (MQ=255)
             attGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAg                                                                                                                                   <  2:335735/148‑1 (MQ=255)
              ttGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGc                                                                                                                                  <  2:241653/148‑1 (MQ=255)
                 ggTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTgg                                                                                                                               >  2:3692640/1‑148 (MQ=255)
                    cATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAg                                                                                                                            <  2:3728185/148‑1 (MQ=255)
                    cATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAg                                                                                                                            <  2:3209841/148‑1 (MQ=255)
                      ttACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAggg                                                                                                                          <  1:3243870/148‑1 (MQ=255)
                         cccAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTg                                                                                                                       >  2:1001764/1‑148 (MQ=255)
                         cccAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTg                                                                                                                       <  2:380718/148‑1 (MQ=255)
                          ccAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGa                                                                                                                      >  2:4540/1‑148 (MQ=255)
                           cAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGaa                                                                                                                     <  1:1961660/148‑1 (MQ=255)
                           cAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGaa                                                                                                                     >  1:338884/1‑148 (MQ=255)
                            aaCTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGaaa                                                                                                                    <  1:3365695/148‑1 (MQ=255)
                             aCTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAat                                                                                                                   >  2:2101745/1‑148 (MQ=255)
                               tCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAatat                                                                                                                 >  2:783449/1‑148 (MQ=255)
                               tCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAatat                                                                                                                 <  2:2870035/148‑1 (MQ=255)
                               tCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAatat                                                                                                                 <  2:3710956/148‑1 (MQ=255)
                                cAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAatata                                                                                                                >  2:899636/1‑148 (MQ=255)
                                cAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAatata                                                                                                                <  2:2588265/148‑1 (MQ=255)
                                  atatTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATc                                                                                                              <  1:748636/148‑1 (MQ=255)
                                  atatTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATc                                                                                                              <  1:2258839/148‑1 (MQ=255)
                                  atatTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATc                                                                                                              <  1:1614249/148‑1 (MQ=255)
                                    atTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCt                                                                                                            >  2:2320717/1‑148 (MQ=255)
                                    atTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCt                                                                                                            >  1:2548024/1‑148 (MQ=255)
                                    atTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCt                                                                                                            >  2:1665153/1‑148 (MQ=255)
                                      tAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTAt                                                                                                          <  1:2695560/148‑1 (MQ=255)
                                        gCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCg                                                                                                        >  1:2962278/1‑148 (MQ=255)
                                         cATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGt                                                                                                       <  1:147644/148‑1 (MQ=255)
                                          aTTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTg                                                                                                      <  1:1309781/148‑1 (MQ=255)
                                           ttAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTgg                                                                                                     >  1:621512/1‑148 (MQ=255)
                                            tAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGt                                                                                                    <  1:2145124/148‑1 (MQ=255)
                                             aaGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGtt                                                                                                   <  2:939554/148‑1 (MQ=255)
                                                      ggCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATaa                                                                                          <  2:3024244/148‑1 (MQ=255)
                                                        cAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAAcc                                                                                        >  2:2757817/1‑148 (MQ=255)
                                                         aaGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCa                                                                                       <  1:2878505/148‑1 (MQ=255)
                                                          aGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAg                                                                                      >  2:503251/1‑148 (MQ=255)
                                                           gggCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAgg                                                                                     >  2:3042653/1‑148 (MQ=255)
                                                            ggCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGa                                                                                    <  1:3623273/148‑1 (MQ=255)
                                                            ggCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGa                                                                                    >  2:1640636/1‑148 (MQ=255)
                                                               gAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACa                                                                                 <  2:643751/148‑1 (MQ=255)
                                                               gAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACa                                                                                 <  2:3420691/148‑1 (MQ=255)
                                                                aaTCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACAt                                                                                <  2:2536608/148‑1 (MQ=255)
                                                                aaTCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACAt                                                                                <  2:1359821/148‑1 (MQ=255)
                                                                  tCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTg                                                                              >  1:3507015/1‑148 (MQ=255)
                                                                  tCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTg                                                                              >  2:1787568/1‑148 (MQ=255)
                                                                   ccGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGc                                                                             >  1:2285661/1‑148 (MQ=255)
                                                                      gaagaaTTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAAc                                                                          >  2:154287/1‑148 (MQ=255)
                                                                       aagaaTTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAAcc                                                                         <  1:2117105/148‑1 (MQ=255)
                                                                       aagaaTTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAAcc                                                                         <  1:3069549/148‑1 (MQ=255)
                                                                         gaaTTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAg                                                                       >  2:2374727/1‑148 (MQ=255)
                                                                            tttGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGc                                                                    <  1:3271149/148‑1 (MQ=255)
                                                                            tttGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGc                                                                    >  2:331775/1‑148 (MQ=255)
                                                                                        gTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTAttt                                                        <  1:721398/148‑1 (MQ=255)
                                                                                         ttACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTc                                                       <  1:3522270/148‑1 (MQ=255)
                                                                                           aCCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCAt                                                     <  2:2623496/148‑1 (MQ=255)
                                                                                             cGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGa                                                   <  2:1908325/148‑1 (MQ=255)
                                                                                              gTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGaa                                                  <  1:2475673/148‑1 (MQ=255)
                                                                                               tCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAAc                                                 <  2:1719721/148‑1 (MQ=255)
                                                                                               tCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAAc                                                 >  1:459079/1‑148 (MQ=255)
                                                                                                cATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACt                                                >  2:1617913/1‑148 (MQ=255)
                                                                                                   tATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGta                                             <  2:1876270/148‑1 (MQ=255)
                                                                                                   tATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGta                                             >  2:2339304/1‑148 (MQ=255)
                                                                                                        gATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTc                                        >  1:2092195/1‑148 (MQ=255)
                                                                                                          tCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCt                                      >  1:3800235/1‑148 (MQ=255)
                                                                                                            aaTATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCa                                    <  1:467674/148‑1 (MQ=255)
                                                                                                             atatAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCaa                                   >  2:1862783/1‑148 (MQ=255)
                                                                                                              tataGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAAc                                  >  2:2222815/1‑148 (MQ=255)
                                                                                                               ataGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACa                                 >  2:2933744/1‑148 (MQ=255)
                                                                                                                 aGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGGACTGTATAGTCCTCAACAgg                               >  2:1320268/1‑148 (MQ=255)
                                                                                                                  gATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGa                              >  2:493474/1‑148 (MQ=255)
                                                                                                                  gATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGa                              <  1:2207315/148‑1 (MQ=255)
                                                                                                                   atatCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGa                              <  2:3611301/147‑1 (MQ=255)
                                                                                                                     atCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCAGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCATCAACAGTATCg                           >  1:1352861/1‑148 (MQ=255)
                                                                                                                      tCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGa                          <  1:805741/148‑1 (MQ=255)
                                                                                                                       cAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGaa                         <  1:3035926/148‑1 (MQ=255)
                                                                                                                       cAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGaa                         <  1:2609215/148‑1 (MQ=255)
                                                                                                                       cAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGaa                         >  1:3085007/1‑148 (MQ=255)
                                                                                                                          tCCTTGGCTGCTCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAgg                      <  2:3494222/148‑1 (MQ=255)
                                                                                                                          tCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAgg                      <  2:3231897/148‑1 (MQ=255)
                                                                                                                          tCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAgg                      >  2:1129358/1‑148 (MQ=255)
                                                                                                                          tCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAgg                      <  1:1854297/148‑1 (MQ=255)
                                                                                                                           ccTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAggg                     <  1:2965274/148‑1 (MQ=255)
                                                                                                                                gctgctCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCATATCGTGGTTCCGGAATAAACATGAAAATAGCAACCAGGTACTAATGTTATTTCATGAACTGTATAGTCCTCAAAATGATCGAAAGGTCTCtt                >  2:2354506/1‑148 (MQ=255)
                                                                                                                                   gctCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCat             >  1:1153895/1‑148 (MQ=255)
                                                                                                                                   gctCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCat             >  2:2842652/1‑148 (MQ=255)
                                                                                                                                    ctCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCata            >  2:2665578/1‑148 (MQ=255)
                                                                                                                                       aTTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCataaat         <  1:1436911/148‑1 (MQ=255)
                                                                                                                                         tGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCataaataa       >  1:2815813/1‑148 (MQ=255)
                                                                                                                                         tGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCataaataa       >  2:3150831/1‑148 (MQ=255)
                                                                                                                                         tGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCataaataa       >  2:2046543/1‑148 (MQ=255)
                                                                                                                                           gTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCATaaataaat     <  1:928556/148‑1 (MQ=255)
                                                                                                                                              cATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCATAAATAAATTTa  >  1:3655613/1‑148 (MQ=255)
                                                                                                                                              cATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCATAAATAAATTTa  >  1:3689957/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TTCTATACCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCATAAATAAATTTA  >  CP009273/4545463‑4545752

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: