Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,427,105 G→C 7.1% intergenic (‑197/‑603) PP_5647 ← / → PP_3917 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,427,1050GC 7.1% 668.0 / 30.0 197intergenic (‑197/‑603)PP_5647/PP_3917hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (94/89);  new base C (8/6);  total (102/95)
Fisher's exact test for biased strand distribution p-value = 7.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 200 total aligned reads displayed.

AGATATCCATTTCGCTTAAATGCGACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCTGCGCGCCTCG  >  NC_002947/4426956‑4427252
                                                                                                                                                     |                                                                                                                                                   
agataTCCATTTCGCTTAAATGCGACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAgcg                                                                                                                                                     >  2:211960/1‑150 (MQ=255)
   tatCCATTTCGCTTAAATGCGACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCgg                                                                                                                                                  <  1:1328280/150‑1 (MQ=255)
          ttCGCTTAAATGCGACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGgggctggg                                                                                                                                           <  1:4309762/150‑1 (MQ=255)
           tCGCTTAAATGCGACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGgggctgggc                                                                                                                                          <  1:2437991/150‑1 (MQ=255)
            cGCTTAAATGCGACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGgggctgggct                                                                                                                                         <  1:1188015/150‑1 (MQ=255)
             gCTTAAATGCGACCGGTGTAAAAAGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGACAGTAGAGAACCGAGCCCGATCACAGAGCCCGTCCAGCAGACAGAAGATCAAGTAGCAAAAAGCCCAGCGCGgggctgggcta                                                                                                                                        >  1:1331400/1‑149 (MQ=255)
             gCTTAAATGCGACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCtt                                                                                                                                        <  1:4443165/150‑1 (MQ=255)
             gCTTAAATGCGACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCtt                                                                                                                                        >  1:3968448/1‑150 (MQ=255)
               ttAAATGCGACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCTGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCt                                                                                                                                      <  2:462629/150‑1 (MQ=255)
               ttAAATGCGACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCt                                                                                                                                      <  1:876415/150‑1 (MQ=255)
               ttAAATGCGACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCt                                                                                                                                      <  1:3477887/150‑1 (MQ=255)
               ttAAATGCGACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCCCCGCGCTGGGCtttt                                                                                                                                      <  1:657293/150‑3 (MQ=255)
               ttAAATGCGACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCCCCGCGCTGGGCtttt                                                                                                                                      <  1:3745650/150‑3 (MQ=255)
                tAAATGCGACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCCCCGCGCTGGGCttttt                                                                                                                                     <  1:3273909/150‑4 (MQ=255)
                    tGCGACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTcag                                                                                                                                 <  1:2915105/150‑1 (MQ=255)
                      cGACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTcagca                                                                                                                               <  2:4198168/150‑1 (MQ=255)
                       gACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTcagcag                                                                                                                              >  2:2392932/1‑150 (MQ=255)
                       gACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTcagcag                                                                                                                              >  1:3165495/1‑150 (MQ=255)
                       gACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTcagcag                                                                                                                              >  2:2339696/1‑150 (MQ=255)
                       gACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTcagcag                                                                                                                              >  2:2042448/1‑150 (MQ=255)
                       gACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTcagcag                                                                                                                              >  2:1845876/1‑150 (MQ=255)
                       gACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTcagcag                                                                                                                              >  1:3715765/1‑150 (MQ=255)
                        aCCGGTGGAAAATGATCGGCCCTAATCCTAGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCACGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCTGGGCTGGGCTTCTGTCAGCaga                                                                                                                             <  2:905576/150‑1 (MQ=255)
                            gTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTc                                                                                                                         <  2:2827176/150‑1 (MQ=255)
                             tGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCg                                                                                                                        <  2:3349894/150‑1 (MQ=255)
                             tGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCg                                                                                                                        >  2:1989009/1‑150 (MQ=255)
                               gAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGAc                                                                                                                      <  1:3281686/150‑1 (MQ=255)
                               gAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGAc                                                                                                                      <  1:210410/150‑1 (MQ=255)
                               gAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGAc                                                                                                                      >  1:1913133/1‑150 (MQ=255)
                                         ggCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGc                                                                                                            >  2:652968/1‑150 (MQ=255)
                                          gCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCg                                                                                                           <  1:252717/150‑1 (MQ=255)
                                              tAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCa                                                                                                       <  2:3001433/150‑1 (MQ=255)
                                                      gCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCt                                                                                               <  1:1681619/150‑1 (MQ=255)
                                                         gatagataTCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCANACAGAAGTTCAAGTAGCAAAAAGCNCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTgg                                                                                            >  2:3719226/1‑150 (MQ=255)
                                                         gatagataTCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTgg                                                                                            >  1:3959786/1‑150 (MQ=255)
                                                          atagataTCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGc                                                                                           <  1:3054811/150‑1 (MQ=255)
                                                          atagataTCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGc                                                                                           <  1:2235169/150‑1 (MQ=255)
                                                           tagataTCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCa                                                                                          >  1:3426812/1‑150 (MQ=255)
                                                            agataTCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAg                                                                                         >  1:4269911/1‑150 (MQ=255)
                                                            agataTCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAg                                                                                         <  2:4300659/150‑1 (MQ=255)
                                                            agataTCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAg                                                                                         >  1:1418985/1‑150 (MQ=255)
                                                             gataTCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGc                                                                                        <  2:1683016/150‑1 (MQ=255)
                                                             gataTCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGc                                                                                        <  2:3385029/150‑1 (MQ=255)
                                                               tatCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTa                                                                                      <  2:96465/150‑1 (MQ=255)
                                                                 tCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACt                                                                                    >  1:1266394/1‑150 (MQ=255)
                                                                 tCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACt                                                                                    >  2:2809571/1‑150 (MQ=255)
                                                                    gCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCa                                                                                 >  1:4275505/1‑150 (MQ=255)
                                                                    gCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCa                                                                                 >  1:415368/1‑150 (MQ=255)
                                                                           aGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGgcg                                                                          >  1:871212/1‑150 (MQ=255)
                                                                                  gTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGa                                                                   <  1:1609137/150‑1 (MQ=255)
                                                                                   tAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGAt                                                                  >  1:512221/1‑150 (MQ=255)
                                                                                    agagAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATa                                                                 <  2:1969493/150‑1 (MQ=255)
                                                                                    agagAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATa                                                                 <  1:3685856/150‑1 (MQ=255)
                                                                                    agagAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCACAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATa                                                                 <  1:2678302/150‑1 (MQ=255)
                                                                                     gagaACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAg                                                                <  1:808137/150‑1 (MQ=255)
                                                                                      agaACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGc                                                               >  1:2472836/1‑150 (MQ=255)
                                                                                          ccGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCtata                                                           <  2:1264271/150‑1 (MQ=255)
                                                                                          ccGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCtata                                                           <  1:1354437/150‑1 (MQ=255)
                                                                                                 cGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCg                                                    >  1:321667/1‑150 (MQ=255)
                                                                                                 cGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCg                                                    <  2:2088251/150‑1 (MQ=255)
                                                                                                      aCTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCAtt                                               <  2:2570932/150‑1 (MQ=255)
                                                                                                      aCTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCAtt                                               <  1:1135163/150‑1 (MQ=255)
                                                                                                        tGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTgg                                             >  1:824018/1‑150 (MQ=255)
                                                                                                        tGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTGGTGATAGCTATAGCGCTCGGCATTgg                                             >  1:4007736/1‑150 (MQ=255)
                                                                                                          aGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGAt                                           <  1:218794/150‑1 (MQ=255)
                                                                                                          aGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGAt                                           <  1:3520623/150‑1 (MQ=255)
                                                                                                          aGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGAt                                           <  2:179938/150‑1 (MQ=255)
                                                                                                          aGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGAt                                           >  2:2171909/1‑150 (MQ=255)
                                                                                                          aGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGAt                                           >  1:2780469/1‑150 (MQ=255)
                                                                                                            cccGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGc                                         >  2:855069/1‑150 (MQ=255)
                                                                                                               gTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCt                                      >  2:940083/1‑150 (MQ=255)
                                                                                                                tCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCGTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTc                                     >  2:4367628/1‑150 (MQ=255)
                                                                                                                 ccAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCg                                    >  1:2381998/1‑150 (MQ=255)
                                                                                                                  cagcagACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGa                                   <  2:3667256/150‑1 (MQ=255)
                                                                                                                   agcagACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGAt                                  <  2:2548746/150‑1 (MQ=255)
                                                                                                                     cagacagaAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATgg                                <  2:3994993/150‑1 (MQ=255)
                                                                                                                     cagacagaAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATgg                                >  1:3871937/1‑150 (MQ=255)
                                                                                                                     cagacagaAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATgg                                <  1:3570697/150‑1 (MQ=255)
                                                                                                                          agaAGTTCAAGTAGCAAAAAGCCCAGCGCGCGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCg                           >  2:2927713/1‑150 (MQ=255)
                                                                                                                           gaAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCgg                          <  2:3548037/150‑1 (MQ=255)
                                                                                                                               ttCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTggg                      >  2:2657608/1‑150 (MQ=255)
                                                                                                                                 cAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGcc                    >  1:1994013/1‑150 (MQ=255)
                                                                                                                                  aaGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCt                                                                                       >  2:729703/1‑82 (MQ=255)
                                                                                                                                  aaGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCt                   >  1:6192/1‑150 (MQ=255)
                                                                                                                                   aGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCtt                  <  1:760610/150‑1 (MQ=255)
                                                                                                                                    gTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTg                 >  1:978535/1‑150 (MQ=255)
                                                                                                                                    gTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTg                 <  1:3790712/150‑1 (MQ=255)
                                                                                                                                       gCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGttt              >  1:689984/1‑150 (MQ=255)
                                                                                                                                         aaaaaGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCt            >  2:4412529/1‑150 (MQ=255)
                                                                                                                                          agaaGCCCAGCCCCGCGCTGGGCTTTTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCTg           >  2:4270264/3‑150 (MQ=255)
                                                                                                                                               cccAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCTgcgcgc      <  1:1989009/150‑1 (MQ=255)
                                                                                                                                               cccAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCTgcgcgc      <  2:1913133/150‑1 (MQ=255)
                                                                                                                                               cccAGCCCCGCGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCTgcgcgc      >  1:3328908/1‑150 (MQ=255)
                                                                                                                                                ccAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCTGCGCGcc     >  1:3224891/1‑150 (MQ=255)
                                                                                                                                                ccAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCTGCGCGcc     >  2:1412702/1‑150 (MQ=255)
                                                                                                                                                ccAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCTGCGCGcc     >  2:1072128/1‑150 (MQ=255)
                                                                                                                                                ccAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTGGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCTGCGCGcc     >  2:3391657/1‑150 (MQ=255)
                                                                                                                                                ccAGCCCCGCGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCTGCGCGcc     >  1:1095975/1‑150 (MQ=255)
                                                                                                                                                 cAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCTGCGCGCCt    >  1:1005216/1‑150 (MQ=255)
                                                                                                                                                  aGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCTGCGCGCCTc   <  2:3939187/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AGATATCCATTTCGCTTAAATGCGACCGGTGGAAAATGATCGGCCCTAATCCTTGCTGATAGATATCAGCGGGACAGCGCCAGTAGAGAACCGAGCCCGATCACTGAGCCCGTCCAGCAGACAGAAGTTCAAGTAGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTGTCAGCAGAGCTCGACTAGACTCAGCGTTCAATTACGCTTGGCAGCTACTCCAATGGGCGCTTGTGATAGCTATAGCGCTCGGCATTGGATGCCCTCGATGGGTCCGGTGGGCCTTGTTTCTGCGCGCCTCG  >  NC_002947/4426956‑4427252

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