Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 3,178,277 G→T E27* (GAA→TAA)  yqiC → BMFP family putative fusogenic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,178,2770GT97.4% 207.0 / ‑3.8 77E27* (GAA→TAA) yqiCBMFP family putative fusogenic protein
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (45/30);  total (45/32)
Fisher's exact test for biased strand distribution p-value = 1.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

TCTGGTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGA  >  CP009273/3178131‑3178422
                                                                                                                                                  |                                                                                                                                                 
tCTGGTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTaa                                                                                                                                                 >  1:590017/1‑149 (MQ=255)
 cTGGTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAg                                                                                                                                                <  2:610769/149‑1 (MQ=255)
  tGGTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACAAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGa                                                                                                                                               >  2:707930/1‑149 (MQ=255)
     tttATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGGTCGGGTAAGAtgt                                                                                                                                            >  1:1314276/1‑149 (MQ=255)
     tttATGACTAACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGGAAGAtgt                                                                                                                                            <  2:631988/149‑1 (MQ=255)
      gtATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGAtgtg                                                                                                                                           <  1:719456/148‑1 (MQ=255)
           aCTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGaa                                                                                                                                      >  1:1895161/1‑149 (MQ=255)
            cTTACCCTTATCGCACTACAATGGCACTCAACACCTAACAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGaaa                                                                                                                                     <  1:483671/149‑1 (MQ=255)
               aCCCTTATCGCACTACAATGGCACTAAACAGCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGCGaaaaaa                                                                                                                                  >  1:1040490/1‑149 (MQ=255)
                cccTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGaaaaaaa                                                                                                                                 >  2:1016338/1‑149 (MQ=255)
                 ccTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGAGGTGGAGCAAAAAAt                                                                                                                                >  1:1863344/1‑149 (MQ=255)
                  cTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATc                                                                                                                               <  1:1200079/149‑1 (MQ=255)
                   ttATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATcc                                                                                                                              >  1:765194/1‑149 (MQ=255)
                   ctATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGCGTAAGATGTGGAGAAAAAAATcc                                                                                                                              <  1:1060979/148‑1 (MQ=255)
                     aTCGCACTACAATGGCACTCAACGCCTATCAGTACANGTAAACTACAATTATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAANGCCTAAAGGAATCAGGGAGTTCGGTTAAGATGTGGATAAAAAAATCCGc                                                                                                                            <  2:1082625/149‑1 (MQ=255)
                      tCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGcc                                                                                                                           >  1:1238256/1‑149 (MQ=255)
                          aCTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAAc                                                                                                                       >  1:1751192/1‑149 (MQ=255)
                            tACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGCATCAGGGCGTTCGGGTAAGATGTGGAGAAAAAAAACAGCCAAAccc                                                                                                                     >  1:1214106/1‑149 (MQ=255)
                            tACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGAGGTGGAGAAAAAAATCCGCCAAAcac                                                                                                                     >  1:752358/1‑147 (MQ=255)
                                aTGGCACTCAACACCTATCAGTACAGGGAACCTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACa                                                                                                                 >  1:1430023/1‑149 (MQ=255)
                                 tGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACaa                                                                                                                >  1:1255050/1‑149 (MQ=255)
                                 tGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACaa                                                                                                                >  2:295549/1‑149 (MQ=255)
                                     ccTCAACACCTATCAGTACAGTGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAgcgc                                                                                                            <  1:1717128/148‑1 (MQ=255)
                                     aCTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCACACCCTACAAgcgc                                                                                                            >  2:1667193/1‑149 (MQ=255)
                                         aaCACCTATCAGTACAGGGAAACTACAATGATTGACCCTAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGAAGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCt                                                                                                        >  2:448097/1‑149 (MQ=255)
                                         aaCACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCt                                                                                                        >  1:2000330/1‑149 (MQ=255)
                                         aaCACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCt                                                                                                        >  1:1901651/1‑149 (MQ=255)
                                             ccTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGAcg                                                                                                    <  1:597241/149‑1 (MQ=255)
                                              cTATCAGTACAGGGAAACTACAGTGATTGACCCGAAAAAAATTGAGCAAATCGCTCACCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCGACAAGCGCAGCTGAcgc                                                                                                   >  2:927025/1‑149 (MQ=255)
                                              cTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAACCCCTACCAGCGCAGCTGAcgc                                                                                                   >  2:116540/1‑149 (MQ=255)
                                               tATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCCGCTGAcgcg                                                                                                  >  2:3527/1‑149 (MQ=255)
                                                aTCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGCCgcgc                                                                                                 >  2:1309325/1‑149 (MQ=255)
                                                 tCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGcc                                                                                                >  1:147052/1‑149 (MQ=255)
                                                        atggaacatACCATGACTGACCCGAAAAAAATTGAGCACATCGCCCGTACAGTCAACGAAGCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCt                                                                                         <  1:565681/141‑1 (MQ=255)
                                                        aGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGGTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCCCGAACt                                                                                         >  1:11458/1‑149 (MQ=255)
                                                          ggAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTg                                                                                       >  1:1303210/1‑149 (MQ=255)
                                                          ggAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGAGCAGCTGACGCGCCTCGACCTTg                                                                                       >  2:858451/1‑149 (MQ=255)
                                                             aaCTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTaa                                                                                    <  1:492936/149‑1 (MQ=255)
                                                              aCTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACCAGCGCAGCTGCCGCGCCCCGGGCTTGTAAg                                                                                   <  2:221800/149‑1 (MQ=255)
                                                                tACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCACAGCGGAGGCGCCTCGACCTTGTAAGcc                                                                                 >  2:1840976/1‑149 (MQ=255)
                                                                   aaTGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGAAGTGGAGCACAAAATCCGCCAAACCCTACAAGAGCAGCTGAAGCGCCTCGCCCATGTAAGCCGTg                                                                              >  1:917541/1‑149 (MQ=255)
                                                                    aTGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGa                                                                             <  1:1605041/149‑1 (MQ=255)
                                                                     tGATTGACCCGACAAAAATTGAGCAAACCGCTCGCCAGGTTCACGAATCAACGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGCGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGATCTTGTAAGCCGTGaa                                                                            <  2:1161977/149‑1 (MQ=255)
                                                                     agATTAACCAGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGaa                                                                            <  2:42501/148‑1 (MQ=255)
                                                                       aTTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAaga                                                                          <  2:881388/149‑1 (MQ=255)
                                                                         tGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGTCTCGACCTTGTAAGCCGTGAAGAGt                                                                        <  1:1038869/149‑1 (MQ=255)
                                                                              cGAAAAAAATTGAGCAAATCGCTCGCCAGGATCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGAc                                                                   >  1:1043450/1‑149 (MQ=255)
                                                                               gAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACg                                                                  >  1:32252/1‑149 (MQ=255)
                                                                                 aaaaaaTTGAGCAAATCGCTCGCCAGGTTCACGAATCCATGCCCAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTc                                                                <  1:1868849/149‑1 (MQ=255)
                                                                                     aaTTGAGCAAATCGCTCGCCAGGTCCACGAATCAATGCCTAAAGGAAGCAGGGAGTTCGGGTAAGATGTGGAGAAATTAAGCAACTTTGATGGACTGGGTCAGCTGACGCGCCTCGACCTTGTAAGCCGTGATGAGTTCGACGTCCaaa                                                            <  2:1305254/149‑1 (MQ=255)
                                                                                        tGAGCACATCGCTCGCCCGGTTCGCGAATCAACGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGAGCCATCCTCGCGTAAGCCGTGAAGAGTTCGACGTCCAAACGc                                                         <  1:487281/149‑1 (MQ=255)
                                                                                          aGCAAATCGCTCGCCAGGTTCACGAATCAATGCATAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAGTCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGCCCTTGTAAGCCGGGAGGAGTTCGACGTCCCAACGCaa                                                       >  2:413386/1‑149 (MQ=255)
                                                                                            cAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAGATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGGGAAGAGTTCGACGGCCACACGCAAgc                                                     >  2:774571/1‑148 (MQ=255)
                                                                                             aaaTCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTc                                                    >  1:703422/1‑149 (MQ=255)
                                                                                                 cgctcgcCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGt                                                <  1:742638/149‑1 (MQ=255)
                                                                                                   ctcgcCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTa                                              >  1:228348/1‑149 (MQ=255)
                                                                                                       ccAGGTGCACGAATCAACGCCTGAAGGGATCAGGGAGTTCGAGTAAGATGGTGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGa                                          <  2:1375513/149‑1 (MQ=255)
                                                                                                          ggTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTc                                       <  2:1863286/149‑1 (MQ=255)
                                                                                                          ggTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTc                                       <  2:1358130/149‑1 (MQ=255)
                                                                                                           gatCACGGATCAGTGCCTAGAGAAGACATGCATGTCAGGTAAGATGGGGAGAAAAAAAACCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCg                                      >  2:1554601/3‑149 (MQ=255)
                                                                                                                     cAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTg                            <  1:321209/149‑1 (MQ=255)
                                                                                                                      aaTGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCCCCTCGGCCTTGGAAGCCGGGAAGGGTTCGACGTCCGATCGCAAGGCCTGTTACGGACTCGTGAAAAACTgg                           <  2:1761746/149‑1 (MQ=255)
                                                                                                                      aaTGCCTAAAGGAAGCAGGGAGTTCGGGGAAGGTGGGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTgg                           <  2:338285/149‑1 (MQ=255)
                                                                                                                       aTGCCTAGAGGAATCAGGGAGTTCGGGTAAGATGTGGAGGAAAAACTCCGCCAAACCCTACAAGCGCGGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGAAGTCCAAACCAACGTCCCGTTCTGGTCTCGTgctttacttgc                          <  2:890528/149‑11 (MQ=255)
                                                                                                                       aTGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGgc                          >  2:1738301/1‑149 (MQ=255)
                                                                                                                          ccTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGgcgat                       >  1:1309902/1‑147 (MQ=255)
                                                                                                                          ccTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGGCCAAACGCAAGTCCTGTTACGGACTCGCGAAAAACTGGCcct                       >  1:605925/1‑149 (MQ=255)
                                                                                                                            tAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAACCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCgctgc                     <  2:1565122/149‑1 (MQ=255)
                                                                                                                              aaGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCtt                   <  1:1603549/149‑1 (MQ=255)
                                                                                                                               aGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAACCCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACCCACGTCCTGTTACGGACTCGTGAAAACCTGGCGCTCCTTg                  >  2:536023/1‑149 (MQ=255)
                                                                                                                                  actCAGGGAGTTCGGGTAAGAAGAAAAAAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAAc               <  2:1076964/147‑1 (MQ=255)
                                                                                                                                   aTCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGGGCTGATTGAACa              >  2:1134292/1‑149 (MQ=255)
                                                                                                                                    tCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACaa             >  1:1588709/1‑149 (MQ=255)
                                                                                                                                    tCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACaa             >  2:1491495/1‑149 (MQ=255)
                                                                                                                                    tCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGAACTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGGGAAAAAATGGCGCTGCTTGAACac             >  2:324663/1‑148 (MQ=255)
                                                                                                                                     cAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACaaa            <  2:538937/149‑1 (MQ=255)
                                                                                                                                        ggAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGaa         >  2:809014/1‑149 (MQ=255)
                                                                                                                                               gggtAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTATAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACCCGTGAAAAACAGGCGCGGCTTGAACAAAGAAACACCgg  >  2:342040/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TCTGGTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGA  >  CP009273/3178131‑3178422

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: