Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 602,109 T→C 6.3% L94P (CTG→CCG)  ribE → riboflavin synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947602,1090TC 6.3% 236.7 / 4.7 79L94P (CTG→CCG) ribEriboflavin synthase
Reads supporting (aligned to +/- strand):  ref base T (37/37);  new base C (2/3);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.70e-01

AACGGTGTATGCCTCACGGCCGTAGAACTGCCGGGCGACGGCTTCTGGGCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCCGCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGCCCACAAGGGTTCGATCACCGTCGACG  >  NC_002947/601961‑602255
                                                                                                                                                    |                                                                                                                                                  
aaCGGTGTATGCCTCACGGCCGTAGAACTGCCGGGCGACGGCTTCTGGGCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCt                                                                                                                                                    >  4:305210/1‑149 (MQ=255)
aaCGGTGTATGCCTCACGGCCGTAGAACTGACGGGCGACGGCTTCTGGGCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCt                                                                                                                                                    >  1:219580/1‑149 (MQ=255)
   ggTGTATGCCTCACGGCCGTAGAACTGCCGGGCGACGGCTTCTGGGCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTggg                                                                                                                                                 <  2:229852/149‑1 (MQ=255)
     tgtATGCCTCACGGCCGTAGAACTGCCGGGCGACGGCTTCTGGGCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTgggcg                                                                                                                                               >  7:192588/1‑149 (MQ=255)
         tGCCTCACGGCCGTAGAACTGCCGGGCGACGGCTTCTGGGCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCCaggcgggta                                                                                                                                           <  2:118714/149‑3 (MQ=255)
         tGCCGCACGGCCGTAGAACTGCCGGGCGACGGCGTCTGGGCCGATGTCAGCGTCGAAACCCTCAAGCCCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCa                                                                                                                                           <  4:209348/149‑1 (MQ=255)
          gCCTCACGGCCGTAGAACTGCCGGGCGACGGCTTCTGGGCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCAc                                                                                                                                          >  5:113760/1‑149 (MQ=255)
           ccTCACGGCCGTAGAACTGCCGGGCGACGGCTTCTGGGCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCAcc                                                                                                                                         <  1:72394/149‑1 (MQ=255)
           ccTCACGGCCGTAGAACTGCCGGGCGACGGCTTCTGGGCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCAc                                                                                                                                          <  8:333544/148‑1 (MQ=255)
            cTCACGGCCGTAGAACTGCCGGGCGACGGCTTCTGGGCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCt                                                                                                                                        >  8:75968/1‑149 (MQ=255)
            cTCACGGCCGTAGAACTGCCGGGCGACGGCTTCTGGGCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCt                                                                                                                                        <  8:192588/149‑1 (MQ=255)
               aCGGCCGTAGAACTGCCGGGCGACGGCTTCTGGGCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCCaggcgggtagtgg                                                                                                                                       <  8:320513/147‑7 (MQ=255)
                  gCCGTAGAACTGCCGGGCGACGGCTTCTGGGCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAg                                                                                                                                  <  6:113760/149‑1 (MQ=255)
                    cGTAGAACTGCCGGGCGACGGCTTCTGGGCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCACGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCg                                                                                                                                >  8:82106/1‑149 (MQ=255)
                    cGTAGAACTGCCGGGCGACGGCTTCTGGGCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCg                                                                                                                                >  2:212310/1‑149 (MQ=255)
                         aaCTGCCGGGCGACGGCTTCTGGGCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCAc                                                                                                                           <  6:295629/149‑1 (MQ=255)
                         aaCTGCCGGGCGACGGCTTCTGGGCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCAc                                                                                                                           <  2:219580/149‑1 (MQ=255)
                           cTGCCGGGCGACGGCTTCTGGGCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGt                                                                                                                         <  2:176556/149‑1 (MQ=255)
                               cGGGCTACTGCTTCTGGGCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTTGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGAc                                                                                                                     <  1:135943/149‑1 (MQ=255)
                               cGGGCGACGGCTTCTGGGCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGAc                                                                                                                     >  2:138901/1‑149 (MQ=255)
                                     aCGGCTTCTGGGCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCTAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGtt                                                                                                               <  1:116883/149‑1 (MQ=255)
                                     aCGGCTTCTGGGCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGtt                                                                                                               >  1:36834/1‑149 (MQ=255)
                                               ggCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTCTGCGGGCACCTGTTCAGCGGCCTCGTCGACTGGGTTGGCGaaa                                                                                                        <  2:98594/146‑1 (MQ=255)
                                                  cGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCAttt                                                                                                  >  6:12950/1‑149 (MQ=255)
                                                  cGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCAttt                                                                                                  <  6:359974/149‑1 (MQ=255)
                                                       tCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTcgcg                                                                                              <  5:180356/148‑1 (MQ=255)
                                                       tCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTcgcg                                                                                              >  5:366507/1‑148 (MQ=255)
                                                             tCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGAt                                                                                       <  6:366507/149‑1 (MQ=255)
                                                             tCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGAt                                                                                       <  5:43866/149‑1 (MQ=255)
                                                                  aCCCTCAAGCGCACCGCCTTCATCGACCGCAAGAGCGGCAGCAAGGTCAACCTGGCAAACGCCCTGACCCCCACTACCCGCCTGGGCGGGCCCCTGGTCAGCGGCCACGTCGACGGGGGTGGCGAAATCATTTCGCGCAGCGCTACCGc                                                                                  >  6:172476/1‑149 (MQ=255)
                                                                     cTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCcg                                                                               >  5:362526/1‑149 (MQ=255)
                                                                       cAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCcgcg                                                                             >  1:146988/1‑149 (MQ=255)
                                                                        aaGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCgcgc                                                                            >  6:106388/1‑149 (MQ=255)
                                                                        aaGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCgcgc                                                                            >  2:200852/1‑149 (MQ=255)
                                                                          gcgcACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCa                                                                          <  7:75968/149‑1 (MQ=255)
                                                                            gcACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGtt                                                                                                               <  2:76144/110‑1 (MQ=255)
                                                                            gcACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATc                                                                        <  2:36834/149‑1 (MQ=255)
                                                                              aCCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATc                                                                        >  2:45184/1‑147 (MQ=255)
                                                                              aCCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATc                                                                        >  7:113567/1‑147 (MQ=255)
                                                                                  ccTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGtac                                                                  >  1:132763/1‑147 (MQ=255)
                                                                                         cGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCCgcgtgc                                                           <  3:268795/149‑1 (MQ=255)
                                                                                          gACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCCgcgtgc                                                           <  1:212310/148‑1 (MQ=255)
                                                                                              tCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCCGCGTGCGTGCa                                                      >  7:324797/1‑149 (MQ=255)
                                                                                                aaGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAg                                                                                           >  4:182761/1‑110 (MQ=255)
                                                                                                aaGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAg                                                                                           <  3:182761/110‑1 (MQ=255)
                                                                                                aaGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCCGCGTGCGTGCa                                                      <  1:45184/147‑1 (MQ=255)
                                                                                                aaGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCCGCGTGCGTGCAc                                                     <  1:200852/148‑1 (MQ=255)
                                                                                                  gagCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCCGCGTGCGTGCACCg                                                   >  2:12864/1‑148 (MQ=255)
                                                                                                  gagCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCCGCGTGCGTGCACCg                                                   >  1:114512/1‑148 (MQ=255)
                                                                                                  gagCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCCGCGTGCGTGCACCg                                                   >  7:215908/1‑148 (MQ=255)
                                                                                                   agCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCCGCGTGCGTGCACCGaa                                                 <  5:36799/149‑1 (MQ=255)
                                                                                                   agCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCCGCGTGCGTGCACCGaa                                                 <  3:297885/149‑1 (MQ=255)
                                                                                                       gcagcaAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTCGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGTCCGCGCCATCCAGTTCCGCGTGCGTGCACCGAAGGaa                                             <  1:251138/149‑1 (MQ=255)
                                                                                                        cagcaAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCCGCGTGCGTGCACCGAAGGAAc                                            <  1:138901/149‑1 (MQ=255)
                                                                                                         agcaAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCCGCGTGCGTGCACCGAAGGAACt                                           >  8:35461/1‑149 (MQ=255)
                                                                                                          gcaAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCCGCGTGCGTGCACCGAAGGAACTg                                          <  7:82106/149‑1 (MQ=255)
                                                                                                              ggTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCAAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCCGCGTGCGTGCACCGAAGGAACTGGcc                                       >  5:227453/1‑148 (MQ=255)
                                                                                                                  aaCCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCCGCGTGCGTGCACCGAAGGAACTGGCCAAGTa                                  >  7:320686/1‑149 (MQ=255)
                                                                                                                     cTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCCGCGTGCGTGCACCGAAGGAACTGGCCAAGTACAt                               >  8:187595/1‑149 (MQ=255)
                                                                                                                        gAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCCGCGTGCGTGCACCGAAGGAACTGGCCAAGTACATc                              <  2:132763/147‑1 (MQ=255)
                                                                                                                        gAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCCGCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGc                            <  5:12950/149‑1 (MQ=255)
                                                                                                                         aaaaaGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCCGCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGcc                           <  7:35461/149‑1 (MQ=255)
                                                                                                                              gCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCCGCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGCCcac                        >  8:264554/1‑147 (MQ=255)
                                                                                                                               cccTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCCGCGTGCGTGCACCGAAGGAACTGGCCAAGt                                   >  7:174841/1‑135 (MQ=255)
                                                                                                                               cccTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCCGCGTGCGTGCACCGAAGGAACTGGCCAAGt                                   <  8:174841/135‑1 (MQ=255)
                                                                                                                                     cccccACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCCGCGTGCGTGCAc                                                     <  3:32014/111‑1 (MQ=255)
                                                                                                                                     cccccACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCCGCGTGCGTGCAc                                                     >  4:32014/1‑111 (MQ=255)
                                                                                                                                      ccccACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCCGCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGCCCACAAGGGTTCGa              >  2:33620/1‑149 (MQ=255)
                                                                                                                                      ccccACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCCGCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGCCCACAAGGGTTCGa              <  8:324797/149‑1 (MQ=255)
                                                                                                                                       cccACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCCGCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGCCCACAAGGGTTCGAt             <  2:114512/149‑1 (MQ=255)
                                                                                                                                        ccACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCCGCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGCCCACAAGGGTTCGATc            <  8:215908/149‑1 (MQ=255)
                                                                                                                                        ccACTACCCGCCTGGGCGGGCACCTGGTCAGCGACCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCCGCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGCCCACAAGGGTTCGAt             >  8:246901/1‑148 (MQ=255)
                                                                                                                                          ggtgcccGCCCAGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTc                                                                  >  7:139805/5‑93 (MQ=255)
                                                                                                                                          ggtgcccGCCCAGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTc                                                                  <  8:139805/89‑1 (MQ=255)
                                                                                                                                          ggtgcccGCCCAGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGAGATAACGCCCGCGCCATCCAGTTCCGCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGCCAACAAGCGTTAGATCAc          >  6:245098/5‑149 (MQ=255)
                                                                                                                                          aCTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCCGCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGCCCACAAGGGTTCGATCAc          >  6:133880/1‑149 (MQ=255)
                                                                                                                                            tACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCCGCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGCCCACAAGGGTTCGATCACCg        >  6:137295/1‑149 (MQ=255)
                                                                                                                                               ccGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCCGCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGCCCACAAGGGTTCGATCACCGTc      >  3:86348/1‑148 (MQ=255)
                                                                                                                                               ccGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTCGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCCGCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGCCCACAAGGGTTCGATCACCGTc      >  1:305780/1‑148 (MQ=255)
                                                                                                                                                 gCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCCGCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGCCCACAAGGGTTCGATCACCGTcagg   >  5:207795/1‑146 (MQ=255)
                                                                                                                                                  ccTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCCGCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGCCCACAAGGGTTCGATCACCGTCGACg  >  5:145243/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AACGGTGTATGCCTCACGGCCGTAGAACTGCCGGGCGACGGCTTCTGGGCCGATGTCAGCGTCGAAACCCTCAAGCGCACCGCCTTCATCGACCTCAAGAGCGGCAGCAAGGTCAACCTGGAAAAAGCCCTGACCCCCACTACCCGCCTGGGCGGGCACCTGGTCAGCGGCCACGTCGACGGGGTTGGCGAAATCATTTCGCGCAGCGATAACGCCCGCGCCATCCAGTTCCGCGTGCGTGCACCGAAGGAACTGGCCAAGTACATCGCCCACAAGGGTTCGATCACCGTCGACG  >  NC_002947/601961‑602255

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: ‑