Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 821,682 G→C 10.8% M350I (ATG→ATC PP_0707 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947821,6820GC10.8% 224.0 / 15.9 75M350I (ATG→ATCPP_0707hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (35/31);  new base C (4/4);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

GTGGTGTTCAGTGCCTTCATCACAGGCCTGGGCCGCGCCATGCTGATGCTGCAGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGACCTGGAGCG  >  NC_002947/821536‑821831
                                                                                                                                                  |                                                                                                                                                     
gtggtgTTCAGTGCCTTCATCACAGGCCTGGGCCGCGCCATGCTGATGCTGCAGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCatgatg                                                                                                                                                    >  2:773966/1‑150 (MQ=255)
       tCAGTGCCTTCATCACAGGCCTGGGCCGCGCCATGCTGATGCTGCAGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATCATGACCAtca                                                                                                                                             <  1:699216/150‑3 (MQ=255)
         aGTGCCTTCATCACAGGCCTGGGCCGCGCCATGCTGATGCTGCAGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAg                                                                                                                                           <  1:213512/150‑1 (MQ=255)
             ccTTCATCACAGGCCTGGGCCGCGCCATGCTGATGCTGCAGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCa                                                                                                                                       <  2:1031533/150‑1 (MQ=255)
              cTTCATCACAGGCCTGGGCCGCGCCATGCTGATGCTGCAGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAgc                                                                                                                                          <  1:536622/146‑1 (MQ=255)
              cTTCATCACAGGCCTGGGCCGCGCCATGCTGATGCTGCAGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAgc                                                                                                                                          >  2:536622/1‑146 (MQ=255)
               ttCATCACAGGCCTGGGCCGCGCCATGCTGATGCTGCAGCGTCCGTCCTGGCGGTTGCCGGCGGTCGACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATg                                                                                                                                     >  2:106041/1‑150 (MQ=255)
               ttCATCACAGGCCTGGGCCGCGCCATGCTGATGCTGCAGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATg                                                                                                                                     >  2:1025597/1‑150 (MQ=255)
                   tcaCAGGCCTGGGCCGCGCCATGCTGATGCTGCAGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACa                                                                                                                                 >  2:677471/1‑150 (MQ=255)
                     acaGGCGTGGGCCGCGGCATGCTGATGCTGCAGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGTCGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATTATGACCATGGAGCGCATGAACAGc                                                                                                                               <  2:440986/150‑1 (MQ=255)
                     acaGGCCTGGGCCGCGCCATGCTGATGCTGCAGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGc                                                                                                                               >  2:507848/1‑150 (MQ=255)
                              ggCCGCGCCATGCTGATGCTGCAGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCGGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCAAGGATCGCGTGAACAGtgtggttctt                                                                                                                      >  2:561286/1‑144 (MQ=255)
                                ccGCGCCATGCTGATGCTGCAGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGc                                                                                                                    >  1:1018874/1‑150 (MQ=255)
                                    gcCATGCTGATGCTGCAGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGc                                                                                                                >  2:231298/1‑150 (MQ=255)
                                       atgctgatgctgCAGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGATGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTg                                                                                                             >  1:243427/1‑150 (MQ=255)
                                         gctgatgctgCAGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGgc                                                                                                           <  1:657606/150‑1 (MQ=255)
                                            gatgctgCAGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCt                                                                                                        >  1:512738/1‑150 (MQ=255)
                                            gatgctgCAGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCt                                                                                                        <  1:424537/150‑1 (MQ=255)
                                                ctgcAGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACg                                                                                                    >  1:430021/1‑150 (MQ=255)
                                                  gcAGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGt                                                                                                  <  1:561286/150‑1 (MQ=255)
                                                    aGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCTCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTgg                                                                                                <  1:812763/150‑1 (MQ=255)
                                                             ccTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAATAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACg                                                                                       >  1:703191/1‑150 (MQ=255)
                                                              cTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACgg                                                                                      >  2:411805/1‑150 (MQ=255)
                                                               tGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACTATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACggg                                                                                     <  2:292633/150‑1 (MQ=255)
                                                               tGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACggg                                                                                     <  1:150577/150‑1 (MQ=255)
                                                                 gCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCt                                                                                   >  1:509938/1‑150 (MQ=255)
                                                                         cGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCgcgc                                                                           <  2:512738/150‑1 (MQ=255)
                                                                          ggCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCt                                                                          >  2:448065/1‑150 (MQ=255)
                                                                             gATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCtggt                                                                       >  1:399019/1‑150 (MQ=255)
                                                                             gATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCtggt                                                                       >  1:1200080/1‑150 (MQ=255)
                                                                                 cACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCtggtggtg                                                                   >  1:9233/1‑150 (MQ=255)
                                                                                  acgacgAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCAt                                                                                                                                               >  2:536080/1‑73 (MQ=255)
                                                                                  acgacgAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCAt                                                                                                                                               <  1:536080/73‑1 (MQ=255)
                                                                                  acgacgAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCtggtggtgg                                                                  <  2:681577/150‑1 (MQ=255)
                                                                                    gactaAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGgcg                                                                >  2:903301/1‑150 (MQ=255)
                                                                                     acgaAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGgcgc                                                               <  2:430021/150‑1 (MQ=255)
                                                                                           tCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCt                                                         <  2:509938/150‑1 (MQ=255)
                                                                                              ccAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACa                                                                                                                                 <  1:747154/75‑1 (MQ=255)
                                                                                              ccAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACa                                                                                                                                 >  2:747154/1‑75 (MQ=255)
                                                                                                aGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCa                                                    <  2:159998/150‑1 (MQ=255)
                                                                                                 gTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAg                                                   <  1:773966/150‑1 (MQ=255)
                                                                                                    cgcTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGgcg                                                <  2:703191/150‑1 (MQ=255)
                                                                                                     gctgggctggTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGgcgc                                               <  1:411805/150‑1 (MQ=255)
                                                                                                       tgggctggTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCGGTGCAGgcgcgc                                             <  2:1200080/150‑1 (MQ=255)
                                                                                                        gggctggTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCgct                                            >  1:1088318/1‑150 (MQ=255)
                                                                                                           ctggTTCCCGAAGGTGCTGGCATTGTCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCgctgct                                         <  1:903301/150‑1 (MQ=255)
                                                                                                             ggTTCCCGAAGGTGCTGGCATTGACACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCtgt                                       >  2:195517/1‑150 (MQ=255)
                                                                                                              gTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCtgtg                                      >  2:39590/1‑150 (MQ=255)
                                                                                                               ttCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCTgtgt                                     >  1:538178/1‑150 (MQ=255)
                                                                                                                   cGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCTGTGtt                                    <  2:851683/147‑1 (MQ=255)
                                                                                                                   cGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCTGTGtt                                    >  1:851683/1‑147 (MQ=255)
                                                                                                                        gTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAAc                                                                                        <  2:760908/90‑1 (MQ=255)
                                                                                                                        gTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAAc                                                                                        >  1:760908/1‑90 (MQ=255)
                                                                                                                           cTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCAcc                         <  2:1018874/150‑1 (MQ=255)
                                                                                                                            tGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCtgt                                       >  2:224604/1‑135 (MQ=255)
                                                                                                                            tGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCtgt                                       <  1:224604/135‑1 (MQ=255)
                                                                                                                             ggCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGc                                                     >  1:449138/1‑120 (MQ=255)
                                                                                                                             ggCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGc                                                     <  2:449138/120‑1 (MQ=255)
                                                                                                                                aTTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTgcg                    >  2:654211/1‑150 (MQ=255)
                                                                                                                                    gCACTGATGGTCATGATGACCATGGGGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCAGCACCCTGCGCAAg                >  1:56950/1‑150 (MQ=255)
                                                                                                                                    gCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAg                >  2:63771/1‑150 (MQ=255)
                                                                                                                                    gCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAg                <  2:9233/150‑1 (MQ=255)
                                                                                                                                     gctccaTGGTCATCATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCtgt                                       >  2:580598/6‑126 (MQ=255)
                                                                                                                                     gctccaTGGTCATCATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCtgt                                       <  1:580598/121‑1 (MQ=255)
                                                                                                                                     gctccaTGGTCATCATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGc               >  2:158318/6‑150 (MQ=255)
                                                                                                                                      aCTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCtggtggtg                                                                   >  1:873201/1‑97 (MQ=255)
                                                                                                                                      aCTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCtggtggtg                                                                   <  2:873201/97‑1 (MQ=255)
                                                                                                                                      aCTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCa              >  1:1106201/1‑150 (MQ=255)
                                                                                                                                       tccaTGGTCATCATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACggg                                                                                     >  2:215255/4‑78 (MQ=255)
                                                                                                                                       tccaTGGTCATCATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACggg                                                                                     <  1:215255/75‑1 (MQ=255)
                                                                                                                                       cTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCAt             <  1:507848/150‑1 (MQ=255)
                                                                                                                                          aTGGTCATCATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAg                                                   >  2:1187263/1‑109 (MQ=255)
                                                                                                                                          aTGGTCATCATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAg                                                   <  1:1187263/109‑1 (MQ=255)
                                                                                                                                             gTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCTGCGTTGGTGGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTAACGCCGCACCCTGCGCAAGCATGACCTg       >  1:866825/1‑150 (MQ=255)
                                                                                                                                             gTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGACCTg       <  1:106041/150‑1 (MQ=255)
                                                                                                                                               cATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCCAAGGGCGGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGACCTGGa     <  1:677471/150‑1 (MQ=255)
                                                                                                                                                  gatgaCCATGTAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGAGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGACCTGGAgcg  >  2:868417/1‑150 (MQ=255)
                                                                                                                                                  gatgaCCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGACCTGGAgcg  >  1:976255/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GTGGTGTTCAGTGCCTTCATCACAGGCCTGGGCCGCGCCATGCTGATGCTGCAGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTGGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGACCTGGAGCG  >  NC_002947/821536‑821831

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: