Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,554,551 A→C 8.4% V140G (GTC→GGC)  PP_4041 ← glycoside hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,554,5510AC 8.4% 177.3 / 7.5 70V140G (GTC→GGC) PP_4041glycoside hydrolase
Reads supporting (aligned to +/- strand):  ref base A (28/36);  new base C (2/4);  total (30/40)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

AGTCAGGCTGCACCGGCAAGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCACCACCGAATTGACTTCGCCCAGGGG  >  NC_002947/4554402‑4554696
                                                                                                                                                     |                                                                                                                                                 
aGTCAGGCTGCACCGGCAAGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGACCGGCGa                                                                                                                                                   >  1:947501/1‑150 (MQ=255)
   cAGGCTGCACCGGCAAGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACgg                                                                                                                                                <  1:639114/150‑1 (MQ=255)
       cTGCACCGGCAAGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGGGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGAGGGCACt                                                                                                                                            <  2:422683/150‑1 (MQ=255)
       cTGCACCGGCAAGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGTGCCGTCGCCGGCCCt                                                                                                                                            <  1:175747/150‑1 (MQ=255)
       cTGCACCGGCAAGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACt                                                                                                                                            >  2:1797387/1‑150 (MQ=255)
       cTGCACCGGCAAGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTCTGGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACt                                                                                                                                            <  1:497495/150‑1 (MQ=255)
         gctccGGCAAGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGTGCCGTCGCCGGCCCt                                                                                                                                            <  1:1019643/145‑1 (MQ=255)
         gCACCGGCAAGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGTGCCGTCGCCGGCCCt                                                                                                                                            >  2:1019643/1‑148 (MQ=255)
         gCACCGGCAAGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGCGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGa                                                                                                                                          <  2:155947/150‑1 (MQ=255)
          cACCGGCAAGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAg                                                                                                                                         <  1:741634/150‑1 (MQ=255)
              ggCAAGTGCGAGGGGTGCTAACACCGGCAGTGGAGTTGGGGCGTGCCCGTGCTGACGCGAAAGCTGCACTACATAGGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGTGCCGTCGCCGGCCCTGaccggt                                                                                                                                     <  1:329431/150‑6 (MQ=255)
               gCAAGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCa                                                                                                                                    >  1:271388/1‑150 (MQ=255)
                  aGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTg                                                                                                                                 <  1:1046015/150‑1 (MQ=255)
                  aGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTg                                                                                                                                 <  1:1797387/150‑1 (MQ=255)
                      cGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTc                                                                                                                             <  2:271388/150‑1 (MQ=255)
                        aGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGa                                                                                                                           <  1:862543/150‑1 (MQ=255)
                        aGGGCTGGAGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGa                                                                                                                           >  2:1613086/1‑150 (MQ=255)
                         gggCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAg                                                                                                                          >  2:994697/1‑150 (MQ=255)
                                    gCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCAccc                                                                                                               >  2:1006628/1‑150 (MQ=255)
                                     ccggttGAAGACCTGGCGCTGGCCCGTGGTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCa                                                                                                              <  2:180892/145‑1 (MQ=255)
                                     cAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTCAGCTGGCACGGTGCCGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCa                                                                                                              <  1:797720/150‑1 (MQ=255)
                                       ggCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGTTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAgg                                                                                                            <  1:1018999/150‑1 (MQ=255)
                                              gACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGt                                                                                                                      <  2:722739/133‑1 (MQ=255)
                                              gACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGt                                                                                                                      >  1:722739/1‑133 (MQ=255)
                                                      gcTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGt                                                                                             >  2:1218443/1‑150 (MQ=255)
                                                       cTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTa                                                                                            <  2:1624036/150‑1 (MQ=255)
                                                         ggCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACGGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCCGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGt                                                                                          >  2:506554/1‑150 (MQ=255)
                                                          gCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGcc                                                                                               <  2:733601/144‑1 (MQ=255)
                                                          gCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGcc                                                                                               >  1:733601/1‑144 (MQ=255)
                                                            ccGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGa                                                                                       >  1:136149/1‑150 (MQ=255)
                                                                     cgcgAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGa                                                                              <  1:749233/150‑1 (MQ=255)
                                                                      gcgAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCTGGGGCACTGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGAc                                                                             >  2:1133380/1‑150 (MQ=255)
                                                                       cgAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGTCGTAGTCGAAGCGCATGAcc                                                                            <  1:1613086/150‑1 (MQ=255)
                                                                         aaaGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAg                                                                          >  1:265249/1‑150 (MQ=255)
                                                                         aaaGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAg                                                                          <  1:994697/150‑1 (MQ=255)
                                                                           aGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGt                                                                        <  2:739754/150‑1 (MQ=255)
                                                                              ggCAGCCGGTATGGTGGTCGCGGGCTGGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGCGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTccc                                                                     >  2:817584/1‑149 (MQ=255)
                                                                                      gTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGaaa                                                             <  2:1412387/150‑1 (MQ=255)
                                                                                      gTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGaaa                                                             <  2:6362/150‑1 (MQ=255)
                                                                                      gTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGAACTGGGTCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGaaa                                                             <  1:1453500/150‑1 (MQ=255)
                                                                                         tggtggTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATgg                                                          <  1:361020/150‑1 (MQ=255)
                                                                                            tggtCGCGGGCTTGCGGGTGCACGGTGCGGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGtt                                                       >  1:449604/1‑150 (MQ=255)
                                                                                              gtCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTc                                                     <  2:136149/150‑1 (MQ=255)
                                                                                               tCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCg                                                    <  1:865408/150‑1 (MQ=255)
                                                                                                 gcgGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCCACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGcc                                                  <  1:1870340/150‑1 (MQ=255)
                                                                                                  ggggCTTGCGGGGGGACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGATGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCa                                                 <  2:140060/149‑1 (MQ=255)
                                                                                                  cgGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCa                                                 <  2:265249/150‑1 (MQ=255)
                                                                                                    ggCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGa                                               <  1:252030/150‑1 (MQ=255)
                                                                                                      cTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACa                                             >  1:682807/1‑150 (MQ=255)
                                                                                                         gCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGc                                          <  1:1218443/150‑1 (MQ=255)
                                                                                                           gggggCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCt                                        >  1:278898/1‑150 (MQ=255)
                                                                                                           gggggCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCt                                        >  1:273578/1‑150 (MQ=255)
                                                                                                                   ggTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGAt                                >  1:1299888/1‑150 (MQ=255)
                                                                                                                    gTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATg                               <  2:273578/150‑1 (MQ=255)
                                                                                                                     tgctgcACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATgc                              <  1:652383/150‑1 (MQ=255)
                                                                                                                          cacaGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGGCAGACAAGCCTTCGACGATGCGcacc                         >  2:747728/1‑150 (MQ=255)
                                                                                                                          cacaGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGcacc                         >  2:749062/1‑150 (MQ=255)
                                                                                                                          cacaGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGcacc                         >  2:1477631/1‑150 (MQ=255)
                                                                                                                           acaGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGcacca                        >  1:1088069/1‑150 (MQ=255)
                                                                                                                            caGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGcaccac                       >  1:444964/1‑150 (MQ=255)
                                                                                                                             aGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAACCGGTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGGCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCaccacc                      >  1:347345/1‑150 (MQ=255)
                                                                                                                                 tCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCACCACCGAAt                  >  1:869783/1‑150 (MQ=255)
                                                                                                                                     ccGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCACCACCGAATTGAc              >  1:358407/1‑150 (MQ=255)
                                                                                                                                        gTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCACCACCGAATTGACTTc           <  1:1133380/150‑1 (MQ=255)
                                                                                                                                        ctCAGTGCCGTCGCCGGCCCTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCa                                                 >  2:1004502/2‑112 (MQ=255)
                                                                                                                                        ctCAGTGCCGTCGCCGGCCCTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCa                                                 <  1:1004502/111‑1 (MQ=255)
                                                                                                                                          cAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCACCACCGAATTGACTTCGc         <  1:1189353/150‑1 (MQ=255)
                                                                                                                                          cAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCACCACCGAATTGACTTCGc         <  1:1189365/150‑1 (MQ=255)
                                                                                                                                          cAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCACCACCGAATTGACTTCGc         <  2:1577192/150‑1 (MQ=255)
                                                                                                                                             ggCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGGAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCACCACCGAATTGACTTCGCCCa      >  2:167865/1‑150 (MQ=255)
                                                                                                                                                 ggCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGAAGATGCGCACCACCGAATTGAATTCGCCCAgggg  >  1:650004/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
AGTCAGGCTGCACCGGCAAGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCACCACCGAATTGACTTCGCCCAGGGG  >  NC_002947/4554402‑4554696

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