Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,554,548 G→T 7.9% A141D (GCC→GAC)  PP_4041 ← glycoside hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,554,5480GT 7.9% 228.1 / 4.6 73A141D (GCC→GAC) PP_4041glycoside hydrolase
Reads supporting (aligned to +/- strand):  ref base G (31/36);  new base T (2/4);  total (33/40)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

GCAGTCAGGCTGCACCGGCAAGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCACCACCGAATTGACTTCGCCCAGGGG  >  NC_002947/4554400‑4554696
                                                                                                                                                    |                                                                                                                                                    
gCANTCAGGCTGCACCGGCAAGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGc                                                                                                                                                     >  2:1189353/1‑150 (MQ=255)
gCANTCAGGCTGCACCGGCAAGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGc                                                                                                                                                     >  2:1189365/1‑150 (MQ=255)
gCAGTCAGGCTGCACCGGCAAGTGCGAGGGCTTGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCGTGCGGGGGCGCGGTCCTGCACAGGATCAACCGGTCAGGGCCGGc                                                                                                                                                     >  2:865408/1‑150 (MQ=255)
  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:1797387/150‑1 (MQ=255)
                    aGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTg                                                                                                                                 <  1:1046015/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                                                                                               >  1:733601/1‑144 (MQ=255)
                                                            gCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGcc                                                                                               <  2:733601/144‑1 (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:273578/1‑150 (MQ=255)
                                                                                                             gggggCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCt                                        >  1:278898/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:1477631/1‑150 (MQ=255)
                                                                                                                            cacaGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGcacc                         >  2:749062/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)
                                                                                                                                                    |                                                                                                                                                    
GCAGTCAGGCTGCACCGGCAAGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCACCACCGAATTGACTTCGCCCAGGGG  >  NC_002947/4554400‑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: ‑