Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 445,110 T→G 18.9% intergenic (+80/‑208) PP_0367 → / → PP_0368 hypothetical protein/acyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947445,1100TG18.9% 410.0 / 80.6 178intergenic (+80/‑208)PP_0367/PP_0368hypothetical protein/acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (76/68);  new base G (8/26);  total (84/94)
Fisher's exact test for biased strand distribution p-value = 2.21e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

Only 100 of 201 total aligned reads displayed.

CGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTT  >  NC_002947/444961‑445259
                                                                                                                                                       |                                                                                                                                                      
cGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtt                                                                                                                                                        <  2:4803577/150‑1 (MQ=255)
cGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtt                                                                                                                                                        <  1:9126814/150‑1 (MQ=255)
  aGCTCGGCGCACAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACTCCTACTCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGTCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                                                      <  2:6943357/150‑1 (MQ=255)
  aGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                                                      <  2:8428280/150‑1 (MQ=255)
  aGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                                                      <  2:6496386/150‑1 (MQ=255)
  aGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                                                      <  1:6378340/150‑1 (MQ=255)
  aGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                                                      <  1:5632262/150‑1 (MQ=255)
  aGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                                                      <  1:5238762/150‑1 (MQ=255)
   gCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGGGTGGCCGCGCAAGCGTCCCCTCAGACACCCAAtttt                                                                                                                                                      >  2:5093121/1‑150 (MQ=255)
     tCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTc                                                                                                                                                  <  1:3579079/150‑1 (MQ=255)
       ggCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTg                                                                                                                                                <  2:3847882/150‑1 (MQ=255)
       ggCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCGTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTg                                                                                                                                                >  2:7537879/1‑150 (MQ=255)
        gcgcAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTgg                                                                                                                                               <  2:8566935/150‑1 (MQ=255)
        gcgcAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAAAATTGTTGTCTga                                                                                                                                               <  2:2390159/150‑2 (MQ=255)
        gcgcAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTga                                                                                                                                               <  2:386023/150‑2 (MQ=25)
          gcAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGACCAGACAAAATTGTTGTCTGATc                                                                                                                                             <  1:4890969/150‑1 (MQ=38)
          gcAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATc                                                                                                                                             <  2:6678256/150‑1 (MQ=25)
           cAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGACCAGACAAAATTGTTGTCTGATCg                                                                                                                                            <  2:9139516/150‑1 (MQ=38)
             aaGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCggg                                                                                                                                          <  2:5695189/150‑1 (MQ=255)
              aGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGc                                                                                                                                         <  1:7717578/150‑1 (MQ=255)
              aGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGc                                                                                                                                         <  1:8444266/150‑1 (MQ=255)
              aGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGACCAGACAAAATTGTTGTCTGATCggac                                                                                                                                         <  1:9028645/150‑3 (MQ=38)
              aGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGCCCGACCAGACAAAATTGTTGTCTGATCggac                                                                                                                                         <  1:784197/150‑3 (MQ=37)
                   gtcAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCtt                                                                                                                                    >  2:8354663/1‑150 (MQ=255)
                      aaGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTgg                                                                                                                                               <  2:6700098/136‑1 (MQ=255)
                        ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCTCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTcgcgt                                                                                                                               <  2:2562934/150‑2 (MQ=255)
                        ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                                               <  2:3465918/150‑1 (MQ=255)
                        ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                               <  1:9099582/150‑1 (MQ=18)
                        ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                               <  1:7031319/150‑1 (MQ=18)
                        ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                               <  1:4546444/150‑1 (MQ=18)
                        ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                               <  1:2370291/150‑1 (MQ=18)
                        ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                               <  2:935717/150‑1 (MQ=18)
                         gCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCggg                                                                                                                              <  2:4886023/150‑1 (MQ=255)
                           aTCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCa                                                                                                                            <  2:7224138/150‑1 (MQ=255)
                              cgacgaAATGCTTGGCCGGATCTTCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGc                                                                                                                         <  1:41504/150‑1 (MQ=255)
                              cgacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGc                                                                                                                         <  1:3619271/150‑1 (MQ=17)
                                  gaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGGATGCGTCCCCTCACCCAACAAGTTTGTCTGGGGGGGGCTCGTCGCGGGCAAGCCCGc                                                                                                                     >  2:8065725/1‑150 (MQ=255)
                                          tGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                                                      >  2:1511749/1‑110 (MQ=255)
                                          tGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                                                      <  1:1511749/110‑1 (MQ=255)
                                          tGGCCGGACCATCGACACCTACGCCTGAATCCTGGCTTCATCCCCTATTTCCGCGCCGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGACCGATCAGACAACAATTTTGTGTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACaa                                                                                                             <  2:7632886/150‑1 (MQ=37)
                                                   catcGACACCTCCGCCTTAATACTTGCGACATCCCCTATTTCAGCGCCTGCATTGTAGCTGCGGGCTTGCCCGCGAATGCGACCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTtgtg                                                                                                    <  2:8205146/150‑1 (MQ=38)
                                                     tcGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCGTGCCCCCGCAAGCGGCCCCTCAGACAACAAAGTTTGCGGGGCGGGCCGCGTCGCGGGCCACCCCCCGCCAACAACGGGGTGtggc                                                                                                  >  2:3748259/1‑150 (MQ=255)
                                                       gACACCACCGCGTTAATAGTTGGGACAGCCCCTATTCCCGCGCGGGCTTTGTAGGCGGGGGGTGGCCGGCTATGGGGCCCGATCACACAAAAATTTGGTGGGGTGGGGCCTCTGCGGGGCCACCCCCGCTCCTACAAGGGTTTGtggcat                                                                                                <  2:8730419/150‑1 (MQ=255)
                                                        acacCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGc                                                                                                                                         >  1:286698/1‑108 (MQ=255)
                                                        acacCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGc                                                                                                                                         <  2:286698/108‑1 (MQ=255)
                                                                    gAATACTTGCGACAGCCCCTTTTTCCGCGCCGGCCTTGTAGGTGGGGGCTTGCCCGCGACTGCGTCCTAACAGACAACAATTTGGTCTGGTCGGGCCTGTTCGCGGGCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGAtt                                                                                   <  2:6640801/150‑1 (MQ=255)
                                                                     atgcggggCGTCACCCCCTATGCCCGGCCCGGCGTTGTCGGCGGGGGGTGCCCCGCGACGGCGCCCTCTCACACAACACTTGTGTGTGGCGGGGCCGCGGGGGGGGCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGAttt                                                                                  <  2:7900112/143‑1 (MQ=255)
                                                                                    cccTATTTCAGCGCCTGCCTTGTAGGAGCGGGGGTGCCCCCCAATTCGGGCGCTCAGACAAAAAATTTTTGTGGTCGGGGCGCGTCCCGGGCCATCCCCCACCAACAATGGTTTTTGGCATGGCATGAAAGATTTTCCAATTTTGGTGGc                                                                   >  2:1414537/1‑150 (MQ=255)
                                                                                     ccTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCCAAATCGTCCCATCAGTCAACAAATTTGTCTGGGAGGGCCGCGTGGCGGGCAAGCCCCCTCCGACAAGGGTGTTTGGCAAGGCATGAATGATTTGGCAATTTTAGTGGCg                                                                  >  2:4554678/1‑150 (MQ=255)
                                                                                     ccTATTTCAGCGCCTGCCTTGTAGGAGCGGGCGTGCCCCCGAATGCGGACGCTCCCACAACAAATTTTTGTGGGGGGGCCTCGTCGCGGGCCAGCCCCCTCCCACAAAGGGGTGTGGGCTGGCATGTCTGAATTGCCAATTTTCCTGGCg                                                                  >  2:3649191/1‑150 (MQ=255)
                                                                                          ttCCGCGCCGGCCTTGTCGGCGGGGGCTGCCCCGCCACTGCGCCCGAACACACAAAATTTTGTTCGGGTCGGGCCTCTCCGGGGCCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAAt                                                             <  2:3296616/150‑1 (MQ=255)
                                                                                                   gcctttgttgGCGCGGGCGGCCCCGCGACTGCGTCCCAACCAACACCATTTTGGGGTGGTGGGGCCTCCCCGCGGGCACCCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATa                                                    <  1:7537879/146‑1 (MQ=255)
                                                                                                           tggAGCGGGCTTGCGCTCTAAGGCTTCCTCTCCTCAAACATTTTGGTCGGGGGGGGCCTCTTGGCGGGAATGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATcc                                            <  2:1968839/149‑1 (MQ=255)
                                                                                                             gcgggggCTGGCCCGCCACTGCGCCCCACCACACAACAATTGTGTGGGGGCGGCCCTCTCCGCGGGCACCCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCtt                                          <  1:6742231/146‑1 (MQ=255)
                                                                                                              gcgggggTTGCCCGGGTATTGCGTCCTATCACACAAAATTTTTGTGGGGCGGGGCCTCTGCGCGGCCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCtt                                          <  1:2986595/147‑1 (MQ=255)
                                                                                                              aGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTACCTCCTTg                                         >  2:5885032/1‑150 (MQ=255)
                                                                                                               gCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTTTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTTCCATTTTTCGTGGCGTAAATTGCACCATAGTCCATCATTGa                                        >  1:7998886/1‑150 (MQ=255)
                                                                                                                     ttGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGaaa                                  >  1:7281357/1‑150 (MQ=255)
                                                                                                                        cccGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                               >  1:5512363/1‑150 (MQ=18)
                                                                                                                        cccGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                               >  2:2994741/1‑150 (MQ=18)
                                                                                                                         ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGGGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                              >  1:4993116/1‑150 (MQ=255)
                                                                                                                         ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACaa                                                                                                             >  2:7185886/1‑71 (MQ=255)
                                                                                                                         ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGTCGTAAAt                                                             <  1:7172161/119‑1 (MQ=255)
                                                                                                                         ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGTCGTAAAt                                                             >  2:7172161/1‑119 (MQ=255)
                                                                                                                         ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                              >  2:3094776/1‑150 (MQ=255)
                                                                                                                         ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                              >  1:6503236/1‑150 (MQ=255)
                                                                                                                         ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                              >  1:1901012/1‑150 (MQ=255)
                                                                                                                         ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                              >  2:6175340/1‑150 (MQ=255)
                                                                                                                         ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                              >  1:444634/1‑150 (MQ=255)
                                                                                                                                    gTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTaa                        >  2:1601660/1‑146 (MQ=255)
                                                                                                                                    gTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGt                    >  2:400243/1‑150 (MQ=255)
                                                                                                                                     tCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTa                   >  1:7493975/1‑150 (MQ=255)
                                                                                                                                      ccGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGa                             >  2:1263837/1‑139 (MQ=255)
                                                                                                                                      ccGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGa                             <  1:1263837/139‑1 (MQ=255)
                                                                                                                                      ccGACCAGACAAAATTGTTGTCTGATCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGAttt                                                                                  <  1:7055673/86‑1 (MQ=18)
                                                                                                                                      ccGACCAGACAAAATTGTTGTCTGATCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGAttt                                                                                  >  2:7055818/1‑86 (MQ=18)
                                                                                                                                       cGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAg                 >  2:1011676/1‑150 (MQ=255)
                                                                                                                                       cGACCAGACAAAATTGTTGTCTGATCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAg                 >  1:7225067/1‑150 (MQ=38)
                                                                                                                                         aTCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGtt               >  1:5121193/1‑150 (MQ=255)
                                                                                                                                          tCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGttt              >  2:4988145/1‑150 (MQ=255)
                                                                                                                                           cAGACAAAATTGTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTc             >  2:5219736/1‑150 (MQ=255)
                                                                                                                                           accacaacaTTTTGGTCTGGGCGCGCCACTCCGGGGCCACCCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTc             <  1:6807959/147‑1 (MQ=255)
                                                                                                                                             gACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaa           >  1:6678504/1‑150 (MQ=255)
                                                                                                                                             gACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaa           >  1:1677932/1‑150 (MQ=255)
                                                                                                                                                aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca        >  2:7286965/1‑150 (MQ=255)
                                                                                                                                                aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca        >  2:7732837/1‑150 (MQ=255)
                                                                                                                                                aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTAGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca        >  2:7671976/1‑150 (MQ=255)
                                                                                                                                                aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGAAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca        >  2:7672240/1‑150 (MQ=255)
                                                                                                                                                aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGAGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca        >  1:5587022/1‑150 (MQ=255)
                                                                                                                                                  acatTTTGGTGTGGTGGGGCCTCTGCGCGGGAACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaacaa       <  1:8139219/150‑1 (MQ=255)
                                                                                                                                                  acaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaacaa       >  1:808252/1‑150 (MQ=255)
                                                                                                                                                  aaatttttGTAGGGTCGGGCCTCCCCGGGGCCAACCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGATAGCTGAGCTAAACGTAAGTTTCAaacaa       <  1:3649191/146‑1 (MQ=255)
                                                                                                                                                   caaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaacaac      >  1:5812002/1‑150 (MQ=255)
                                                                                                                                                    aaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACt     >  2:8170980/1‑150 (MQ=255)
                                                                                                                                                    aaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACt     >  1:7083304/1‑150 (MQ=255)
                                                                                                                                                     aTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTg    >  2:9091468/1‑150 (MQ=255)
                                                                                                                                                      ttttGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGt   >  2:20738/1‑150 (MQ=255)
                                                                                                                                                      ttttGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACGGt   >  2:7853607/1‑150 (MQ=255)
                                                                                                                                                       tttGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGtt  <  2:6768158/150‑1 (MQ=255)
                                                                                                                                                       tttGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGtt  >  1:1407600/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                      
CGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTT  >  NC_002947/444961‑445259

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