Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 806,078 T→G 5.6% M418L (ATG→CTG)  PP_0693 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947806,0780TG 5.6% 221.9 / 4.1 71M418L (ATG→CTG) PP_0693hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (39/28);  new base G (2/2);  total (41/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

CGATGGCCAGGCCTGGGGCCTCGAAGCGTTCGCCCTTGAAGCCGTCGAGCACCGCTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGT  >  NC_002947/805931‑806226
                                                                                                                                                   |                                                                                                                                                    
cGATGGCCAGGCCTGGGGCCTCGAAGCGTTCGCCCTTGAAGCCGTCGAGCACCGCTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATc                                                                                                                                                     <  2:158462/149‑1 (MQ=255)
  aTGGCCAGGCCTGGGGCCTCGACGCGTTCGCCCTTGAAGCCGTCGAGCACCGCTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAg                                                                                                                                                   >  8:183317/1‑149 (MQ=255)
      ccAGGCCTGGGGCCTCGAAGCGTTCGCCCTTGAAGCCGTCGAGCACCGCTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCgg                                                                                                                                               <  8:85045/149‑1 (MQ=255)
      ccAGGCCTGGGGCCTCGAAGCGTTCGCCCTTGAAGCCGTCGAGCACCGCTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCgg                                                                                                                                               >  3:276946/1‑149 (MQ=255)
      ccAGGCCTGGGGCCTCGAAGCGTTCGCCCTTGAAGCCGTCGAGCACCGCTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCgg                                                                                                                                               >  6:183711/1‑149 (MQ=255)
        aGGCCTGGGGCCTCGAAGCGTTCGCCCTTGAAGCCGTCGAGCACCGCTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTc                                                                                                                                             <  5:99978/149‑1 (MQ=255)
          gCCTGGGGCCTCGAAGCGTTCGCCCTTGAAGCCGTCGAGCACCGCTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGa                                                                                                                                           >  3:268941/1‑149 (MQ=255)
            cTGGGGCCTCGAAGCGTTCGCCCTTGAAGCCGTCGAGCACCGCTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACg                                                                                                                                         >  6:54115/1‑149 (MQ=255)
            cTGGGGCCTCGAAGCGTTCGCCCTTGAAGCCGTCGAGCACCGCTTCCAGCGATTTTACACACAGGTCGCTGTAGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACg                                                                                                                                         >  7:186994/1‑149 (MQ=255)
                  ccTCGAAGCGTTCGCCCTTGAAGCCGTCGAGCACCGCTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCctgc                                                                                                                                   >  8:253482/1‑149 (MQ=255)
                     cGAAGCGTTCGCCCTTGAAGCCGTCGAGCACCGCTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGGCGTTGAACAGACGCATCAGGAGGTCGACGTCCtgctgc                                                                                                                                >  8:2716/1‑149 (MQ=255)
                       aaGCGTTCGCCCTTGAAGCCGTCGAGCACCGCTTCCAGCGATTTTACCAACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCAGCAGGACGTCGACCGCCTGATGCgt                                                                                                                              <  4:79181/149‑2 (MQ=255)
                          cGTTCGCCCTTGAAGCCGTCGAGCACCGCTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCATCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGaa                                                                                                                           <  8:141537/149‑1 (MQ=255)
                          cGTTCGCCCTTGAAGCCGTCGAGCACCGCTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGaa                                                                                                                           <  5:229267/149‑1 (MQ=255)
                            ttCGCCCTTGAAGCCGTCGAGCACCGCTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAAct                                                                                                                         >  8:238484/1‑149 (MQ=255)
                                 ccaagAAGCCGTCGAGCACCGCATCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGAGAACTCCCAGCGGCAGGCTGAACAGCGAGCCGTAGAACAGACGCATCAGGCGGTCGACGTCCAGCTGCTAGAACTCTTcg                                                                                                                    >  2:331766/5‑149 (MQ=255)
                                 ccTTGAAGCCGTCGAGCACCGCTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTcg                                                                                                                    >  6:351769/1‑149 (MQ=255)
                                 ccTTGAAGCCGTCGAGCACCGCTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTcg                                                                                                                    <  5:54115/149‑1 (MQ=255)
                                   ttGAAGCCGTCGAGCACCGCTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTc                                                                                                                     <  3:315103/146‑1 (MQ=255)
                                   ttGAAGCCGTCGAGCACCGCTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTc                                                                                                                     >  4:315103/1‑146 (MQ=255)
                                       aGCCGTCGAGCACCGCTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAgg                                                                                                              >  8:151553/1‑149 (MQ=255)
                                          cGTCGAGCACCGCTTCCAGCGATTTTACGCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCgg                                                                                                           >  1:207745/1‑149 (MQ=255)
                                           gTCGAGCACCGCTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCggg                                                                                                          <  8:186994/149‑1 (MQ=255)
                                           gTCGAGCACCGCTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCggg                                                                                                          >  8:267343/1‑149 (MQ=255)
                                               aGCACCGCTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATa                                                                                                      <  5:183711/149‑1 (MQ=255)
                                                      cTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCtgttg                                                                                               >  5:118348/1‑149 (MQ=255)
                                                         ccAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCtgttgtt                                                                                             <  4:276946/148‑1 (MQ=255)
                                                         ccAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCtgttgt                                                                                              <  7:151553/147‑1 (MQ=255)
                                                          cAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCtgttgttgt                                                                                           >  1:294586/1‑149 (MQ=255)
                                                          cAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCtgttgttgt                                                                                           <  6:118348/149‑1 (MQ=255)
                                                          cAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCtgttgttgt                                                                                           >  6:282451/1‑149 (MQ=255)
                                                           aGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGGtgttgttgt                                                                                           >  4:22866/1‑148 (MQ=255)
                                                            gCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCtgttgttgt                                                                                           >  8:320333/1‑147 (MQ=255)
                                                                     cccACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATcc                                                                                <  2:294586/149‑1 (MQ=255)
                                                                        acaGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAgg                                                                                                                                                  <  1:235650/80‑1 (MQ=255)
                                                                        acaGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAgg                                                                                                                                                  >  2:235650/1‑80 (MQ=255)
                                                                               gCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAgcggggcc                                                                       <  2:207745/148‑7 (MQ=255)
                                                                                 tGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTg                                                                                    <  8:132624/133‑1 (MQ=255)
                                                                                 tGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTg                                                                                    >  7:132624/1‑133 (MQ=255)
                                                                                           cTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGAGGTCCTGCTGCGAGAACGCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCttaa                                                          <  7:238484/149‑3 (MQ=255)
                                                                                              gATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAg                                                         <  4:268941/147‑1 (MQ=255)
                                                                                              gATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAg                                                         <  5:282451/147‑1 (MQ=255)
                                                                                                   cgcgATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCtgt                                                                                                 <  6:102113/102‑1 (MQ=255)
                                                                                                   cgcgATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCtgt                                                                                                 >  5:102113/1‑102 (MQ=255)
                                                                                                    gcgATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGttttt                                                 <  5:351769/149‑1 (MQ=255)
                                                                                                      gATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGAGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTc                                                >  2:302911/1‑148 (MQ=255)
                                                                                                        tctcCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTc                                                >  7:305951/1‑146 (MQ=255)
                                                                                                        tctcCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCa                                              <  8:305951/148‑1 (MQ=255)
                                                                                                         ctcCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAgg                                            >  3:67419/1‑149 (MQ=255)
                                                                                                          tcCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTc                                                <  8:53121/144‑1 (MQ=255)
                                                                                                          tcCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTc                                                >  7:53121/1‑144 (MQ=255)
                                                                                                           cccAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCGGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTc                                          >  8:161873/1‑149 (MQ=255)
                                                                                                            ccAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCCGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCAAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTcg                                         >  1:262038/1‑149 (MQ=255)
                                                                                                             cAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTcgc                                        >  3:175341/1‑149 (MQ=255)
                                                                                                             cAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTcgc                                        <  7:320333/149‑1 (MQ=255)
                                                                                                                           aCAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCGGGTc                                          <  7:45070/133‑1 (MQ=255)
                                                                                                                           aCAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCGGGTc                                          >  8:45070/1‑133 (MQ=255)
                                                                                                                            cAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTCGCGTACGCGCtcgtc                           >  3:361123/1‑147 (MQ=255)
                                                                                                                                      ttGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTCGCGTACGCGCTCGTCGAGGTCTTCgg                >  5:360494/1‑148 (MQ=255)
                                                                                                                                        gAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCAc             >  8:345014/1‑149 (MQ=255)
                                                                                                                                         aaCAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCAt                                                                                                       <  7:260198/58‑1 (MQ=255)
                                                                                                                                         aaCAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCAt                                                                                                       >  8:260198/1‑58 (MQ=255)
                                                                                                                                         aaCAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCAc             <  7:267343/148‑1 (MQ=255)
                                                                                                                                            aGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGaa         <  6:360494/149‑1 (MQ=255)
                                                                                                                                            aGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGaa         >  5:110565/1‑149 (MQ=255)
                                                                                                                                              tcGCAGCAGGACGTCGACCGCCTGATGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTCGCGTAc                                    <  4:127418/119‑1 (MQ=255)
                                                                                                                                              tcGCAGCAGGACGTCGACCGCCTGATGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTCGCGTAc                                    >  3:127418/2‑120 (MQ=255)
                                                                                                                                              tcGCAGCAGGACGTCGACCGCCTGATGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACt       >  7:222999/2‑149 (MQ=255)
                                                                                                                                              aCGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACt       >  8:162741/1‑149 (MQ=255)
                                                                                                                                              aCGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACt       >  1:185351/1‑149 (MQ=255)
                                                                                                                                                  aTCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCg   >  3:230009/1‑149 (MQ=255)
                                                                                                                                                   tCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCCTGACCTGCTTCAGGCGTTTTTCCAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGt  >  3:235690/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGATGGCCAGGCCTGGGGCCTCGAAGCGTTCGCCCTTGAAGCCGTCGAGCACCGCTTCCAGCGATTTTACCCACAGGTCGCTGTCGTCCAGCTCGATACCGCGATCTCCCAGCGGCAGGCTGAACAGCGCGCCGTTGAACAGACGCATCAGGCGGTCGACGTCCTGCTGCGAGAACTCTTCGCGCAGGCGGGCATAGCTGTTGTTGTCAGCGTGATCCAGCTGCTGCTTGACCTGCTTCAGGCGTTTTTCCAGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGT  >  NC_002947/805931‑806226

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