Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,479,900 T→C 7.8% V12A (GTA→GCA) ‡ PP_3971 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,479,9000TC 7.8% 180.0 / 4.6 64V12A (GTA→GCA) ‡PP_3971hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (32/27);  new base C (2/3);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

TGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGATGCCGTGAAAGCGGCCTGGGCTCTCCTGATCGTGGCCCTACCGCTTGTTGGCCTGGGGATATGGGGCGTTGCTGGCCCGCGA  >  NC_002947/4479755‑4480048
                                                                                                                                                 |                                                                                                                                                    
tGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGc                                                                                                                                                   <  7:213071/149‑1 (MQ=255)
   gCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGc                                                                                                                                                <  6:289811/149‑1 (MQ=255)
   gCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGACCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCgg                                                                                                                                                 >  1:114240/1‑148 (MQ=255)
      aCAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAAt                                                                                                                                             >  7:141132/1‑149 (MQ=255)
        aaaaTGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAAttc                                                                                                                                           >  2:331713/1‑149 (MQ=255)
           aTGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAAttcttc                                                                                                                                        <  5:148193/149‑1 (MQ=255)
           aTGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAAttcttc                                                                                                                                        <  1:331713/149‑1 (MQ=255)
            tGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGCTACGGCAAtccagaa                                                                                                                                       <  5:115274/149‑7 (MQ=255)
            agTTCGATCAAAAGCAGCATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGCTACGGCAAtccagaa                                                                                                                                       <  7:114542/148‑7 (MQ=255)
                      aaaGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGAc                                                                                                                             <  8:141132/149‑1 (MQ=255)
                      aaaGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGAc                                                                                                                             >  2:137120/1‑149 (MQ=255)
                          cAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCtgt                                                                                                                         >  7:130842/1‑149 (MQ=255)
                            ggATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATtcttct                                                                                                                                       <  8:196063/133‑1 (MQ=255)
                            ggATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATtcttct                                                                                                                                       >  7:196063/1‑133 (MQ=255)
                            ggATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTgg                                                                                                                       >  5:92622/1‑149 (MQ=255)
                             gATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGAATTGCCGCTACGGCAAtccaga                                                                                                                                        >  1:137191/1‑126 (MQ=255)
                             gATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGAATTGCCGCTACGGCAAtccaga                                                                                                                                        <  2:137191/131‑6 (MQ=255)
                               tAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGc                                                                                                                     >  8:155939/1‑148 (MQ=255)
                                aaCTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAAttcttc                                                                                                                                        >  2:240970/1‑128 (MQ=255)
                                aaCTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAAttcttc                                                                                                                                        <  1:240970/128‑1 (MQ=255)
                                  cTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATc                                                                                                                 <  8:150857/149‑1 (MQ=255)
                                     ccGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGc                                                                                                                                     <  8:258757/126‑1 (MQ=255)
                                     ccGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGc                                                                                                                                     >  7:258757/1‑126 (MQ=255)
                                       gCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGt                                                                                                               >  7:220471/1‑146 (MQ=255)
                                       gCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAg                                                                                                            <  8:220471/149‑1 (MQ=255)
                                             tCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCgtgt                                                                                                       >  4:283373/1‑148 (MQ=255)
                                               aTCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGGAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTcc                                                                                                    >  3:1802/1‑149 (MQ=255)
                                               aTCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAAAGATCATGAACGGTCTCGCCAGTTACTTCTGGAATGCCGTAGCGGCAAATCTTAAGCTCATCGACCTGTGGGCCATCGACACCCTCTTcc                                                                                                    >  4:335690/1‑149 (MQ=255)
                                                tCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCg                                                                                                   <  4:332204/149‑1 (MQ=255)
                                                 ccTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTTGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCgg                                                                                                  >  7:134804/1‑149 (MQ=255)
                                                 ccTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCg                                                                                                   >  1:327993/1‑148 (MQ=255)
                                                     gTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGc                                                                                                                     >  7:327223/1‑126 (MQ=255)
                                                     gTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGc                                                                                                                     <  8:327223/126‑1 (MQ=255)
                                                      tGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGa                                                                                             <  6:92622/149‑1 (MQ=255)
                                                           ggTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGa                                                                                        <  1:137120/149‑1 (MQ=255)
                                                               ctgcACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGa                                                                                    <  7:310628/149‑1 (MQ=255)
                                                                tgcACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGAc                                                                                            <  1:154518/140‑1 (MQ=255)
                                                                tgcACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGAc                                                                                            >  2:154518/1‑140 (MQ=255)
                                                                      cGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTggg                                                                                                                      <  6:268803/108‑1 (MQ=255)
                                                                      cGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTggg                                                                                                                      >  5:268803/1‑108 (MQ=255)
                                                                      cGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGATGCCGt                                                                              >  8:120797/1‑148 (MQ=255)
                                                                        gCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGATGCCGTGaa                                                                           <  7:155939/149‑1 (MQ=255)
                                                                           ccTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGATGCCGTGAAAGc                                                                        <  3:335690/149‑1 (MQ=255)
                                                                                  cccTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCGGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGATGCCGTGAAAGCGGCCTgg                                                                 >  6:92341/1‑149 (MQ=255)
                                                                                  cccTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGATGCCGTGAAAGCGGCCTgg                                                                 >  6:127539/1‑149 (MQ=255)
                                                                                  cccTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGATGCCGTGAAAGCGGCCTgg                                                                 >  8:41275/1‑149 (MQ=255)
                                                                                   ccTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGATGCCGTGAAAGCGGCCTggg                                                                <  7:41275/149‑1 (MQ=255)
                                                                                   ccTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGATGCCGTGAAAGCGGCCTggg                                                                <  2:114240/149‑1 (MQ=255)
                                                                                      gACCACCTCTAGAGGAGGAATGATCATGATCGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGATGCCGTGAAAGCGGCCTggg                                                                <  7:159542/146‑1 (MQ=255)
                                                                                      gACCACCTCTAGAGGAGGAATGATCATGAGCGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGATGCCGTGAAAGCGGCCTggg                                                                >  8:159542/1‑146 (MQ=255)
                                                                                                  aggaggAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGATGCCGTGAAAGCGGCCTGGGCTCTCCTGATCGTgg                                                 >  1:135972/1‑149 (MQ=255)
                                                                                                    gaggaATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGATGCCGTGAAAGCGGCCTGGGCTCTCCTGATCGTGGc                                                >  7:204889/1‑148 (MQ=255)
                                                                                                         aTGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGATGCCGTGAAAGCGGCCTGGGCTCTCCTGATCGTGGCCCTAcc                                          >  8:178042/1‑149 (MQ=255)
                                                                                                             tCATGACCGGTCTCGCCAGTTACTTCTGGATGGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGATGCCGTGAAAGCGGCCTGGGCTCTCCTGATCGTGGCCCTACCGCtt                                      <  6:337884/149‑1 (MQ=255)
                                                                                                             tCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGATGCCGTGAAAGCGGCCTGGGCTCTCCTGATCGTGGCCCTACCGCtt                                      >  7:327643/1‑149 (MQ=255)
                                                                                                               aTGAACGGTCTCGCCGGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGATGCCGTGAAAGCGGCCTGGGCTCTCCTGATCGTGGCCCTACCGCttg                                     <  5:92341/148‑1 (MQ=255)
                                                                                                                tGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGATGCCGTGTAAGCGGCCTGGGCTCTCCTGATCGTGGCCCTACCGCttgtt                                   >  5:308325/1‑149 (MQ=255)
                                                                                                                   aCGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGATGCCGTGAAAGCGGCCTGGGCTCTCCTGATCGTGGCCCTACCGCTTGTTGGc                                <  4:1802/149‑1 (MQ=255)
                                                                                                                           gCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGATGCCGTGAAAGCGGCCTGGGCTCTCCTGATCGTGGCCCTACCGCTTGTTGGCCTGGGGa                         >  2:272949/1‑148 (MQ=255)
                                                                                                                                    ttCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGATGCCGTGAAAGCGGCCTGGGCTCTCCTGATCGTGGCCCTACCGCTTGTTGGCCTGGGGATATGGGGCGt               >  2:243211/1‑149 (MQ=255)
                                                                                                                                    agaagaaTTGCCGCTACGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTAAGCGTGTTCAGGAGCGACAAGACCGATGCCGTGAAAGCGGCCTGGGCTCTCCTGATCGTGGCCCTACCGCTTGTTGGcc                               >  6:52911/7‑133 (MQ=255)
                                                                                                                                         gATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGATGCCGTGAAAGCGGCCTGGGCTCTCCTGATCGTGGCCCTACCGCTTGTTGGCCTGGGGATATGGGGCGTTGCt           <  2:327993/148‑1 (MQ=255)
                                                                                                                                            tGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGATGCCGTGAAAGCGGCCTGGGCTCTCCTGATCGTGGCCCTACCGCTTGTTGGCCTGGGGATATGGGGCGTTGCTGGc        <  3:283373/148‑1 (MQ=255)
                                                                                                                                                gTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGATGCCGTGAAAGCGGCCTGGGCTCTCCTGATCGTGGCCCTACCGCTTGTTGGCCTGGGGATATGGGGCGTTGCTGGCCcgcg   >  8:47117/1‑149 (MQ=255)
                                                                                                                                                 tAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGATGCCGTGAAAGCGGCCTGGGCTCTCCTGATCGTGGCCCTACCGCTTGTTGGCCTGGGGATATGGGGCGTTGCTGGCCCGCGa  <  7:120797/149‑1 (MQ=255)
                                                                                                                                                 tAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGATGCCGTGAAAGCGGCCTGGGCTCTCCTGATCGTGGCCCTACCGCTTGTTGGCCTGGGGATATGGGGCGTTGCTGGCCCGCGa  <  1:243211/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGATGCCGTGAAAGCGGCCTGGGCTCTCCTGATCGTGGCCCTACCGCTTGTTGGCCTGGGGATATGGGGCGTTGCTGGCCCGCGA  >  NC_002947/4479755‑4480048

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑