Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 895,213 G→C 5.7% intergenic (+50/+14) PP_0777 → / ← PP_0778 glutathione peroxidase/group 1 family glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947895,2130GC 5.7% 218.4 / 3.8 70intergenic (+50/+14)PP_0777/PP_0778glutathione peroxidase/group 1 family glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base G (29/37);  new base C (2/2);  total (31/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

TGGGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTCGCGCATTCAG  >  NC_002947/895069‑895361
                                                                                                                                                |                                                                                                                                                    
tggGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCCCGcgcggc                                                                                                                                                 <  2:906728/150‑1 (MQ=255)
   gCAAGGATGGCCGGGTGCTGGCGCGTTTTTCGCCGCGCTCGGCGCCGGATGACACAGCTGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCt                                                                                                                                              <  1:21119/150‑1 (MQ=255)
    cAAGGATGGCCGGGTGCTGGCGCGTTTCTAGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGCGCCGCCCGCGCGGCACtc                                                                                                                                             <  1:299817/150‑2 (MQ=255)
     aaGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGAGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGACGCTgg                                                                                                                                            >  2:817246/1‑150 (MQ=255)
      aGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGa                                                                                                                                           <  2:1707795/150‑1 (MQ=255)
      aGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGa                                                                                                                                           <  1:1357811/150‑1 (MQ=255)
      aGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCCAGCGCCGCCCGcgcggcactcga                                                                                                                                           <  1:6544/150‑7 (MQ=255)
       ggATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGCGCCGCCCGCGCGGCACTCGAt                                                                                                                                          <  1:591078/150‑1 (MQ=39)
        gATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGATAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGAtc                                                                                                                                         >  2:1055396/1‑150 (MQ=255)
         aTGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGAtct                                                                                                                                        <  1:1309836/150‑1 (MQ=255)
          tGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGAtctc                                                                                                                                       >  2:277040/1‑150 (MQ=255)
           ggCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTca                                                                                                                                      <  2:924160/150‑1 (MQ=255)
           ggCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTca                                                                                                                                      <  2:1730290/150‑1 (MQ=255)
             ccGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTcaca                                                                                                                                    >  2:193488/1‑150 (MQ=255)
                  tacTGGCGCGTTTCTCGCCGCGCACGGCTCCGGCTGACCCAGCGGTGGTACAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCTCCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGgcgc                                                                                                                               <  2:981862/148‑1 (MQ=255)
                            tttCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAGACAtt                                                                                                                     <  2:617545/150‑1 (MQ=255)
                            tttCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACAtt                                                                                                                     >  2:1144352/1‑150 (MQ=255)
                             ttCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTg                                                                                                                    <  1:1836992/150‑1 (MQ=255)
                                tcGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCgcgg                                                                                                                                                      >  2:1090727/1‑113 (MQ=255)
                                tcGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCgcgg                                                                                                                                                      <  1:1090727/113‑1 (MQ=255)
                                 cgccgcGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCa                                                                                                                <  2:561362/150‑1 (MQ=255)
                                 cgccgcGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCa                                                                                                                <  1:277040/150‑1 (MQ=255)
                                      cgcACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACATGCGCTGCAAACATTGCGCAGAACg                                                                                                           <  2:1283373/150‑1 (MQ=255)
                                        cACGGCGCCGGATGACCCCGCGGTGGTGCAGGCCATCGAGAAACCCCTGGTCTGAAGAGTGCCTGACGTGATGTTTTGGCATCTGTGAGAGCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGcc                                                                                                         <  2:710757/150‑1 (MQ=255)
                                        cACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGcc                                                                                                         <  1:1349683/150‑1 (MQ=255)
                                         aCGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGccc                                                                                                        >  2:1471494/1‑150 (MQ=255)
                                             cgcCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGcc                                                                                                    <  1:193488/150‑1 (MQ=255)
                                                     gACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCAcc                                                                                            <  1:1314560/150‑1 (MQ=255)
                                                     gACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCAcc                                                                                            <  1:436906/150‑1 (MQ=255)
                                                       cccAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAg                                                                                          <  1:1702997/150‑1 (MQ=255)
                                                       cccAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAg                                                                                          <  1:1156727/150‑1 (MQ=255)
                                                        ccAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGTGCAGAACGCCCAGCCGTAGCACCAgg                                                                                         <  1:817246/150‑1 (MQ=255)
                                                        ccAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAgg                                                                                         <  1:998439/150‑1 (MQ=255)
                                                        ccAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGGGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAgg                                                                                         >  2:1707919/1‑150 (MQ=255)
                                                                 gtgCAGGCCATCGAGAAAGTCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCa                                                                                >  2:800177/1‑150 (MQ=255)
                                                                    cAGGCCATCGAGAAAGCCCTGGCATGAAGTGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATgc                                                                             <  1:292930/150‑1 (MQ=255)
                                                                     aGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATgcg                                                                            >  2:1641237/1‑150 (MQ=255)
                                                                      ggCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGGCGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCTCAGCCGTAGCACCAGGTTGCCGGCATgcgc                                                                           >  2:1017266/1‑150 (MQ=255)
                                                                          aTCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATTGAGTGCCGCGCGGGCGGCGCTGGACCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGtt                                                                       >  2:1037758/1‑150 (MQ=255)
                                                                          aTCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGcc                                                                                                         >  1:1131096/1‑116 (MQ=255)
                                                                          aTCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGcc                                                                                                         <  2:1131096/116‑1 (MQ=255)
                                                                           tCGAGAAAGCCCTGGACTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTc                                                                      >  1:64218/1‑150 (MQ=255)
                                                                             gagaAAGCCCTGGCCTGAAGAGTGACTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAg                                                                    >  1:257464/1‑150 (MQ=255)
                                                                               gaAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCCAGCGCCGCCCGcgcggcac                                                                                                                                               <  2:583701/73‑3 (MQ=21)
                                                                               gaAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCCAGCGCCGCCCGcgcggcac                                                                                                                                               >  1:583701/1‑71 (MQ=21)
                                                                                aaaGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATa                                                                 >  1:1001408/1‑150 (MQ=255)
                                                                                aaaGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATa                                                                 >  1:1883618/1‑150 (MQ=255)
                                                                                   gCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATAccc                                                              <  1:1228137/150‑1 (MQ=255)
                                                                                    cccTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATAcccc                                                             >  1:441221/1‑150 (MQ=255)
                                                                                     ccTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCt                                                            >  1:122129/1‑150 (MQ=255)
                                                                                      cTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTc                                                           <  2:204098/150‑1 (MQ=255)
                                                                                      cTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTc                                                           <  1:1055396/150‑1 (MQ=255)
                                                                                          ccTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAAc                                                       <  2:1157752/150‑1 (MQ=255)
                                                                                           cTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCt                                                                                                                              <  2:546143/78‑1 (MQ=255)
                                                                                           cTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCt                                                                                                                              >  1:546143/1‑78 (MQ=255)
                                                                                           cTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAAATACCCGTCGAAcg                                                      >  1:1270565/1‑149 (MQ=255)
                                                                                                      cTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATc                                           <  2:122129/150‑1 (MQ=255)
                                                                                                      cTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATc                                           >  1:861050/1‑150 (MQ=255)
                                                                                                           gTGNTGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGAGACCCCTCGAACTGTTCGACAATCGCCgg                                      >  2:1186315/1‑150 (MQ=255)
                                                                                                                gTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGc                                 >  2:345706/1‑150 (MQ=255)
                                                                                                                gTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGc                                 >  2:597496/1‑150 (MQ=255)
                                                                                                                 ttttGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCGAGcc                                >  2:944263/1‑150 (MQ=255)
                                                                                                                 ttttGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGcc                                >  1:922830/1‑150 (MQ=255)
                                                                                                                      gCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTggc                           <  1:1471494/150‑1 (MQ=255)
                                                                                                                           tgtgAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTg                      <  2:441221/150‑1 (MQ=255)
                                                                                                                           tgtgAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTg                      >  2:1292935/1‑150 (MQ=255)
                                                                                                                             tgAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGc                    <  1:1786248/150‑1 (MQ=255)
                                                                                                                              gagaTCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCa                   <  1:345706/150‑1 (MQ=255)
                                                                                                                                 atccagcgccgccCGCGCGGCACTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGc                >  1:184401/8‑150 (MQ=255)
                                                                                                                                   cGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTg              >  2:213639/1‑150 (MQ=255)
                                                                                                                                     agcgccgccCGCGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTc            >  1:713789/4‑150 (MQ=255)
                                                                                                                                     aGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGTCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTc            <  1:1230727/150‑1 (MQ=255)
                                                                                                                                       tGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGTCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTcgc          >  2:579845/1‑150 (MQ=255)
                                                                                                                                       ggCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATCCCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTcgc          <  1:301957/149‑1 (MQ=255)
                                                                                                                                        gCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACGTTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTcgcg         <  1:800177/150‑1 (MQ=255)
                                                                                                                                            cgcgGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTCGCGCAtt     <  1:1641237/150‑1 (MQ=255)
                                                                                                                                               gggcggCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCCGTCGCGCATTCAg  <  1:1144352/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TGGGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTCGCGCATTCAG  >  NC_002947/895069‑895361

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