Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 445,098 T→C 6.3% intergenic (+68/‑220) 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,0980TC 6.3% 226.0 / 2.6 78intergenic (+68/‑220)PP_0367/PP_0368hypothetical protein/acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (34/39);  new base C (3/2);  total (37/41)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

CAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTT  >  NC_002947/444949‑445247
                                                                                                                                                     |                                                                                                                                                      
cAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGAt                                                                                                                                                        >  2:2134690/1‑150 (MQ=255)
 ccgatcaaacagttgtttgaaacttacgtttagctcagctttctgtcaaggatggactatggtgccatttacgccacgaaaaatggcaaatcagtcatgccatgccacaaacCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATc                                                                                                                                                       >  1:1786030‑M1/113‑150 (MQ=255)
 aaGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGGCACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATc                                                                                                                                                       <  1:1823560/150‑1 (MQ=255)
  aGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCa                                                                                                                                                      <  2:1416338/150‑1 (MQ=255)
      cGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGaca                                                                                                                                                  >  2:2106247/1‑150 (MQ=255)
        aacagttgtttgaaacttacgtttagctcagctttctgtcaaggatggactatggtgccatttacgccacgaaaaatggcaaatcagtcatgccatgccacaaacCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGacaac                                                                                                                                                <  2:2066372‑M1/45‑1 (MQ=255)
           tCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGACATAGCGGGCTTGCCCGCGAATGCGTCCGACCAGacaaaatt                                                                                                                                             >  1:1204917/1‑146 (MQ=255)
            cGAGCTCGGCGCAAAGGTCGTCCAGGCATCCGACGAAATGCTTGGCCGGATCATCGACCCCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtt                                                                                                                                            <  1:117774/150‑1 (MQ=255)
              aGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                                         <  2:2272631/150‑1 (MQ=255)
                tttgaaacttacgtttagctcagctttctgtcaaggatggactatggtgccatttacgccacgaaaaatggcaaatcagtcatgccatgccacaaacCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGt                                                                                                                                       >  1:2066372‑M1/98‑150 (MQ=255)
                 tCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTc                                                                                                                                      <  1:2246658/150‑1 (MQ=255)
                 tCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTc                                                                                                                                      >  1:706098/1‑150 (MQ=37)
                      tcAAAGGTCGTAAAGGAATCAGTCGAACTGTTTTGCCTGATCACCCGCTCCTACTCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTTCGACCAGACAAAATTGTTGTCTGATc                                                                                                                                 <  2:1204917/149‑1 (MQ=19)
                      gcAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTc                                                                                                                                 <  1:1360660/150‑1 (MQ=255)
                        aaaGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCgg                                                                                                                               >  2:750179/1‑150 (MQ=255)
                        aaaGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCgg                                                                                                                               >  1:1376274/1‑150 (MQ=255)
                          aGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGc                                                                                                                             >  1:1555297/1‑150 (MQ=255)
                             tcgtcAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGCCCGACCAGACAAAATTGTTGTCTga                                                                                                                                   >  1:1792324/1‑140 (MQ=38)
                             tcgtcAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGCCCGACCAGACAAAATTGTTGTCTga                                                                                                                                   <  2:1792324/141‑2 (MQ=38)
                                tcAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTc                                                                                                                       <  2:1920825/150‑1 (MQ=255)
                                    ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                                   <  1:750179/150‑1 (MQ=255)
                                    ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAATTTGTAGGGTGGGGATTCTTCGCgg                                                                                                                   >  2:1959672/1‑150 (MQ=255)
                                    cagctttctgtcaaggatggactatggtgccatttacgccacggctaaatataaatcactcatgccatgccacaaacCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                                   >  2:664355‑M1/78‑150 (MQ=255)
                                    cagctttctgtcaaggatggactatggtgccatttacgccacgaaaaatggcaaatcagtcatgccatgccacaaacCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                                   >  2:1701452‑M1/78‑150 (MQ=255)
                                          cgacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACATCAATTTTGGCGGGTCCGGCCTCTTCGCgggtaggt                                                                                                             <  2:1273350/150‑6 (MQ=255)
                                          cgacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGc                                                                                                             <  1:2026551/150‑1 (MQ=255)
                                            acgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGc                                                                                                             >  2:256064/1‑148 (MQ=255)
                                            acgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGc                                                                                                             <  1:256064/148‑1 (MQ=255)
                                              tcaaggatggactatggtgccatttacgccacgaaaaatggcaaatcagtcatgccatgccacaaacCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTc                                                                                                                                 <  2:648840‑M1/59‑1 (MQ=255)
                                              tcaaggatggactatggtgccatttacgccacgaaaaatggcaaatcagtcatgccatgccacaaacCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTc                                                                                                                                 >  1:648840‑M1/68‑126 (MQ=255)
                                                 aTGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCTTGCCCGcgaa                                                                                                      <  2:2253025/150‑4 (MQ=255)
                                                        actagggtgccatttccgccaagaaaaagggcaactcagtcatgccatcccacaaccCCTTGTAGGAGCGGGCTGGCCCGCGAATGCGCCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAgg                                                                                               >  2:1781413‑M1/58‑150 (MQ=255)
                                                          cGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAAAATTGTTGTCTGATCgga                                                                                                                              <  2:1467805/117‑2 (MQ=37)
                                                          cGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAAAATTGTTGTCTGATCgga                                                                                                                              >  1:1467805/1‑116 (MQ=37)
                                                           ggATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGtt                                                                                            >  2:896295/1‑150 (MQ=255)
                                                                 tcGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAg                                                                                                >  1:1250015/1‑140 (MQ=255)
                                                                 tcGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAg                                                                                                <  2:1250015/140‑1 (MQ=255)
                                                                    acacCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGtggcatg                                                                                   <  1:896295/150‑1 (MQ=255)
                                                                    acacCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGtggcatg                                                                                   <  1:1959672/150‑1 (MQ=255)
                                                                      acCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGTTCAGACAACAATTTTGTCTCTACCTGGCGTGGCGCGGGCAAGCCCGCTCCTACAAGGGTTTGtggcatggc                                                                                 <  1:2254507/150‑1 (MQ=255)
                                                                        cTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGa                                                                                                                                                    >  1:1176733/1‑82 (MQ=255)
                                                                        cTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGa                                                                                                                                                    <  2:1176733/82‑1 (MQ=255)
                                                                        cTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTggcatggcat                                                                               >  1:2162242/1‑150 (MQ=255)
                                                                           tgCCTGAATACTTGCGTCCTCCCATATTTTGGCGCCGTCCTGGTAGCAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGAc                                                                            <  2:623437/149‑1 (MQ=38)
                                                                            gCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGtggcatggc                                                                                 >  1:133806/1‑144 (MQ=255)
                                                                            gCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGtggcatggc                                                                                 <  2:133806/144‑1 (MQ=255)
                                                                                gAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGtggcat                                                                                    >  2:1853428/1‑137 (MQ=255)
                                                                                gAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGtggcat                                                                                    <  1:1853428/137‑1 (MQ=255)
                                                                                 aaatggcaaatcagtcatgccatgccacaaacCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCaa                                                                                                               <  2:2054810‑M1/77‑1 (MQ=255)
                                                                                 aaatggcaaatcagtcatgccatgccacaaacCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCaa                                                                                                               >  1:2054810‑M1/33‑109 (MQ=255)
                                                                                              tCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACt                                                                           >  2:1854437/1‑132 (MQ=255)
                                                                                              tCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACt                                                                           <  1:1854437/132‑1 (MQ=255)
                                                                                                cccTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTTGc                                                       >  1:2416417/1‑150 (MQ=255)
                                                                                                cccTATTTCAGCGCCTGCATTGTAGGAGCGGGCTTGCCCGCGAATGGGTTCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCTTGGGATGACTGGTTTGCCATTTTTCGTGGc                                                       >  2:527940/1‑150 (MQ=38)
                                                                                                 ccTATTTCAGCGCCTGCCTTGTAGGAGCGGCCGCGCCCGCCAATGAGGTCGCTTGGTCTACCATTTTGTCTTGACGGGCCTCTTTACGGGCAACCCCCCTCCTTCAAGGGTTTTTGGGATCGCCTTTCTGATTTGCCAatttttttttcg                                                      >  2:1759525/1‑143 (MQ=255)
                                                                                                          gcgcCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCa                                             >  1:2320393/1‑150 (MQ=255)
                                                                                                          gcgcCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCa                                             <  2:2162242/150‑1 (MQ=255)
                                                                                                             cctgcctTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGtggcatg                                                                                   <  2:1695999/109‑1 (MQ=255)
                                                                                                             cctgcctTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGtggcatg                                                                                   >  1:1695999/1‑109 (MQ=255)
                                                                                                               tgcctTGTAGGAGCGGGCTTGCCCGCGACTGCGTCCGATCAGCCAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATa                                        <  1:1620541/150‑1 (MQ=255)
                                                                                                                gcctTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATaa                                       >  2:1981859/1‑149 (MQ=255)
                                                                                                                        GGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTNTCTGNTCGGGCCTCTTc                                                                                                                       <  1:345457‑M1/62‑1 (MQ=255)
                                                                                                                        GGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTc                                                                                                                       >  2:345457‑M1/1‑62 (MQ=255)
                                                                                                                             GGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACa                          <  2:2281003‑M1/150‑1 (MQ=255)
                                                                                                                              GGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTGTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAg                         <  1:642903‑M1/150‑1 (MQ=255)
                                                                                                                                  TGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAg                     >  1:922021‑M1/1‑150 (MQ=255)
                                                                                                                                     CCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGc                                                                                                                 >  1:1360800‑M1/1‑55 (MQ=255)
                                                                                                                                     CCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGc                                                                                                                 <  2:1360800‑M1/55‑1 (MQ=255)
                                                                                                                                     CCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCCCCATAGTCCATCCTTGACAGAAAGCTg                  <  2:288261‑M1/150‑1 (MQ=255)
                                                                                                                                     CCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                  <  1:1039513‑M1/150‑1 (MQ=255)
                                                                                                                                     CCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                  <  2:1617182‑M1/150‑1 (MQ=255)
                                                                                                                                     CCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                  <  1:2340438‑M1/150‑1 (MQ=255)
                                                                                                                                     CCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                  <  1:448337‑M1/150‑1 (MQ=255)
                                                                                                                                     CCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCCTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGa                                                                             >  2:2022639‑M1/1‑91 (MQ=255)
                                                                                                                                     CCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCCTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGa                                                                             <  1:2022639‑M1/91‑1 (MQ=255)
                                                                                                                                      CGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGa                 >  1:1769468‑M1/1‑150 (MQ=255)
                                                                                                                                       GCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAg                >  2:2340438‑M1/1‑150 (MQ=255)
                                                                                                                                             gCGTCGGCTCAGTCAACAATGTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGGTCCTGCACGGGTTTATGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAAc          <  2:722388/150‑1 (MQ=255)
                                                                                                                                               gTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGt        >  1:103352/1‑150 (MQ=255)
                                                                                                                                                     tCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGttt  >  1:1304398/1‑150 (MQ=255)
                                                                                                                                                     tCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGttt  >  1:772296/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTT  >  NC_002947/444949‑445247

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