Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,376,178 C→A 100% E1004* (GAA→TAA)  mscM ← mechanosensitive channel protein, miniconductance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,376,1780CA98.4% 590.5 / 0.6 187E1004* (GAA→TAA) mscMmechanosensitive channel protein, miniconductance
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base A (90/93);  total (90/97)
Fisher's exact test for biased strand distribution p-value = 2.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 193 total aligned reads displayed.

CGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGATGGTGGTGG  >  CP009273/4376033‑4376321
                                                                                                                                                 |                                                                                                                                               
cGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGc                                                                                                                                               <  1:3295162/148‑1 (MQ=255)
 ggTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGGGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCt                                                                                                                                              >  2:1779218/1‑148 (MQ=255)
  gTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTa                                                                                                                                             <  1:2331843/148‑1 (MQ=255)
   tGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTAt                                                                                                                                            <  2:581364/148‑1 (MQ=255)
    gACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTAtt                                                                                                                                           <  1:1878573/148‑1 (MQ=255)
      cccATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTgg                                                                                                                                         >  1:360052/1‑148 (MQ=255)
      cccATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTgg                                                                                                                                         >  1:360085/1‑148 (MQ=255)
       ccATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGc                                                                                                                                        <  1:1610980/148‑1 (MQ=255)
       ccATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGc                                                                                                                                        >  1:170036/1‑148 (MQ=255)
        cATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCa                                                                                                                                       >  2:743953/1‑148 (MQ=255)
        cATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCa                                                                                                                                       >  2:2745606/1‑148 (MQ=255)
        cATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCa                                                                                                                                       >  2:2745271/1‑148 (MQ=255)
        cATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCa                                                                                                                                       >  2:2504199/1‑148 (MQ=255)
          tctcAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATc                                                                                                                                     >  2:3183242/1‑148 (MQ=255)
             cAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGc                                                                                                                                  <  2:2280086/148‑1 (MQ=255)
               gcggcgTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAg                                                                                                                                <  1:322485/148‑1 (MQ=255)
                   cgTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGc                                                                                                                            <  2:1814183/148‑1 (MQ=255)
                      aaaTACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCgg                                                                                                                         >  2:3493587/1‑148 (MQ=255)
                        aTACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTa                                                                                                                       >  2:2841802/1‑148 (MQ=255)
                        aTACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTa                                                                                                                       >  1:2949089/1‑148 (MQ=255)
                        aTACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTa                                                                                                                       >  1:283966/1‑148 (MQ=255)
                         tccGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTAt                                                                                                                      <  2:1043988/146‑1 (MQ=255)
                         tACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTAt                                                                                                                      >  1:3029408/1‑148 (MQ=255)
                         tACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTAt                                                                                                                      <  1:2700809/148‑1 (MQ=255)
                            gCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATAGt                                                                                                                   <  2:2653686/148‑1 (MQ=255)
                                    aaaaTCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAAcacca                                                                                                           <  1:3028144/148‑1 (MQ=255)
                                     aaaTCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAAcaccac                                                                                                          >  1:2457159/1‑148 (MQ=255)
                                         cTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCAcgcgc                                                                                                      >  2:2885670/1‑148 (MQ=255)
                                          tGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCAcgcgcg                                                                                                     <  1:921882/148‑1 (MQ=255)
                                            aaTCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTg                                                                                                   >  1:1767990/1‑148 (MQ=255)
                                              tCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGAc                                                                                                 >  1:523051/1‑148 (MQ=255)
                                                 ccTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGa                                                                                              >  1:1345452/1‑148 (MQ=255)
                                                       gCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCaga                                                                                        <  1:3213955/148‑1 (MQ=255)
                                                        cAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGACTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCagag                                                                                       <  1:2071833/148‑1 (MQ=255)
                                                         aGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCagaga                                                                                      <  2:594904/148‑1 (MQ=255)
                                                          gATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCagagag                                                                                     <  1:3008554/148‑1 (MQ=255)
                                                          gATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCagagag                                                                                     >  1:1144819/1‑148 (MQ=255)
                                                             ccaccaGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGa                                                                                  >  2:3089307/1‑148 (MQ=255)
                                                             ccaccaGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGa                                                                                  <  2:1600094/148‑1 (MQ=255)
                                                              caccaGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGAc                                                                                 <  2:1331647/148‑1 (MQ=255)
                                                              caccaGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGAc                                                                                 <  1:3427696/148‑1 (MQ=255)
                                                               accaGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTTTTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGAcc                                                                                <  1:1122500/148‑1 (MQ=255)
                                                                   ggAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGtt                                                                            >  1:715185/1‑148 (MQ=255)
                                                                    gAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTg                                                                           >  1:356737/1‑148 (MQ=255)
                                                                    gAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTg                                                                           >  1:1308165/1‑148 (MQ=255)
                                                                    gAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTg                                                                           <  1:2178387/148‑1 (MQ=255)
                                                                     aaGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGa                                                                          <  1:764910/148‑1 (MQ=255)
                                                                      aGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGAt                                                                         <  1:1181247/148‑1 (MQ=255)
                                                                      aGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGAt                                                                         <  1:1475398/148‑1 (MQ=255)
                                                                      aGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGAt                                                                         <  2:804727/148‑1 (MQ=255)
                                                                       gACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATa                                                                        >  2:1320605/1‑148 (MQ=255)
                                                                       gACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCATGACCGAGTCAGAGAGCGACCAGTTGATa                                                                        >  1:1666658/1‑148 (MQ=255)
                                                                        aCTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATaa                                                                       <  2:2624141/148‑1 (MQ=255)
                                                                         cTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCTACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATaaa                                                                      >  2:264542/1‑148 (MQ=255)
                                                                         cTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATaaa                                                                      >  2:264877/1‑148 (MQ=255)
                                                                         cTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATaaa                                                                      <  2:2792038/148‑1 (MQ=255)
                                                                           tCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACt                                                                    >  2:1077278/1‑148 (MQ=255)
                                                                            ccGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTg                                                                   <  2:2531259/148‑1 (MQ=255)
                                                                                tGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCg                                                               >  2:1360619/1‑148 (MQ=255)
                                                                                tGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCg                                                               <  2:838732/148‑1 (MQ=255)
                                                                                 gCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCgg                                                              <  2:2447095/148‑1 (MQ=255)
                                                                                  ccGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGAAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGt                                                             >  1:1826470/1‑148 (MQ=255)
                                                                                     ggTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGtaat                                                          >  2:2873021/1‑148 (MQ=255)
                                                                                        tGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATaaa                                                       <  1:2002699/148‑1 (MQ=255)
                                                                                         gTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAAc                                                      <  2:1883493/148‑1 (MQ=255)
                                                                                         gTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAAc                                                      <  2:2033798/148‑1 (MQ=255)
                                                                                          tCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACg                                                     >  2:877099/1‑148 (MQ=255)
                                                                                             aTCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCt                                                  <  2:3107524/148‑1 (MQ=255)
                                                                                               ccccAGCGTGCCGCGACGCGCTGCGTTGAGCAGTTTTTCCGTAAATTCTTTGCTATTGGCTTCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGCCCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTg                                                <  1:1214215/146‑1 (MQ=255)
                                                                                               cACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTg                                                >  1:3063586/1‑148 (MQ=255)
                                                                                               cACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTg                                                >  2:1435802/1‑148 (MQ=255)
                                                                                                aCCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGt                                               <  1:1039468/148‑1 (MQ=255)
                                                                                                aCCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGt                                               <  1:1034833/148‑1 (MQ=255)
                                                                                                aCCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGt                                               <  1:3138319/148‑1 (MQ=255)
                                                                                                        gcagcGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACg                                       <  2:1678846/148‑1 (MQ=255)
                                                                                                         cagcGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGa                                      <  1:274564/148‑1 (MQ=255)
                                                                                                         cagcGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGa                                      <  1:274470/148‑1 (MQ=255)
                                                                                                          agcGACGCGCTGCGGTGAGCAAGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGAt                                     <  2:2566168/148‑1 (MQ=255)
                                                                                                           gcGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGAtt                                    >  1:2038235/1‑148 (MQ=255)
                                                                                                            cGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTa                                   <  2:144405/148‑1 (MQ=255)
                                                                                                            cGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTa                                   >  2:443897/1‑148 (MQ=255)
                                                                                                                  gcTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCtt                             <  2:1687980/148‑1 (MQ=255)
                                                                                                                   cTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCttt                            <  2:487682/148‑1 (MQ=255)
                                                                                                                   cTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATGGGCATCGGCAGGGGCCGGCATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCttt                            <  1:2660097/148‑1 (MQ=255)
                                                                                                                    tGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTa                           <  1:2098255/148‑1 (MQ=255)
                                                                                                                     gCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTAc                          <  2:2588871/148‑1 (MQ=255)
                                                                                                                      cGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACg                         >  1:2127853/1‑148 (MQ=255)
                                                                                                                                gATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCg               >  2:936829/1‑148 (MQ=255)
                                                                                                                                gATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCg               >  2:683270/1‑148 (MQ=255)
                                                                                                                                gATTTCCGTCACTTCTTAGCTATCGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCg               <  2:2161447/148‑1 (MQ=255)
                                                                                                                                   ttCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTg            >  2:1954170/1‑148 (MQ=255)
                                                                                                                                     ccGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGAt          >  2:638114/1‑148 (MQ=255)
                                                                                                                                        tCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAACCGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGAtggt       <  2:3063876/148‑1 (MQ=255)
                                                                                                                                          aCTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGAtggtgg     <  2:1678109/148‑1 (MQ=255)
                                                                                                                                           cttcttAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGAtggtggt    >  2:1973951/1‑148 (MQ=255)
                                                                                                                                           cttcttAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGAtggtggt    >  2:1235036/1‑148 (MQ=255)
                                                                                                                                           cttcttAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGAtggtggt    >  1:1918000/1‑148 (MQ=255)
                                                                                                                                            ttcttaGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGAtggtggtg   >  2:3103238/1‑148 (MQ=255)
                                                                                                                                            ttcttaGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGAtggtggtg   >  1:2586993/1‑148 (MQ=255)
                                                                                                                                             tcttaGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGAtggtggtgg  >  1:1511371/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGATGGTGGTGG  >  CP009273/4376033‑4376321

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: