Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,161,249 A→T 5.1% Q26L (CAG→CTG)  PP_2776 → homocysteine S‑methyltransferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,161,2490AT 5.1% 233.4 / 4.0 79Q26L (CAG→CTG) PP_2776homocysteine S‑methyltransferase family protein
Reads supporting (aligned to +/- strand):  ref base A (28/45);  new base T (2/2);  total (32/47)
Fisher's exact test for biased strand distribution p-value = 6.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.66e-01

GGTCGCCGCCAATCAGGCGACCGTCGATCGTTGGTTGAGTGCCACTGGCACCCAGGTGGCTGCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGTTTCACATTGGCGAGGAACGCTTTGCGC  >  NC_002947/3161107‑3161395
                                                                                                                                              |                                                                                                                                                  
ggTCGCCGCCAATCAGGCGACCGTCGATCGTTGGTTGAGTGCCACTGGCACCCAGGTGGCTGCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGaa                                                                                                                                             <  1:1739282/150‑1 (MQ=255)
  tCGCCGCCAATTAGGCGACCGTCGATCGTTGGTTGAGTGCCACTGGCACCCAGGTGGCTGCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATg                                                                                                                                           >  1:443507/1‑150 (MQ=255)
  tCGCCGCCAATCAGGCGACCGTCGATCGTTGGTTGAGTGCCACTGGCACCCAGGTGGCTGCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATg                                                                                                                                           <  2:1443575/150‑1 (MQ=255)
     ccgccCATCAGCCGACCGTCGATCGTTGGTTGCGTACCACTGGCACCCAGCTGGCTGCCCAACGACTCGGCGTGCCTCGACTTTGCTCGCTGCCAGGCCGGGACGAGAACTCTATCGGATCGGGGTACCATTCCCGCCGCCGGCATGGTc                                                                                                                                        >  2:592653/1‑150 (MQ=255)
     ccgccAATCAGGCGACCGTCGATCGTTGGTTGAGTGCCACTGGCACCCAGGTGGCTGCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCTGCCGGaag                                                                                                                                            >  1:666564/1‑145 (MQ=255)
     ccgccAATCAGGCGACCGTCGATCGTTGGTTGAGTGCCACTGGCACCCAGGTGGCTGCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCTGCCGGaag                                                                                                                                            <  2:666564/146‑2 (MQ=255)
     ccgccAATCAGGCGACCGTCGATCGTTGGTTGAGTGCCACTGGCACCCAGGTGGCTGCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTc                                                                                                                                        <  1:537256/150‑1 (MQ=255)
     ccgccAATCAGGCGACCGTCGATCGTTGGTTGAGTGCCACTGGCACCCAGGTGGCTGCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTc                                                                                                                                        <  2:1674572/150‑1 (MQ=255)
        ccAATCAGGCGACCGTCGATCGTTGGTTGAGTGCCACTGGCACCCAGGTGGCTGCCCAATGACTCGGCGTGCATGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGc                                                                                                                                     >  1:1564324/1‑150 (MQ=255)
            tCAGGCGACCGTCGATCGTTGGTTGAGTGCCACTGGCACCCAGGTTGCTGCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGTCAGCCGGAATGGTCGGCACTg                                                                                                                                 <  2:166950/150‑1 (MQ=255)
            tCAGGCGACCGTCGATCGTTGGTTGAGTGCCACTGGCACCCAGGTGGCTGCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTg                                                                                                                                 >  1:1879097/1‑150 (MQ=255)
                 cGACCGTCGATCGTTGGTTGAGTGCCACTGGCACCCAGGTGGCTGCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCt                                                                                                                            <  1:215409/150‑1 (MQ=255)
                   aCCGTCGATCGTTGGTTGAGTGCCACTGGCACCCAGGTGGCTGCCCAATGACTCGGCGTGCCTGTACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGt                                                                                                                          <  1:646104/150‑1 (MQ=255)
                     cGTCGATCGTTGGTTGAGTGCCACTGGCACCCAGGTGGCTGCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCGGCCGGAATGGTCGGCACTGGCGCTGTTc                                                                                                                        >  1:1716524/1‑150 (MQ=255)
                     cGTCGATCGTTGGTTGAGTGCCACTGGCACCCAGGTGGCTGCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTc                                                                                                                        >  1:1428389/1‑150 (MQ=255)
                       tCGATCGTTGGTTGAGTGCCACTGGCACCCAGGTGGCTGCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGtt                                                                                                                         >  1:1531651/1‑147 (MQ=255)
                       tCGATCGTTGGTTGAGTGCCACTGGCACCCAGGTGGCTGCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGtt                                                                                                                         <  2:1531651/147‑1 (MQ=255)
                       tCGATCGTTGGTTGAGTGCCACTGGCACCCAGGTGGCTGCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGa                                                                                                                      >  2:750070/1‑150 (MQ=255)
                        cGATCGTTGGTTGAGTGCCACTGGCACCCAGGTGGCTGCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGaa                                                                                                                     <  2:616412/150‑1 (MQ=255)
                                ggttgAGTGCCACTGGCACCCAGGTGGCTGCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGtt                                                                                                                         >  1:1012032/1‑138 (MQ=255)
                                ggttgAGTGCCACTGGCACCCAGGTGGCTGCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGtt                                                                                                                         <  2:1012032/138‑1 (MQ=255)
                                  ttgAGTGCCACTGGCACCCAGGTGGCTGCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAAt                                                                                                           >  2:1890964/1‑150 (MQ=255)
                                  ctgggTGCCAGTGGCACCCAGGTGGCTGCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTTGCGCTGTTCGGAGCGCCCGAAt                                                                                                           >  2:750985/5‑150 (MQ=255)
                                   tgAGTGCCACTGGCACCCAGGTGGCTGCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACCATGCCTCAAGCGGATGGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATa                                                                                                          <  1:705045/150‑1 (MQ=255)
                                      gTGCCACTGGCACCCAGGTGGCTGCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGt                                                                                                       <  1:178242/150‑1 (MQ=255)
                                       tGCCACTGGCACCCAGGTGGCTGCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTg                                                                                                      <  1:1412753/150‑1 (MQ=255)
                                        gCCACTGGCACCCAGGTGGCTGCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGa                                                                                                     <  2:1823177/150‑1 (MQ=255)
                                         ccACTGGCACCCAGGTGGCTGCCCAGTGGCTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGAc                                                                                                    >  1:257811/1‑150 (MQ=255)
                                          cACTGGCACCCAGGTGGCTGCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACg                                                                                                   >  2:1871350/1‑150 (MQ=255)
                                          cACTGGCACCCAGGTGGCTGCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACg                                                                                                   >  2:115156/1‑150 (MQ=255)
                                            cTGGCACCCAGGTGGCTGCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCa                                                                                                 >  2:866699/1‑150 (MQ=255)
                                             tGGCACCCAGGTGGCTGCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGt                                                                                                       >  1:501335/1‑143 (MQ=255)
                                             tGGCACCCAGGTGGCTGCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGt                                                                                                       <  2:501335/143‑1 (MQ=255)
                                              ggCACCCAGGTGGCTGCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAgg                                                                                               <  1:552272/150‑1 (MQ=255)
                                                cACCCAGGTGGCTGCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGtt                                                                                             <  1:866699/150‑1 (MQ=255)
                                                 aCCCAGGTGGCTGCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTc                                                                                                                        >  2:213691/1‑122 (MQ=255)
                                                 aCCCAGGTGGCTGCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTc                                                                                                                        <  1:213692/122‑1 (MQ=255)
                                                 aCCCAGGTGGCTGCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTc                                                                                                                        <  1:213691/122‑1 (MQ=255)
                                                 aCCCAGGTGGCTGCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTc                                                                                                                        >  2:213692/1‑122 (MQ=255)
                                                      ggTGGCTGCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCCGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGa                                                                                       <  1:803035/150‑1 (MQ=255)
                                                        tGGCTGCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGc                                                                                     <  2:1879097/150‑1 (MQ=255)
                                                             gCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGGTATAAGCGCCCGCATACGTGACGACAGTTCACCGAGCGTTCa                                                                                >  1:144/1‑150 (MQ=255)
                                                               ccAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCAtt                                                                              <  1:782871/150‑1 (MQ=255)
                                                                cAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCAtcg                                                                             >  1:1596954/1‑148 (MQ=255)
                                                                    gACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTgccgc                                                                         <  2:1629186/150‑1 (MQ=255)
                                                                       tCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCgg                                                                      >  2:1419939/1‑150 (MQ=255)
                                                                         ggCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCgggg                                                                    <  2:1284279/150‑1 (MQ=255)
                                                                              gCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAgg                                                               <  1:115156/150‑1 (MQ=255)
                                                                                cTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGt                                                             >  2:455324/1‑150 (MQ=255)
                                                                                       ttGCTCGATGGTGGGATGGGACGAGAACTCAACCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATcacc                                                      <  1:1637681/150‑1 (MQ=255)
                                                                                            cGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCaa                                                 <  2:83030/150‑1 (MQ=255)
                                                                                             gATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAAc                                                <  1:1189736/150‑1 (MQ=255)
                                                                                                 gtggGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCt                                            >  2:1516981/1‑150 (MQ=255)
                                                                                                 gtggGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCt                                            >  1:981601/1‑150 (MQ=255)
                                                                                                          gaCGAGAACTCAAGCGGATCGGCGCACCCTTCCGATATCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCggtgg                                   <  1:489048/150‑1 (MQ=255)
                                                                                                           aCGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCggtggt                                  <  1:1286879/150‑1 (MQ=255)
                                                                                                             gagaACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGgcgc                                                                                                                             <  2:130522/57‑1 (MQ=255)
                                                                                                             gagaACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGgcgc                                                                                                                             >  1:130522/1‑57 (MQ=255)
                                                                                                             gagaACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGc                                <  2:68825/150‑1 (MQ=255)
                                                                                                                aaatCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGGACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGt                             <  2:370757/147‑1 (MQ=255)
                                                                                                                  cTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGttt                           <  1:236411/150‑1 (MQ=255)
                                                                                                                  cTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGttt                           <  1:1871350/150‑1 (MQ=255)
                                                                                                                   tCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCt                                            >  1:414644/1‑132 (MQ=255)
                                                                                                                   tCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCt                                            <  2:414644/132‑1 (MQ=255)
                                                                                                                     aaGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTa                                                            <  2:471634/114‑1 (MQ=255)
                                                                                                                     aaGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTa                                                            >  1:471634/1‑114 (MQ=255)
                                                                                                                          gATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGTTTCACATTgg                   >  1:454845/1‑150 (MQ=255)
                                                                                                                             cGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGTTTCACATTGGCGa                <  2:335440/150‑1 (MQ=255)
                                                                                                                             cGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGTTTCACATTGGCGa                <  1:455324/150‑1 (MQ=255)
                                                                                                                             cGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGTTTCACATTGGCGa                <  1:1240697/150‑1 (MQ=255)
                                                                                                                             cGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGTTTCACATTGGCGa                <  2:1398832/150‑1 (MQ=255)
                                                                                                                             cGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTCCCAAGCGACCGAATATGTAAGACGGTCTCACCGCGCCGACATTTCCGCCTTGTCCATGTTAAGTACCACCAACAGCTACTCGGTGGTGCCGTTTCACATTGGCGa                <  2:1161872/150‑1 (MQ=255)
                                                                                                                                   accattCCGGCTGCCGGAAGGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGtt                            <  2:82281/128‑1 (MQ=255)
                                                                                                                                   aCCATTCCGGCTGCCGGAAGGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGtt                            >  1:82281/1‑132 (MQ=255)
                                                                                                                                       ttCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGTTTCACATTGGCGAGGAACGCttt      >  1:212104/1‑150 (MQ=255)
                                                                                                                                       ttCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGTTTCACATTGGCGAGGAACGCttt      >  1:1776558/1‑150 (MQ=255)
                                                                                                                                       ttCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGTTTCACATTGGCGAGGAACGCttt      <  2:1776558/150‑1 (MQ=255)
                                                                                                                                        tCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGTTTCACATTGGCGAGGAACGCTTTg     <  2:443507/150‑1 (MQ=255)
                                                                                                                                           ggCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATCCGTGCCGCAGGTTCACCGAGCGTTCATTTCCGCTGGAGCCAGCGAAATCGCGACGAACGGCTACGGGGTTGTGCCGTTTCACATTGGCGAGGAACGCTTTgcgc  <  1:976421/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
GGTCGCCGCCAATCAGGCGACCGTCGATCGTTGGTTGAGTGCCACTGGCACCCAGGTGGCTGCCCAATGACTCGGCGTGCCTGGACTTTGCTCGATGGTGGGATGGGACGAGAACTCAAGCGGATCGGCGCACCCTTCCGGCAGCCGGAATGGTCGGCACTGGCGCTGTTCGAAGCGCCCGAATACGTGACGCAGGTTCACCGAGCGTTCATTGCCGCCGGGGCCAGGGTAATCACCACCAACAGCTACGCGGTGGTGCCGTTTCACATTGGCGAGGAACGCTTTGCGC  >  NC_002947/3161107‑3161395

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: