Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 5,610,048 C→A S173I (AGC→ATC)  PP_4929 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,610,0480CA91.1% 228.9 / 13.6 90S173I (AGC→ATC) PP_4929LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (5/3);  new base A (43/39);  total (48/42)
Fisher's exact test for biased strand distribution p-value = 7.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.67e-01

CTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGC  >  NC_002947/5609903‑5610192
                                                                                                                                                 |                                                                                                                                                
cTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGcc                                                                                                                                              >  2:311788/1‑150 (MQ=255)
cTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGcc                                                                                                                                              >  2:263766/1‑150 (MQ=255)
 tGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGccc                                                                                                                                             >  1:1457945/1‑150 (MQ=255)
     ggCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACg                                                                                                                                         >  1:165910/1‑150 (MQ=255)
        tGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCAc                                                                                                                                      >  1:1602464/1‑150 (MQ=255)
              cctcgccGCCGGGGTAACAGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGgcgcc                                                                                                                                >  2:322980/1‑149 (MQ=255)
                      ccgGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACGGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCAAACGCACGGCGCTCCAAGGaa                                                                                                                        <  2:690186/150‑1 (MQ=255)
                          ggTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCCCGGCGCTCCAAGGAAGCCg                                                                                                                    >  2:755263/1‑150 (MQ=255)
                           gTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCCACTGTTGGCGGGTCACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCgg                                                                                                                   <  2:1357963/150‑1 (MQ=255)
                             aaCCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGa                                                                                                                 <  2:339772/150‑1 (MQ=255)
                              aCCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGAt                                                                                                                <  2:1457945/150‑1 (MQ=255)
                                cGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGCAGCCgttattc                                                                                                              >  2:1433119/1‑144 (MQ=255)
                                cGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGc                                                                                                              >  1:1357207/1‑150 (MQ=255)
                                  cTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCGATTGTTGGGGGGTGACCCCGCTTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGAACATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCtt                                                                                                            <  2:468477/150‑1 (MQ=255)
                                  cTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCtt                                                                                                            <  1:785251/150‑1 (MQ=255)
                                   tttGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCttt                                                                                                           >  2:535194/1‑150 (MQ=255)
                                      gttgtGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGGCCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGt                                                                                                        <  2:619456/150‑1 (MQ=255)
                                      gttgtGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGt                                                                                                        <  1:1284904/150‑1 (MQ=255)
                                      ggtgtgGGGTGATTAGCAACAGCCGGTGGCGGGCCAATTGTTGGCGGGGGACCCGGCCTTAGCCCGCGAGCGAGGGCACGACCGCGCACACCGTGACCATCTCCACGATGCCCACCGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGt                                                                                                        <  2:500939/148‑1 (MQ=255)
                                       ttgtGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGGGGGGGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTg                                                                                                       >  1:1634940/1‑150 (MQ=255)
                                         gtgGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTgga                                                                                                     <  1:1235489/150‑1 (MQ=255)
                                         gtgGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTgga                                                                                                     <  1:548044/150‑1 (MQ=255)
                                          tgGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTggag                                                                                                    >  1:831906/1‑150 (MQ=255)
                                          tgGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTggag                                                                                                    <  1:1162060/150‑1 (MQ=255)
                                              gTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCGACGCACGGCGCTGCAAGGAAGCCGGGATGCTTTCGTGGAGgaag                                                                                                <  2:1209164/150‑1 (MQ=255)
                                                  tCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACa                                                                                            >  2:1458776/1‑150 (MQ=255)
                                                  tCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACa                                                                                            >  1:843675/1‑150 (MQ=255)
                                                  tCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACa                                                                                            >  1:1251618/1‑150 (MQ=255)
                                                   cAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGcc                                                                                                                     >  2:44663/1‑124 (MQ=255)
                                                   cAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGcc                                                                                                                     <  1:44663/124‑1 (MQ=255)
                                                   cAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAg                                                                                           <  2:1602464/150‑1 (MQ=255)
                                                    aGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAgg                                                                                          <  1:1635007/150‑1 (MQ=255)
                                                     gCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGc                                                                                         >  2:252279/1‑150 (MQ=255)
                                                     gCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGc                                                                                         >  2:924357/1‑150 (MQ=255)
                                                       aaCTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGc                                                                                         >  1:511324/1‑148 (MQ=255)
                                                       aaCTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGc                                                                                         <  2:511324/148‑1 (MQ=255)
                                                          tGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCACCGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGa                                                                                    >  2:1097469/1‑150 (MQ=255)
                                                            ccGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGTGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCTCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATc                                                                                  <  1:486140/150‑1 (MQ=255)
                                                            ccGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCa                                                                                                                                       >  2:1266262/1‑97 (MQ=255)
                                                            ccGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCa                                                                                                                                       <  1:1266262/97‑1 (MQ=255)
                                                              ggTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCgg                                                                                <  2:1015380/150‑1 (MQ=255)
                                                               gTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGgc                                                                               <  1:924357/150‑1 (MQ=255)
                                                                   cGGGCCAATTGTTGGAGCGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGgcgcgc                                                                           >  1:60010/1‑150 (MQ=255)
                                                                        cAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTc                                                                      >  1:514955/1‑150 (MQ=255)
                                                                        cAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTc                                                                      >  2:210390/1‑150 (MQ=255)
                                                                        cAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAAATCGGCGCGCCGCTc                                                                      >  2:400690/1‑150 (MQ=255)
                                                                         aaTTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTcc                                                                     >  1:420068/1‑150 (MQ=255)
                                                                         aaTTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCAc                                                                                                                                      <  2:1173554/85‑1 (MQ=255)
                                                                         aaTTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCAc                                                                                                                                      >  1:1173554/1‑85 (MQ=255)
                                                                         aaTTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTcc                                                                     >  1:1003296/1‑150 (MQ=255)
                                                                         aaTTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTcc                                                                     >  2:1551528/1‑150 (MQ=255)
                                                                         aaTTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACAAAGTAACCATCTCCACAACTCCCAACGGAGGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTcc                                                                     >  1:211937/1‑150 (MQ=255)
                                                                                ggCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGACCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCt                                                              >  1:535645/1‑150 (MQ=255)
                                                                                 gCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTggag                                                                                                    >  2:719878/1‑111 (MQ=255)
                                                                                 gCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTggag                                                                                                    <  1:719878/111‑1 (MQ=255)
                                                                                     gggACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCCGCTTgcgc                                                         <  2:385585/148‑1 (MQ=255)
                                                                                            ggCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACAt                                                  >  2:393610/1‑150 (MQ=255)
                                                                                              ccTTCGCCCGCCCGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGTCGCTCCACCAGCTTGCGCGCCACATCg                                                <  1:977265/150‑1 (MQ=255)
                                                                                                ttCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGcc                                              >  1:341320/1‑150 (MQ=255)
                                                                                                     ccgccAGCGGGTGGCCCACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTggg                                         <  2:314568/150‑1 (MQ=255)
                                                                                                             gggTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTg                                 <  1:393610/150‑1 (MQ=255)
                                                                                                              ggTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTgg                                <  2:1251618/150‑1 (MQ=255)
                                                                                                               gTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGc                               <  2:843675/150‑1 (MQ=255)
                                                                                                                tGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGcc                              <  2:1256435/150‑1 (MQ=255)
                                                                                                                 ggCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCa                             <  2:165910/150‑1 (MQ=255)
                                                                                                                   ccGACCGCGCACACCGTAACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGt                           <  2:607886/150‑1 (MQ=255)
                                                                                                                     gACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTg                         >  2:683965/1‑150 (MQ=255)
                                                                                                                      aCCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTCCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGc                        >  1:552678/1‑150 (MQ=255)
                                                                                                                       ccGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGc                        >  1:49338/1‑149 (MQ=255)
                                                                                                                       ccGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGc                        <  2:49338/149‑1 (MQ=255)
                                                                                                                        cgcgCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGccc                                                                                       <  1:1004016/85‑1 (MQ=255)
                                                                                                                        cgcgCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGccc                                                                                       >  2:1004016/1‑85 (MQ=255)
                                                                                                                             acacCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCGGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCa                 >  1:942384/1‑150 (MQ=255)
                                                                                                                             acacCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCa                 <  1:1551528/150‑1 (MQ=255)
                                                                                                                             acacCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCa                 <  1:683965/150‑1 (MQ=255)
                                                                                                                                 cGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACCGGTCGGGATCGGCGCGCCGCTCCACCAGCTGGAGCGCCACATCGCTCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGaa             >  1:827270/1‑150 (MQ=255)
                                                                                                                                 cGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGcc                                              <  2:1338240/117‑1 (MQ=255)
                                                                                                                                 cGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGcc                                              >  1:1338240/1‑117 (MQ=255)
                                                                                                                                 cGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGaa             >  1:137925/1‑150 (MQ=255)
                                                                                                                                  gTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCGGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAAc            <  2:942384/150‑1 (MQ=255)
                                                                                                                                  gTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAAc            >  2:261028/1‑150 (MQ=255)
                                                                                                                                  gTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAAc            <  1:1311007/150‑1 (MQ=255)
                                                                                                                                      ccATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGt        <  1:688252/150‑1 (MQ=255)
                                                                                                                                       cATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTa       >  2:651659/1‑150 (MQ=255)
                                                                                                                                        aTCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCg                                                                                 <  2:788279/75‑1 (MQ=255)
                                                                                                                                        aTCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCg                                                                                 >  1:788279/1‑75 (MQ=255)
                                                                                                                                         tctcCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGc     >  2:320473/1‑150 (MQ=255)
                                                                                                                                           tcCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCg                                                                                 >  1:316201/1‑72 (MQ=255)
                                                                                                                                           tcCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCg                                                                                 <  2:316201/72‑1 (MQ=255)
                                                                                                                                            ccACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGc  >  1:885762/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CTGCGGGCTGATCGCCTCGCCGCCGGGGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGC  >  NC_002947/5609903‑5610192

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: