Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,376,178 C→A 92.7% 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,1780CA92.7% 617.6 / 40.1 228E1004* (GAA→TAA) mscMmechanosensitive channel protein, miniconductance
Reads supporting (aligned to +/- strand):  ref base C (14/3);  new base A (103/108);  total (117/111)
Fisher's exact test for biased strand distribution p-value = 1.02e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

Only 100 of 231 total aligned reads displayed.

GGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGATGGTGGTG  >  CP009273/4376034‑4376320
                                                                                                                                                |                                                                                                                                              
ggTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCt                                                                                                                                             >  1:1872706/1‑148 (MQ=255)
 gTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTa                                                                                                                                            <  2:563436/148‑1 (MQ=255)
  tGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTAt                                                                                                                                           <  2:993579/148‑1 (MQ=255)
  tGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTAt                                                                                                                                           >  1:1783851/1‑148 (MQ=255)
      ccATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGc                                                                                                                                       <  1:3168348/148‑1 (MQ=255)
       cATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCa                                                                                                                                      >  1:44821/1‑148 (MQ=255)
       cATCTCAGCGGAGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCTGGATTTCCGTCACTTCTTAGCTATTGGCa                                                                                                                                      <  2:1405798/148‑1 (MQ=255)
        aTCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCAt                                                                                                                                     <  1:2849967/148‑1 (MQ=255)
         tctcAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATc                                                                                                                                    >  1:650462/1‑148 (MQ=255)
          ctcAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCg                                                                                                                                   >  1:1991524/1‑148 (MQ=255)
              gcggcgTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAg                                                                                                                               <  2:2979388/148‑1 (MQ=255)
                  cgTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGc                                                                                                                           >  1:1898169/1‑148 (MQ=255)
                  cgTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGc                                                                                                                           <  1:3511302/148‑1 (MQ=255)
                    tAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCg                                                                                                                         >  1:3127429/1‑148 (MQ=255)
                    tAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCg                                                                                                                         <  2:3549780/148‑1 (MQ=255)
                     aaaTACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCgg                                                                                                                        >  2:1260839/1‑148 (MQ=255)
                     aaaTACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCgg                                                                                                                        >  1:291062/1‑148 (MQ=255)
                          cGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCg                                                                                                                   >  1:2891216/1‑148 (MQ=255)
                            cAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTc                                                                                                                 >  2:131550/1‑148 (MQ=255)
                              gCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCaa                                                                                                               >  1:48956/1‑148 (MQ=255)
                              gCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCaa                                                                                                               >  2:2534444/1‑148 (MQ=255)
                                tCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAaca                                                                                                             <  1:315344/148‑1 (MQ=255)
                                 cGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAacac                                                                                                            <  2:2161963/148‑1 (MQ=255)
                                  gAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAAcacc                                                                                                           >  2:1048297/1‑148 (MQ=255)
                                     aaTCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCAcg                                                                                                        <  1:306757/148‑1 (MQ=255)
                                       tCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCAcgcg                                                                                                      >  1:2370074/1‑148 (MQ=255)
                                        cTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCAcgcgc                                                                                                     >  1:3478876/1‑148 (MQ=255)
                                        cTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCAcgcgc                                                                                                     <  1:2114188/148‑1 (MQ=255)
                                           aaTCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTg                                                                                                  >  1:2313169/1‑148 (MQ=255)
                                           aaTCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTg                                                                                                  >  2:38282/1‑148 (MQ=255)
                                              ccccTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGAcc                                                                                               >  2:2266307/1‑148 (MQ=255)
                                                   gTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTc                                                                                          >  1:2167052/1‑148 (MQ=255)
                                                   gTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTc                                                                                          >  2:2123543/1‑148 (MQ=255)
                                                     tGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCag                                                                                        <  1:978764/148‑1 (MQ=255)
                                                     tGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCag                                                                                        <  2:2329694/148‑1 (MQ=255)
                                                     tGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCag                                                                                        >  1:1623023/1‑148 (MQ=255)
                                                       cAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCagag                                                                                      <  2:1872724/148‑1 (MQ=255)
                                                        aGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCagaga                                                                                     >  2:2802949/1‑148 (MQ=255)
                                                          aTCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGc                                                                                   >  2:2596686/1‑148 (MQ=255)
                                                          aTCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGc                                                                                   <  1:2199325/148‑1 (MQ=255)
                                                            ccaccaGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGa                                                                                 <  1:246112/148‑1 (MQ=255)
                                                             caccaGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGAc                                                                                <  1:1038496/148‑1 (MQ=255)
                                                              accaGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGAcc                                                                               >  2:2441308/1‑148 (MQ=255)
                                                               ccaGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCa                                                                              >  1:1939701/1‑148 (MQ=255)
                                                                caGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCTGGGGCCGGTATCGTCAACACCACGCGCGTGACCTAGTCAGAGAGCGACCAg                                                                             <  1:1644505/148‑1 (MQ=255)
                                                                 aGGAAGACTTGGGGTGCGGGGTTGTCGATCACCAGCGAGCAGCGACGCGGTGGGTTGAGCTGGATTTCCGTCACTTATTAGATATTGGAATCGGCTGGGGCCGGTATCTTAAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGt                                                                            <  1:1231018/148‑1 (MQ=255)
                                                                  ggAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGtt                                                                           >  1:1864369/1‑148 (MQ=255)
                                                                   gAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTg                                                                          >  1:3432397/1‑148 (MQ=255)
                                                                   gAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTg                                                                          <  1:2811155/148‑1 (MQ=255)
                                                                   gAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTg                                                                          >  1:1718953/1‑148 (MQ=255)
                                                                   gAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTg                                                                          >  1:2036359/1‑148 (MQ=255)
                                                                   gAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTg                                                                          <  2:3127429/148‑1 (MQ=255)
                                                                      gACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCTGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATa                                                                       >  1:3275689/1‑148 (MQ=255)
                                                                            cGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCCCGCGCGGGACCGAGTCAGAGAGCGACCAGTTGATAAACTGc                                                                 <  2:1205979/148‑1 (MQ=255)
                                                                            cGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGc                                                                 <  2:3028613/148‑1 (MQ=255)
                                                                                     gTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGtaata                                                        <  1:565716/148‑1 (MQ=255)
                                                                                      ttGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGtaataa                                                       <  1:1266106/148‑1 (MQ=255)
                                                                                       tGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATaaa                                                      <  1:304746/148‑1 (MQ=255)
                                                                                       tGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATGTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATaaa                                                      <  2:394867/148‑1 (MQ=255)
                                                                                          cGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGCGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGc                                                   <  1:2687342/148‑1 (MQ=255)
                                                                                          cGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGc                                                   <  2:2987960/148‑1 (MQ=255)
                                                                                           gATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGcc                                                  <  2:650462/148‑1 (MQ=255)
                                                                                           gATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGcc                                                  <  1:366253/148‑1 (MQ=255)
                                                                                           gATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGcc                                                  <  1:3414407/148‑1 (MQ=255)
                                                                                              cACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTg                                               >  1:865800/1‑148 (MQ=255)
                                                                                              cACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTg                                               >  1:820432/1‑148 (MQ=255)
                                                                                               aCCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGt                                              >  2:1371305/1‑148 (MQ=255)
                                                                                               aCCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGt                                              <  2:1783851/148‑1 (MQ=255)
                                                                                                 cAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTc                                            >  2:1681545/1‑148 (MQ=255)
                                                                                                   gCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCgg                                          >  2:1120179/1‑148 (MQ=255)
                                                                                                    cGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGc                                         <  2:1145618/148‑1 (MQ=255)
                                                                                                    cGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGc                                         <  1:1182604/148‑1 (MQ=255)
                                                                                                    cGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGc                                         >  2:2537501/1‑148 (MQ=255)
                                                                                                     gcgcagcGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCATAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCa                                        <  1:2413772/146‑1 (MQ=255)
                                                                                                       gcagcGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACg                                      <  2:3153320/148‑1 (MQ=255)
                                                                                                       gcagcGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGAAACGTAAACACCACCCGCGTGACCGAGTCAGAGAGCGACCACTTGATAAACTGCTCGGAAATAAACGCCTTGTTCGGCACg                                      >  2:3371987/1‑148 (MQ=255)
                                                                                                         agcGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGAt                                    <  1:2204279/148‑1 (MQ=255)
                                                                                                            gACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTAt                                 <  2:2952511/148‑1 (MQ=255)
                                                                                                             aCGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTAtc                                <  1:259422/148‑1 (MQ=255)
                                                                                                             aCGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTAtc                                <  1:2055462/148‑1 (MQ=255)
                                                                                                              cgcgCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTAtct                               >  2:3308458/1‑148 (MQ=255)
                                                                                                              cgcgCTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTAtct                               >  1:1983404/1‑148 (MQ=255)
                                                                                                                  cTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCttt                           <  1:2864666/148‑1 (MQ=255)
                                                                                                                  cTGCGGTGAGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCttt                           <  1:2640334/148‑1 (MQ=255)
                                                                                                                    gCGGTGAGCAGGATTTCCGTCACTTCTGAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTAc                         <  2:264952/148‑1 (MQ=255)
                                                                                                                          aGCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGGGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTccc                   <  1:1397653/148‑1 (MQ=255)
                                                                                                                           gCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCa                  >  2:1530305/1‑148 (MQ=255)
                                                                                                                           gCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCa                  <  2:614086/148‑1 (MQ=255)
                                                                                                                           gCAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCa                  >  2:923157/1‑148 (MQ=255)
                                                                                                                            cAGGATTTCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAg                 <  1:1400264/148‑1 (MQ=255)
                                                                                                                                  ttCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTg           >  1:2971804/1‑148 (MQ=255)
                                                                                                                                  ttCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTg           <  2:1623023/148‑1 (MQ=255)
                                                                                                                                  ttCCGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTg           >  2:837481/1‑148 (MQ=255)
                                                                                                                                    ccGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGAt         <  2:2891216/148‑1 (MQ=255)
                                                                                                                                     cGTCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGAtg        >  1:900772/1‑148 (MQ=255)
                                                                                                                                       tCACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGAtggt      >  1:1449199/1‑148 (MQ=255)
                                                                                                                                        cACTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGAtggtg     <  2:3316417/148‑1 (MQ=255)
                                                                                                                                         aCTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGAtggtgg    >  1:966115/1‑148 (MQ=255)
                                                                                                                                         aCTTCTTAGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGAtggtgg    >  1:420219/1‑148 (MQ=255)
                                                                                                                                           ttcttaGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGAtggtggtg  <  2:1939701/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                              
GGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGCAGATCCACCAGGAAGACTTCCGGTGCCGGGTTGTCGATCACCAGCGAGCAGCGACGCGCTGCGGTGAGCAGGATTTCCGTCACTTCTTCGCTATTGGCATCGGCAGGGGCCGGTATCGTCAACACCACGCGCGTGACCGAGTCAGAGAGCGACCAGTTGATAAACTGCTCGGTAATAAACGCCTTGTTCGGCACGATTATCTCTTTACGGTCCCAGTCGCTGATGGTGGTG  >  CP009273/4376034‑4376320

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: