Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,554,543 G→T 8.3% P143T (CCT→ACT)  PP_4041 ← glycoside hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,554,5430GT 8.3% 214.7 / 5.5 72P143T (CCT→ACT) PP_4041glycoside hydrolase
Reads supporting (aligned to +/- strand):  ref base G (30/36);  new base T (2/4);  total (32/40)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

GCAGTCAGGCTGCACCGGCAAGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCACCACCGAATTGACTTCGCCCA  >  NC_002947/4554400‑4554692
                                                                                                                                               |                                                                                                                                                     
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                                                                                           <  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:6362/150‑1 (MQ=255)
                                                                                        gTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGaaa                                                         <  2:1412387/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)
                                                                                                                                               |                                                                                                                                                     
GCAGTCAGGCTGCACCGGCAAGTGCGAGGGCTGGTGACGCAGGCTGAAGACCTGGCGCTGGCCCGTGCTGACGCGAAAGCGGCAGCCGGTATGGTGGTCGCGGGCTTGCGGGGGCACGGTGCTGCACAGGATCAACCGGTCAGGGCCGGCGACGGCACTGAGGGTCAGTGGGTCGAGCCGTTCCACCCAGGGCACGCTGCGGCCGTAGTCGAAGCGCATGACCAGGTCCATTTCGAAATGGGTTTCGCCAGACAAGCCTTCGACGATGCGCACCACCGAATTGACTTCGCCCA  >  NC_002947/4554400‑4554692

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