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,2770GT96.7% 162.3 / ‑5.9 64E27* (GAA→TAA) yqiCBMFP family putative fusogenic protein
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (33/26);  total (36/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

TGGTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGC  >  CP009273/3178133‑3178420
                                                                                                                                                |                                                                                                                                               
tGGTTTATGACTTACCCTTATCGCGCTACAATGGCACTCAACTCCTGTCGCTCTACGTTAACTACATTGATTGTCCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGa                                                                                                                                             <  2:28874/149‑1 (MQ=255)
tGGTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACCCCTGTCAGTACAGGGAAACTACAATGATTGACCCGAAAAAACTTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGa                                                                                                                                             >  2:403345/1‑149 (MQ=255)
tGGTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGTGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGa                                                                                                                                             >  1:672910/1‑149 (MQ=255)
tGGTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGCAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTCAAGGAATCAGGGAGTTCGAGTAAGa                                                                                                                                             >  2:403426/1‑149 (MQ=255)
    ttATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGAtgtg                                                                                                                                         <  1:693781/149‑1 (MQ=255)
    ttATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAAGGGAGCAAATCGCTCGCCAGGTTCACGAAGCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGAtgtg                                                                                                                                         <  2:143554/149‑1 (MQ=255)
     tATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGTCCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTgg                                                                                                                                        <  2:14180/149‑1 (MQ=255)
        gACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCTGGGAGTTCGGGTAAGATGTGgaga                                                                                                                                     <  2:768850/149‑1 (MQ=255)
          cTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGCCCCGCAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCACTGCCTAACGGAATCAGGGAGTTCGGGTAAGATGTGGAGaag                                                                                                                                   >  1:422857/1‑148 (MQ=255)
             aCCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCAAAAGGAATCAGGGAGATCGGGAAAGATGTGGAGaaaaaa                                                                                                                                >  2:226815/1‑149 (MQ=255)
             aCCCTTACCGCACAACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGaaaaaa                                                                                                                                <  2:827648/149‑1 (MQ=255)
               ccTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAAt                                                                                                                              >  2:228984/1‑149 (MQ=255)
                cTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGGGTTCGGGTAAGATGTGGAGAAAAAAATc                                                                                                                             >  2:13053/1‑149 (MQ=255)
                cTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTAAACAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATc                                                                                                                             >  2:300515/1‑149 (MQ=255)
                    tCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGTCCCGATTTAAATTGTGCTATTCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGcc                                                                                                                         <  2:63459/149‑1 (MQ=255)
                     cacACTACAATGTCCCTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGTATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCa                                                                                                                        <  1:77657/147‑1 (MQ=255)
                      gCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAAACAATGCCTAAAGGAATCAGGGAGTTCGGGAAAAATGGGGAAAAAAAAATCCGCCaa                                                                                                                       >  2:134500/1‑149 (MQ=255)
                          tACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGAAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAAccc                                                                                                                   >  2:128394/1‑149 (MQ=255)
                          tACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAAATGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAAccc                                                                                                                   >  2:806610/1‑149 (MQ=255)
                             aaTGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCGGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGACTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCCCCAAACCCTAc                                                                                                                >  1:681954/1‑149 (MQ=255)
                                    cTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCa                                                                                                         <  2:150363/149‑1 (MQ=255)
                                    cTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCCACAAGCGCa                                                                                                         >  2:206001/1‑149 (MQ=255)
                                         cacCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCACACGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGa                                                                                                    >  2:919389/1‑149 (MQ=255)
                                          acCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCCCAAAGCGCACCTGAc                                                                                                   >  1:21553/1‑149 (MQ=255)
                                          acCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGATCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCTGCCAAACCATACAAGCGCAGCtatg                                                                                                   >  1:271365/1‑146 (MQ=255)
                                           ccTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGAcg                                                                                                  <  1:936135/149‑1 (MQ=255)
                                             tATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGATCAGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGAcgcg                                                                                                >  2:885740/1‑149 (MQ=255)
                                              aTCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACgcgc                                                                                               >  2:123727/1‑149 (MQ=255)
                                                        ggAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTg                                                                                     <  2:102989/149‑1 (MQ=255)
                                                        ggAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGATCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTg                                                                                     >  2:348931/1‑149 (MQ=255)
                                                           aaCTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCCCGCCAGCTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTaa                                                                                  <  1:300515/149‑1 (MQ=255)
                                                                  aTGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGATGCGCCTCGACCTTGTAAGCCGTGa                                                                           <  2:999379/149‑1 (MQ=255)
                                                                   tGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCTACCTTCTAAGCCGTGaa                                                                          >  2:507120/1‑149 (MQ=255)
                                                                       tGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAgcg                                                                                                           >  2:744198/1‑112 (MQ=255)
                                                                       tGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAgcg                                                                                                           <  1:744198/112‑1 (MQ=255)
                                                                          cccGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCg                                                                   >  2:512928/1‑149 (MQ=255)
                                                                           ccGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGATCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAACCCTGAAAAGTTCGa                                                                  >  2:171398/1‑149 (MQ=255)
                                                                            cGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGAc                                                                 >  2:314753/1‑149 (MQ=255)
                                                                                aaaaaTTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCACACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTcc                                                             >  2:93008/1‑149 (MQ=255)
                                                                                   aaTTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCaaa                                                          <  2:952039/149‑1 (MQ=255)
                                                                                    aTTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAAc                                                         >  2:941657/1‑149 (MQ=255)
                                                                                        aGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCaa                                                     <  2:875214/149‑1 (MQ=255)
                                                                                          cAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGt                                                   >  2:976813/1‑149 (MQ=255)
                                                                                            aaTCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTcc                                                 >  1:659534/1‑149 (MQ=255)
                                                                                               cgctcgcCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGt                                              <  1:209820/149‑1 (MQ=255)
                                                                                                gctcgcCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGATCGGGGGGAAGATGAGGAGAAAAAAATCCGCCAAACCCTACAAGCGCACCTGACGCGGCTCGACCTTGTAAGCCGTGAAGAGTTCAAAGTCCAAACGCAAGTCCCGta                                             >  2:247523/1‑148 (MQ=255)
                                                                                                   cgcCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGCGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCGGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGACCTCGACGTCCAAACGCAAGTCCTGTTACg                                          <  2:672910/149‑1 (MQ=255)
                                                                                                     ccAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCCCCTCGACCATGTAAGCCGTGAAGAGTTCAACGTCCAAAAGCAAGTCCTATTACGGa                                        >  2:44786/1‑149 (MQ=255)
                                                                                                           tCACGAATCAATACCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCACACCCTACGAGCGCAGCTGACGCGCCTCGACCTTATAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGaatgctg                                  >  2:743832/1‑143 (MQ=255)
                                                                                                           tCACGAATCAAGGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGATTACGACTTCCAAACGCAAGTCCTGTTACGGACTCGTg                                  <  1:478082/149‑1 (MQ=255)
                                                                                                              cGAATCAATGCGTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGaaa                               <  2:1032549/149‑1 (MQ=255)
                                                                                                               gAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGaaaa                              >  2:449116/1‑149 (MQ=255)
                                                                                                               gAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCAACCTTGTAAGCCGAGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGaaaa                              >  2:540703/1‑149 (MQ=255)
                                                                                                                  tCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGGGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACt                           <  1:690153/149‑1 (MQ=255)
                                                                                                                   cAATGCCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAACCCGCCAAACCCTACAAGCGCAGCTGACTCGCCTCTACCTTGTAAGTCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTg                          <  1:171398/149‑1 (MQ=255)
                                                                                                                     aTGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGgc                        >  2:510902/1‑149 (MQ=255)
                                                                                                                          tAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCgctgc                   >  2:1033900/1‑149 (MQ=255)
                                                                                                                                   cAGGGAGTTCGTGTAAGATGTGGAGAAAAATATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACaaa          <  2:412758/149‑1 (MQ=255)
                                                                                                                                    atggAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTTGTCGTCCAAACGCAAGTCCTGTAAAGAAAAAATGAAAAACTGGCGCTGCTTGAACAAAg         <  1:885740/147‑1 (MQ=255)
                                                                                                                                     gggAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGa        <  1:81971/149‑1 (MQ=255)
                                                                                                                                       gAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGGTTGAACAAAGAAt      <  1:123727/149‑1 (MQ=255)
                                                                                                                                          ttCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAg   <  1:980195/149‑1 (MQ=255)
                                                                                                                                          ttCGGGCAAGATGTGGAAAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAATTGGCGCTGCTTGAACAAAGAATCAg   <  2:514824/149‑1 (MQ=255)
                                                                                                                                           tCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAAAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTAAAAAACTGGCGCTGCTTGAACAAAAAATCAGc  >  2:122213/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                               
TGGTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGC  >  CP009273/3178133‑3178420

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