Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 892,897 A→G 6.7% V273A (GTC→GCC)  pta ← phosphate acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947892,8970AG 6.7% 213.9 / 4.6 74V273A (GTC→GCC) ptaphosphate acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (36/33);  new base G (2/3);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.54e-01

GGTGTCGTACGAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCGCCGGCATTGAGCACCTGG  >  NC_002947/892749‑893043
                                                                                                                                                    |                                                                                                                                                  
ggTGTCGTACGAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGa                                                                                                                                                    >  4:119411/1‑149 (MQ=255)
       tACGAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAg                                                                                                                                             >  4:177563/1‑149 (MQ=255)
         cGAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGATGCCATTGATCGCTGCCAGGCTGACAGCGAGga                                                                                                                                           <  1:36748/149‑1 (MQ=255)
         cGAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGga                                                                                                                                           >  7:255187/1‑149 (MQ=255)
          gAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGgat                                                                                                                                          >  2:314933/1‑149 (MQ=255)
            gCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGTCAGCCTGGCAGCGatcaatg                                                                                                                                        <  8:296538/149‑7 (MQ=255)
            gCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGgatga                                                                                                                                        >  2:231515/1‑149 (MQ=255)
               ggTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATgt                                                                                                                                     <  8:255187/149‑1 (MQ=255)
                      aCCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTcgcg                                                                                                                              >  6:307146/1‑149 (MQ=255)
                        cGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGTCAGCCTGGCAGCGat                                                                                                                                             >  5:143169/1‑131 (MQ=255)
                        cGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGTCAGCCTGGCAGCGat                                                                                                                                             <  6:143169/132‑2 (MQ=255)
                         gAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGa                                                                                                                                              <  8:28042/130‑1 (MQ=255)
                          aaaGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCg                                                                                                                          >  5:34123/1‑149 (MQ=255)
                           aaGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGc                                                                                                                         >  3:105946/1‑149 (MQ=255)
                             ggATCGGCAGGCCGGCCTGCATGGCACCGCGGCACAGGCCAAGGATCCGCACGTCAGGCTTGCTGTAGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGccc                                                                                                                       >  5:105355/1‑149 (MQ=255)
                             ggATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGccc                                                                                                                       >  1:48335/1‑149 (MQ=255)
                               aTCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCgg                                                                                                                     >  5:258316/1‑149 (MQ=255)
                                tCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCggg                                                                                                                    <  1:231515/149‑1 (MQ=255)
                                     aGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTac                                                                                                               <  3:177563/149‑1 (MQ=255)
                                      ggccggccTTCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGTCTCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTacc                                                                                                              <  7:205043/149‑1 (MQ=255)
                                      ggccggccTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTacc                                                                                                              >  2:101400/1‑149 (MQ=255)
                                            ccTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAgg                                                                                                        <  1:101400/149‑1 (MQ=255)
                                            ccTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAgg                                                                                                        <  6:34123/149‑1 (MQ=255)
                                                      aCCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAgg                                                                                              >  6:174683/1‑149 (MQ=255)
                                                       ccGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCggg                                                                                                                    >  5:43399/1‑126 (MQ=255)
                                                       ccGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCggg                                                                                                                    <  6:43399/126‑1 (MQ=255)
                                                         gcgGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCa                                                                                           <  3:273094/149‑1 (MQ=255)
                                                            gCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTa                                                                                        <  3:119411/149‑1 (MQ=255)
                                                                aGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCgg                                                                                    <  5:307146/149‑1 (MQ=255)
                                                                        gATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCgtgt                                                                            >  2:71224/1‑149 (MQ=255)
                                                                          tACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTg                                                                          <  2:48335/149‑1 (MQ=255)
                                                                          tACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTg                                                                          <  1:314933/149‑1 (MQ=255)
                                                                             gCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCa                                                                                           <  7:299936/129‑1 (MQ=255)
                                                                             gCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCa                                                                                           >  8:299936/1‑129 (MQ=255)
                                                                             gCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGc                                                                        >  2:186301/1‑148 (MQ=255)
                                                                               aCGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCAc                                                                     >  5:292034/1‑149 (MQ=255)
                                                                               aCGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCAc                                                                     >  5:90299/1‑149 (MQ=255)
                                                                                    aGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACcgtgc                                                                <  8:264808/149‑1 (MQ=255)
                                                                                        ttGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACcgtgcgtg                                                             >  1:24580/1‑148 (MQ=255)
                                                                                          gCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGca                                                           <  6:258316/148‑1 (MQ=255)
                                                                                             gTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCg                                                                                     >  4:193023/1‑119 (MQ=255)
                                                                                             gTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCg                                                                                     <  3:193023/119‑1 (MQ=255)
                                                                                             gTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAga                                                       <  4:105946/149‑1 (MQ=255)
                                                                                                gCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAgatga                                                    >  5:128385/1‑149 (MQ=255)
                                                                                                 cTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTgg                                                                         <  1:146092/127‑1 (MQ=255)
                                                                                                 cTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTgg                                                                         >  2:146092/1‑127 (MQ=255)
                                                                                                   ggTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAgatgatga                                                 >  6:320454/1‑149 (MQ=255)
                                                                                                    gTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAgatgatgat                                                >  2:12062/1‑149 (MQ=255)
                                                                                                     tCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATc                                               <  6:292034/149‑1 (MQ=255)
                                                                                                        gcagcaAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCtt                                             >  6:2527/1‑148 (MQ=255)
                                                                                                        gcagcaAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCtt                                             <  5:2527/148‑1 (MQ=255)
                                                                                                        gcagcaAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTg                                            <  1:186301/149‑1 (MQ=255)
                                                                                                            caAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTg                                                                          <  4:113280/115‑1 (MQ=255)
                                                                                                            caAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTg                                                                          >  3:113280/1‑115 (MQ=255)
                                                                                                                     cAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCgct                               >  7:188637/1‑149 (MQ=255)
                                                                                                                         gggACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTc                           <  8:188637/149‑1 (MQ=255)
                                                                                                                         gggACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTc                           <  6:90299/149‑1 (MQ=255)
                                                                                                                                cATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCTCCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCg                          <  5:171601/143‑1 (MQ=255)
                                                                                                                                cATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACcgt                                                                  >  7:144297/1‑103 (MQ=255)
                                                                                                                                cATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACcgt                                                                  <  8:144297/103‑1 (MQ=255)
                                                                                                                                cATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCg                          >  6:171601/1‑143 (MQ=255)
                                                                                                                                cATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCg                    >  6:102104/1‑149 (MQ=255)
                                                                                                                                 aTTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCGc                   >  8:200773/1‑149 (MQ=255)
                                                                                                                                  ttGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCCCAGATGATGATCTTGCTCATGCGCCGCTGCTCGTACTCGcc                  >  6:305346/1‑149 (MQ=255)
                                                                                                                                  ttGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCGcc                  <  6:128385/149‑1 (MQ=255)
                                                                                                                                    gATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGc                                                                                      >  7:157670/1‑79 (MQ=255)
                                                                                                                                    gATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGc                                                                                      <  8:157670/79‑1 (MQ=255)
                                                                                                                                      tCGCTGTCAGCCTGGCAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTc                           <  6:63978/136‑1 (MQ=255)
                                                                                                                                      tCGCTGTCAGCCTGGCAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTc                           >  5:63978/1‑136 (MQ=255)
                                                                                                                                      tCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCGCCGGCa              >  7:123876/1‑149 (MQ=255)
                                                                                                                                           gCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCGCCGGCATTGAg         <  2:24580/149‑1 (MQ=255)
                                                                                                                                           gCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCGCCGGCATTGAg         >  4:173139/1‑149 (MQ=255)
                                                                                                                                            ccAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCg                                                                                     <  4:28536/72‑1 (MQ=255)
                                                                                                                                            ccAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCg                                                                                     >  3:28536/1‑72 (MQ=255)
                                                                                                                                            ccAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCGCCGGCATTGAg         <  5:174683/148‑1 (MQ=255)
                                                                                                                                                    aCAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCGCCGGCATTGAGCACCTgg  <  2:91216/147‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GGTGTCGTACGAGCCGGTGCTGACCGAAAGGATCGGCAGGCCGGCCTGCAGGGCACCGCGGCACAGGCCAAGGATACGCACGTCAGGCTTGCTGTCGCTGGTCAGCAGCAAGCCGGCCAGGGGGACGCCATTGATCGCTGCCAGGCTGACAGCGAGGATGATGTCGTCGCGGTCGCCCGGGGTTACCACCAGGGTGCCCGAGGTCAGTAGCGGGACCGTGTTGGCCACCGTGCGTGCACAGATGATGATCTTGCTCATGCGCCGCTGCTCGTAATCGCCGGCATTGAGCACCTGG  >  NC_002947/892749‑893043

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