Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,789,896 C→G 5.7% G1315A (GGA→GCA)  pvdI ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,789,8960CG 5.7% 230.2 / 2.8 70G1315A (GGA→GCA) pvdInon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base C (34/32);  new base G (2/2);  total (36/34)
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.29e-01

GCCGATGTTGATCTTCTCGTTGCCGCTGTAGCGCGCCAGCAGC‑AGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTGGCCTTCGCTCAGGCGCTCGCGTTGCCAGCTGGCATAG  >  NC_002947/4789748‑4790042
                                                                                                                                                     |                                                                                                                                                  
gCCGATGTTGATCTTCTCGTTGCCGCTGTAGCGCGCCAGCAGC‑AGGCCGAACGAGGC‑AGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTcc                                                                                                                                                   >  5:109195/1‑149 (MQ=255)
 ccGATGTTGATCTTCTCGTTGCCGCTGTAGCGCGCCAGCAGC‑AGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTcc                                                                                                                                                   <  3:105263/149‑1 (MQ=255)
 ccGATGTTGATCTTCTCGTTGCCGCTGTAGCGCGCCAGCAGC‑AGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTcc                                                                                                                                                   >  3:308289/1‑149 (MQ=255)
  cGATGTTGATCTTCTCGTTGCCGCTGTAGCGCGCCAGCAGC‑AGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCg                                                                                                                                                  >  7:137479/1‑149 (MQ=255)
     tGTTGATCTTCTCGTTGCCGCTGTAGCGCGCCAGCAGC‑AGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGccc                                                                                                                                               <  8:212278/149‑1 (MQ=255)
      gTTGATCTTCTCGTTGCCGCTGTAGCGCGCCAGCAGCTAGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGccc                                                                                                                                               >  1:68406/1‑149 (MQ=255)
         gATCTTCTCGTTGCCGCTGTAGCGCGCCAGCAGC‑AGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGa                                                                                                                                            >  2:35624/1‑148 (MQ=255)
          aTCTTCTCGTTGCCGCTGTAGCGCGCCAGCAGC‑AGGCCGAACGAGGCTAGGATGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCCCGTGCCGACGTCCGCCCCGa                                                                                                                                            >  8:254759/1‑147 (MQ=255)
          aTCTTCTCGTTGCCGCTGTAGCGCGCCAGCAGC‑AGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGAt                                                                                                                                           <  2:135647/148‑1 (MQ=255)
              tctcGTTGCCGCTGTAGCGCGCCAGCAGC‑AGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCtt                                                                                                                                      >  4:34211/1‑149 (MQ=255)
               ctcGTTGCCGCTGTAGCGCGCCAGCAGC‑AGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCt                                                                                                                                       <  8:196529/147‑1 (MQ=255)
               ctcGTTGCCGCTGTAGCGCGCCAGCAGC‑AGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCt                                                                                                                                       >  7:196529/1‑147 (MQ=255)
                  gTTGCCGCTGTAGCGCGCCAGCAGC‑AGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGc                                                                                                                                   <  5:48695/148‑1 (MQ=255)
                   ctgctgcTGGCGCGCTACAGCGGC‑AGGCCGAACGGGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAg                                                                                                                                 >  2:54670/2‑149 (MQ=255)
                            tAGCGCGCCAGCAGC‑AGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGc                                                                                                                          <  4:308289/147‑1 (MQ=255)
                            tAGCGCGCCAGCAGC‑AGGCCGAACGAGGCTAGGAAGACGGGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCg                                                                                                                        <  4:262065/149‑1 (MQ=255)
                              gcgcgcCAGCAGC‑AGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGTGCCGt                                                                                                                       <  8:147232/148‑1 (MQ=255)
                              gcgcgcCAGCAGC‑AGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGt                                                                                                                       <  8:137479/148‑1 (MQ=255)
                                gcgcCAGCAGC‑AGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGGACGTCTGGACGTGCGcct                                                                                                                                               >  2:69213/1‑121 (MQ=255)
                                gcgcCAGCAGC‑AGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGGACGTCTGGACGTGCGcct                                                                                                                                               <  1:69213/122‑2 (MQ=255)
                                 cgcCAGCAGC‑AGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGc                                                                                                                   >  1:282021/1‑149 (MQ=255)
                                   ccAGCAGC‑AGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCgg                                                                                                                 >  2:144768/1‑149 (MQ=255)
                                         gc‑agGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCa                                                                                                           <  7:65129/149‑1 (MQ=255)
                                              gCCGAATGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGt                                                                                                        >  4:263841/1‑148 (MQ=255)
                                                   aCGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTgc                                                                                                  <  6:62568/149‑1 (MQ=255)
                                                     gAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTgcag                                                                                                >  4:266577/1‑149 (MQ=255)
                                                     gAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTgcag                                                                                                <  1:35624/149‑1 (MQ=255)
                                                        gCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCAc                                                                                             >  3:77877/1‑149 (MQ=255)
                                                        gCTAGGAAGACGTGGAACAGCGTGGCGTTACAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCAc                                                                                             >  1:28826/1‑149 (MQ=255)
                                                            ggAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGc                                                                                          <  8:15609/148‑1 (MQ=255)
                                                               aGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGaaa                                                                                      <  1:121273/149‑1 (MQ=255)
                                                               aGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGaaa                                                                                      >  1:296111/1‑149 (MQ=255)
                                                               aGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGaaa                                                                                      <  2:282021/149‑1 (MQ=255)
                                                                 aCGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAat                                                                                     >  3:183071/1‑148 (MQ=255)
                                                                    tGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAatcatc                                                                                 <  3:263841/149‑1 (MQ=255)
                                                                      gAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCtt                                                                               <  6:75695/149‑1 (MQ=255)
                                                                       aaCAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTc                                                                              <  6:109195/149‑1 (MQ=255)
                                                                          aGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCTGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAg                                                                           <  2:170760/149‑1 (MQ=255)
                                                                           gCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAgg                                                                          <  3:34211/149‑1 (MQ=255)
                                                                            cGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGt                                                                         >  5:294516/1‑149 (MQ=255)
                                                                            cGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGt                                                                         >  4:324172/1‑149 (MQ=255)
                                                                             gTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTg                                                                        >  5:146728/1‑149 (MQ=255)
                                                                              tGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTgc                                                                       >  2:90969/1‑149 (MQ=255)
                                                                                     gCAGCGCACTGCCAGGCTCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGc                                                                 <  1:90969/148‑1 (MQ=255)
                                                                                     gCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGc                                                                 >  6:43467/1‑148 (MQ=255)
                                                                                      cAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGcc                                                                <  4:183071/148‑1 (MQ=255)
                                                                                      cAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGc                                                                 >  3:270719/1‑147 (MQ=255)
                                                                                              tGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGcc                                                         >  2:108819/1‑147 (MQ=255)
                                                                                                ccAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGc                                                     >  6:313960/1‑149 (MQ=255)
                                                                                                ccAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGc                                                     >  3:97055/1‑149 (MQ=255)
                                                                                                 cAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCt                                                    <  4:77877/149‑1 (MQ=255)
                                                                                                 cAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCt                                                    <  2:296111/149‑1 (MQ=255)
                                                                                                    gcgcCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCGGCCAGTAGCCGAGCTggg                                                 <  2:28826/149‑1 (MQ=255)
                                                                                                      gctgcAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCt                                               <  5:24402/146‑1 (MQ=255)
                                                                                                      gccgcAGGTCCGCCGTCACCTGCGCGGGGAGGCGCACGTCCAGACGTCCGCGCCTATAGCTTTGCAGCTCTGGCCGTGGGCGGTCGGc                                                                                                            <  5:197527/88‑1 (MQ=255)
                                                                                                      gacgcAGGTCAGACGTCAACTGCGCGGGAAGGCGAAAGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGc                                                                                                            >  6:197527/3‑88 (MQ=255)
                                                                                                       ccgcAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCtt                                              <  1:144768/149‑1 (MQ=255)
                                                                                                                cgccgTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTgg                                      <  2:68406/148‑1 (MQ=255)
                                                                                                                       aaCTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTGGCCTTCGCt                              >  4:345794/1‑149 (MQ=255)
                                                                                                                          tGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTGGCCTTCGCTCAg                           >  7:165450/1‑149 (MQ=255)
                                                                                                                               gggCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTGGCCTTCGCTCAGGCGCt                      >  3:100573/1‑149 (MQ=255)
                                                                                                                                   gggCGGACGTCTGGACGTGCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGc                                            >  7:328883/2‑123 (MQ=255)
                                                                                                                                   gggCGGACGTCTGGACGTGCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGc                                            <  8:328883/122‑1 (MQ=255)
                                                                                                                                   aGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTGGCCTTCGCTCAGGCGCTcgcg                  >  1:236175/1‑149 (MQ=255)
                                                                                                                                   aGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGAAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTGGCCTTCGCTCAGGCGCTcgcg                  >  1:109475/1‑149 (MQ=255)
                                                                                                                                     gcgcACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATc                                         >  5:78837/1‑124 (MQ=255)
                                                                                                                                      cgcACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTGGCCTTCGCTCAGGCGCTCGCGTTg               <  5:43467/149‑1 (MQ=255)
                                                                                                                                       gcACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTGGCCTTCGCTCAGGCGCTCGCGTTGc              >  6:283910/1‑149 (MQ=255)
                                                                                                                                       gcACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTGGCCTTCGCTCAGGCGCTCGCGTTGc              <  1:54670/149‑1 (MQ=255)
                                                                                                                                          cGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTGGCCTTCGCTCAGGCGCTCGCGTTGCCAg           <  4:97055/149‑1 (MQ=255)
                                                                                                                                           gTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTGGCCTTCGCTCAGGCGCTCGCGTTGCCAGc          >  2:87784/1‑149 (MQ=255)
                                                                                                                                                  cGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTGGCCTTCGCTCAGGCGCTCGCGTTGCCAGCTGGCATa   <  2:236175/149‑1 (MQ=255)
                                                                                                                                                   gTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTGGCCTTCGCTCAGGCGCTCGCGTTGCCAGCTGGCATAg  >  4:186107/1‑149 (MQ=255)
                                                                                                                                                     ccGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAg                                                           <  3:168083/90‑1 (MQ=255)
                                                                                                                                                     ccGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAg                                                           >  4:168083/1‑90 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GCCGATGTTGATCTTCTCGTTGCCGCTGTAGCGCGCCAGCAGC‑AGGCCGAACGAGGCTAGGAAGACGTGGAACAGCGTGGCGTTGCAGCGCACTGCCAGGCGCCGCAGGTCCGCCGTCAACTGCGCGGGCAGGCGCACGTCCAGACGTCCGCCCCGATAGCTTTGCAGCTCGGGCCGTGGGCGGTCGGCAGGTAACTGCAGCACTGCGAAATCATCTTCCAGGTGCGCCTGCCAGTAGCCGAGCTGGGCTTGCATCTGGCCTTCGCTCAGGCGCTCGCGTTGCCAGCTGGCATAG  >  NC_002947/4789748‑4790042

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