Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,050,543 C→G 8.2% P450R (CCG→CGG)  spoT → bifunctional (p)ppGpp synthase/hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,050,5430CG 8.2% 222.7 / 10.6 73P450R (CCG→CGG) spoTbifunctional (p)ppGpp synthase/hydrolase
Reads supporting (aligned to +/- strand):  ref base C (30/37);  new base G (2/4);  total (32/41)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.76e-01

GCCCAAAGGCTCCACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCCTGCTGAACAAGGTACTCACT  >  NC_002947/6050397‑6050682
                                                                                                                                                  |                                                                                                                                           
gCCCAAAGGCTCCACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCggg                                                                                                                                          >  2:977120/1‑150 (MQ=255)
      aGGCTCCACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCGGGGCGCtga                                                                                                                                    <  2:831595/150‑3 (MQ=255)
       ggCTCCACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCgg                                                                                                                                   <  1:1040799/150‑1 (MQ=255)
        gCTCCACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCGGGGCGCtgacg                                                                                                                                  <  1:1137156/150‑5 (MQ=255)
         cTCCACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGcc                                                                                                                                 <  1:128404/150‑1 (MQ=255)
           ccACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCaa                                                                                                                               >  2:1377436/1‑150 (MQ=255)
            cACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCaaa                                                                                                                              >  2:248510/1‑150 (MQ=255)
                gCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCg                                                                                                                          <  2:1114395/150‑1 (MQ=255)
                   gTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCa                                                                                                                       <  1:1030690/150‑1 (MQ=255)
                      gACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATgg                                                                                                                    <  2:104419/150‑1 (MQ=255)
                        cTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGcg                                                                                                                  >  2:1016878/1‑149 (MQ=255)
                         ttCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTc                                                                                                                 >  1:1058333/1‑150 (MQ=255)
                            gCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAAc                                                                                                              <  2:829879/150‑1 (MQ=255)
                               tACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACttt                                                                                                           >  1:116972/1‑150 (MQ=255)
                                 cgcgGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAGGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTtgt                                                                                                         <  2:217703/150‑1 (MQ=255)
                                 cgcgGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTtgt                                                                                                         <  2:252910/150‑1 (MQ=255)
                                 cgcgGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTtgt                                                                                                         <  1:1377436/150‑1 (MQ=255)
                                 cgcgGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTtgt                                                                                                         <  2:1204698/150‑1 (MQ=255)
                                  gcgGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATTGCTCAACTTtgtg                                                                                                        <  1:409706/150‑1 (MQ=255)
                                  gcgGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTtgtg                                                                                                        <  1:1353587/150‑1 (MQ=255)
                                  gcgGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTtgtg                                                                                                        <  2:1384565/150‑1 (MQ=255)
                                  gcgGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTtgtg                                                                                                        <  1:640768/150‑1 (MQ=255)
                                        cacaCCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGcccc                                                                                                                                             >  1:209869/1‑107 (MQ=255)
                                        cacaCCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGcccc                                                                                                                                             <  2:209869/107‑1 (MQ=255)
                                        cacaCCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCg                                                                                                  >  2:1111809/1‑150 (MQ=255)
                                         acacCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCgg                                                                                                 <  1:453273/150‑1 (MQ=255)
                                         acacCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCgg                                                                                                 >  1:200344/1‑150 (MQ=255)
                                          cacCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCggg                                                                                                <  1:977120/150‑1 (MQ=255)
                                           acCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCgg                                                                                                                         <  1:185060/124‑1 (MQ=255)
                                           acCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCgg                                                                                                                         >  2:185060/1‑124 (MQ=255)
                                           acCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCGCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCAGGGCTCAACGTTGTGGTCGCgggg                                                                                               >  1:864207/1‑149 (MQ=255)
                                                cGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCgg                                                                                                 <  2:1231177/143‑1 (MQ=255)
                                                cGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCgg                                                                                                 >  1:1231177/1‑143 (MQ=255)
                                                  tCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGgcac                                                                                        >  2:1425257/1‑150 (MQ=255)
                                                     gCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCGGGGCGCtga                                                                                                                                    <  1:846086/103‑3 (MQ=255)
                                                     gCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCGGGGCGCtga                                                                                                                                    >  2:846086/1‑101 (MQ=255)
                                                      cAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCgg                                                                                                 <  1:912086/137‑1 (MQ=255)
                                                      cAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCgg                                                                                                 >  2:912086/1‑137 (MQ=255)
                                                      cAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGgcacgcac                                                                                    <  2:116972/150‑1 (MQ=255)
                                                            ttGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTc                                                                                                                                     <  1:205722/95‑1 (MQ=255)
                                                            ttGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTc                                                                                                                                     >  2:205722/1‑95 (MQ=255)
                                                              gCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATcc                                                                            <  2:1273565/150‑1 (MQ=255)
                                                               cATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCaa                                                                                                               <  1:965131/114‑1 (MQ=255)
                                                               cATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCaa                                                                                                               >  2:965131/1‑114 (MQ=255)
                                                                aTCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGTAAGGCACGCACGAATATCCGc                                                                          >  1:949481/1‑150 (MQ=255)
                                                                aTCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGc                                                                          <  2:200344/150‑1 (MQ=255)
                                                                   gCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCaa                                                                                                               >  2:778524/1‑110 (MQ=255)
                                                                   gCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCaa                                                                                                               <  1:778524/110‑1 (MQ=255)
                                                                    cTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCgg                                                                                                                         <  2:1257361/99‑1 (MQ=255)
                                                                    cTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCgg                                                                                                                         >  1:1257361/1‑99 (MQ=255)
                                                                         cGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTc                                                                 >  1:802550/1‑150 (MQ=255)
                                                                         cGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTc                                                                 >  1:322839/1‑150 (MQ=255)
                                                                          gCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTtgtg                                                                                                        <  2:439953/110‑1 (MQ=255)
                                                                          gCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTtgtg                                                                                                        >  1:439953/1‑110 (MQ=255)
                                                                                aCCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCGGGGCGCtga                                                                                                                                    <  2:585350/76‑3 (MQ=255)
                                                                                aCCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCGGGGCGCtga                                                                                                                                    >  1:585350/1‑74 (MQ=255)
                                                                                aCCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTGTGTGTTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAAc                                                          <  2:1058333/150‑1 (MQ=255)
                                                                                 ccgccgCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACa                                                         >  1:126986/1‑150 (MQ=255)
                                                                                  cgccgcCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAACCCCGGCATGGCTCAACATTGTGGGCTCGGTCTATGCACGCACGAATATCCGCCACGCGCTCAAGCAACAg                                                        <  1:248510/150‑1 (MQ=255)
                                                                                   gccgccTGGCGCCGCTGTCCGAACCGCTACAAAGCGACTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAgg                                                                                           >  2:712530/1‑114 (MQ=255)
                                                                                   gccgccTGGCGCCGCTGTCCGAACCGCTACAAAGCGACTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAgg                                                                                           <  1:712530/114‑1 (MQ=255)
                                                                                             gccgcTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGa                                             <  2:126986/150‑1 (MQ=255)
                                                                                               cgcTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGcgccgc                                                  <  2:175020/143‑1 (MQ=255)
                                                                                               cgcTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGcgccgc                                                  >  1:175020/1‑143 (MQ=255)
                                                                                                       gAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCgg                                                                                                                         <  1:436809/64‑1 (MQ=255)
                                                                                                       gAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCgg                                                                                                                         >  2:436809/1‑64 (MQ=255)
                                                                                                                   aGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAgcgc                       >  1:689064/1‑150 (MQ=255)
                                                                                                                            aCAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGAGCGGCTCGGAGTCCATCAGCCTGGGCGAGCGCCTGCTGAAc              >  1:785772/1‑150 (MQ=255)
                                                                                                                                tGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCCTGCTGAACAAgg          >  2:1386325/1‑150 (MQ=255)
                                                                                                                                gggAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCCTGCTGAACAAgg          <  2:1093883/149‑1 (MQ=255)
                                                                                                                                  gAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCCTGCTGAACAAGGTa        >  2:1177965/1‑150 (MQ=255)
                                                                                                                                    aaTCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCaa                                                                                             <  2:337594/63‑1 (MQ=255)
                                                                                                                                    aaTCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCaa                                                                                             >  1:337594/1‑63 (MQ=255)
                                                                                                                                       cGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTCTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCCTGCTGAACAAGGTACTCAc   <  1:321162/150‑1 (MQ=255)
                                                                                                                                        gTCCGCGCCCCGGGCGCTCGGCGTACCCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGGACGTCCCCATATCCGCGACGCGCTCACCCAACCGCGCCGCTCGGAGTCCATCAGCCTGGGCGTGTGCCTGCGGAACAAGGTCCTCACt  <  2:864207/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                           
GCCCAAAGGCTCCACGGCCGTCGACTTCGCCTACGCGGTCCACACCGACGTCGGCAACAGTTGCATCGCTTGCCGCATCAACCGCCGCCTGGCGCCGCTGTCCGAACCGCTACAAAGCGGCTCGACAGTGGAAATCGTCAGCGCCCCGGGCGCTCGGCCAAACCCGGCATGGCTCAACTTTGTGGTCTCGGGCAAGGCACGCACGAATATCCGCCACGCGCTCAAGCAACAGCGCCGCTCGGAGTCCATCAGCCTGGGCGAGCGCCTGCTGAACAAGGTACTCACT  >  NC_002947/6050397‑6050682

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