Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,564,931 A→T 7.8% L133Q (CTG→CAG)  pitB ← phosphate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,564,9310AT 7.8% 189.2 / 4.7 64L133Q (CTG→CAG) pitBphosphate transporter
Reads supporting (aligned to +/- strand):  ref base A (30/29);  new base T (3/2);  total (33/31)
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.54e-01

GGCGCCACCATTTCAGGCCGATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTGGCCGGTGTTCACATTGATCAGCAG  >  NC_002947/1564785‑1565079
                                                                                                                                                   |                                                                                                                                                    
ggCGCCACCATTTCAGGCCGATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTggccagg                                                                                                                                                    <  8:259674/149‑1 (MQ=255)
 gcgcCACCATTTCAGGCCGATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTggccaggc                                                                                                                                                   <  1:104617/149‑1 (MQ=255)
 gcgcCACCATTTCAGGCCGATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTggccaggc                                                                                                                                                   >  2:104617/1‑149 (MQ=255)
     caccaTTTCAGGCCGATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAac                                                                                                                                               <  1:50645/149‑1 (MQ=255)
            tCAGGCCGATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAg                                                                                                                                         <  2:262123/148‑1 (MQ=255)
                   gATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAAc                                                                                                                                 <  8:328309/149‑1 (MQ=255)
                        gCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGAt                                                                                                                             >  1:167696/1‑148 (MQ=255)
                          aCCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAg                                                                                                                          >  8:126087/1‑149 (MQ=255)
                              gggCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAg                                                                                                                                         <  7:192528/130‑1 (MQ=255)
                              gggCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAg                                                                                                                                         >  8:192528/1‑130 (MQ=255)
                               ggCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCgtgt                                                                                                                     <  5:337647/149‑1 (MQ=255)
                                 cagcaACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGc                                                                                                                   <  4:171254/149‑1 (MQ=255)
                                         gCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGa                                                                                                           <  7:230052/149‑1 (MQ=255)
                                                 ggCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTCGATCAGGGCGTTGGCCTGGCCAacgc                                                                                                                                             <  1:4330/108‑3 (MQ=255)
                                                 ggCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCTGGCCAacgc                                                                                                                                             >  2:4330/1‑105 (MQ=255)
                                                     aTTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGc                                                                                                >  5:169752/1‑148 (MQ=255)
                                                      ttGGCGAAACCACAAAAGACATAACGATATCAATCGCTTTCTACCAATTGCCACCGTCGCACAGCGGAATGTCGTT‑GATCAGGGCGTTGGCAAGCCCAACACAGAGGATGGAACCGATAAGCGTGTGCGAACTGGAGGCGGGGATGcc                                                                                               >  8:52952/1‑148 (MQ=255)
                                                         gCGATACCACCAGCGAAATGTCGATATCGATCCTTTTCTGCCAGTTGACTTCGTCGCCGAGCGGAATGTTGTT‑TATAAGGGCGTTGTCAAGGACAACACCGATGTTGGAACCTATCTTTGTTTGCTAACTGGAGGCGGGGATGCCaaa                                                                                            <  7:52952/148‑1 (MQ=255)
                                                         gCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCgg                                                                                                      >  3:156102/1‑138 (MQ=255)
                                                         gCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCgg                                                                                                      <  4:156102/138‑1 (MQ=255)
                                                         gCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCaaa                                                                                            <  7:126087/148‑1 (MQ=255)
                                                           gagaCCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTa                                                                                         >  2:274055/1‑149 (MQ=255)
                                                           gagaCCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTa                                                                                         <  3:181157/149‑1 (MQ=255)
                                                            agaCCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTAc                                                                                        <  2:167696/149‑1 (MQ=255)
                                                                  ccaGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTg                                                                                  >  7:143214/1‑149 (MQ=255)
                                                                  ccaGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTg                                                                                  >  5:139390/1‑149 (MQ=255)
                                                                                         gCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGc                                                           <  1:274055/149‑1 (MQ=255)
                                                                                          cTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCg                                                            >  1:326375/1‑147 (MQ=255)
                                                                                           tttCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTc                                                                         >  8:254035/1‑133 (MQ=255)
                                                                                           tttCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTc                                                                         <  7:254035/133‑1 (MQ=255)
                                                                                           tttCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGcc                                                             >  8:143717/1‑145 (MQ=255)
                                                                                           tttCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGcc                                                             <  7:143717/145‑1 (MQ=255)
                                                                                           tttCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGcag                                                         >  2:198688/1‑149 (MQ=255)
                                                                                           tttCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGcag                                                         >  2:29758/1‑149 (MQ=255)
                                                                                           tttCTGCCAGTTGACGCCGTCGCCCAGCGGAATATCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGcag                                                         >  1:49704/1‑149 (MQ=255)
                                                                                                  cAGTTGACGCCGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCaa                                                  <  1:198688/149‑1 (MQ=255)
                                                                                                   aGTTGACGCCGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCGGGTTCCAGGTGATGGCCGCAGCCAGCacc                                                 >  4:56387/1‑147 (MQ=255)
                                                                                                   aGTTGACGCCGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAAc                                                 >  7:316075/1‑149 (MQ=255)
                                                                                                   aGTTGACGCCGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAAc                                                 >  6:146928/1‑149 (MQ=255)
                                                                                                    gTTGACGCCGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAAc                                                 >  5:354585/1‑148 (MQ=255)
                                                                                                       gACGCCGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACgag                                              <  8:316075/148‑1 (MQ=255)
                                                                                                         cGCCGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACAGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCAAGGTTCCAGGTGATGGCCGCAGCCAGCAACgag                                              >  6:112858/1‑146 (MQ=255)
                                                                                                           ccGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGAGCGAACTGGAGCCGGGGATGCCAAACTACCAGGTGCCCAGGTTCCAGGTGATGCCCGCAGCCAGCAACGAGAAGAc                                         >  6:226517/1‑149 (MQ=255)
                                                                                                             gTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCa                                       >  7:264910/1‑149 (MQ=255)
                                                                                                             gTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCa                                       >  5:228425/1‑149 (MQ=255)
                                                                                                              tCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCAt                                      <  2:326375/149‑1 (MQ=255)
                                                                                                              tCGCCCAGCGGAATATCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCAt                                      <  2:49704/149‑1 (MQ=255)
                                                                                                                 cccAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCg                                                            >  4:242396/1‑124 (MQ=255)
                                                                                                                 cccAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCg                                                            <  3:242396/124‑1 (MQ=255)
                                                                                                                 cccAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGc                                   <  6:131111/149‑1 (MQ=255)
                                                                                                                 cccAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGc                                   >  7:216255/1‑149 (MQ=255)
                                                                                                                                tt‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTGGCCGgt                    >  6:314841/1‑149 (MQ=255)
                                                                                                                                   gATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTGGCCGgtg                   <  6:169752/148‑1 (MQ=255)
                                                                                                                                   gATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCCGGGATGCCAAAGTAACAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTGGCCGgagt                  >  7:232590/1‑149 (MQ=255)
                                                                                                                                    aTCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTGGCCGgtgt                  <  8:143214/148‑1 (MQ=255)
                                                                                                                                    aTCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTGGCCGGTGtt                 <  6:139390/149‑1 (MQ=255)
                                                                                                                                     tCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTGGCCGGTGtt                 >  3:320602/1‑148 (MQ=255)
                                                                                                                                     tCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGACAAGACCATGGCCAGGCCGTGGCCGGTGTTc                <  1:29758/149‑1 (MQ=255)
                                                                                                                                       tcggtgTTGGCCTGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTGGCCGGTGTTca               >  8:35674/6‑148 (MQ=255)
                                                                                                                                       aGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACgag                                              >  1:121946/1‑117 (MQ=255)
                                                                                                                                       aGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACgag                                              <  2:121946/117‑1 (MQ=255)
                                                                                                                                       aGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTGGCCGGTGTTcac              <  6:228425/149‑1 (MQ=255)
                                                                                                                                         ggtgTTGGCCTGGCCAACGCCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGcc                                  >  8:200789/4‑127 (MQ=255)
                                                                                                                                         ggtgTTGGCCTGGCCAACGCCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGcc                                  <  7:200789/124‑1 (MQ=255)
                                                                                                                                         ggCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTGGCCGGTGTTCatat            <  8:216255/149‑1 (MQ=255)
                                                                                                                                            gTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTGGCCGGTGTTCACATTGa         <  4:320602/149‑1 (MQ=255)
                                                                                                                                             ttGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTGGCCGGTGTTCACATTGAt        >  3:154109/1‑149 (MQ=255)
                                                                                                                                                 ccaggtcaACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTGGCCGGTGTTCACATTGATcag     >  4:297853/1‑148 (MQ=255)
                                                                                                                                                  caggccaACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCACGGCCTGGCCATggtctt                     >  2:34417/1‑127 (MQ=255)
                                                                                                                                                  caggccaACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCACGGCCTGGCCATggtctt                     <  1:34417/131‑5 (MQ=255)
                                                                                                                                                   aggccaACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTGGCCGGTGTTCACATTGATcagcag  <  6:354585/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GGCGCCACCATTTCAGGCCGATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTT‑GATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTGGCCGGTGTTCACATTGATCAGCAG  >  NC_002947/1564785‑1565079

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