Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,878,022 A→G 5.8% intergenic (‑37/+54) PP_5151 ← / ← PP_5152 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,878,0220AG 5.8% 252.6 / 3.0 85intergenic (‑37/+54)PP_5151/PP_5152hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (41/39);  new base G (3/2);  total (44/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGATTGAACCAGCCACGCTCGAACAGGCGAAACTCCTGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTCTTTTTGAATCGATAGAAAAGCGCCGGCTCGCTGACCATGTACAAGTTGCCCCTGTCGTCCATGGCCACCCCTTCAGCCCGCGGAATG  >  NC_002947/5877879‑5878169
                                                                                                                                               |                                                                                                                                                     
agaTTGAACCAGCCACGCTCGAACAGGCGAAACTCCTGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGaa                                                                                                                                                   <  7:118623/148‑1 (MQ=255)
  aTTGAACCAGCCACGCTCGAACAGGCGAAACTCCTGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATgg                                                                                                                                               <  5:274346/149‑1 (MQ=255)
  aTTGAACCAGCCACGCTCGAACAGGCGAAACTCCTGCCCAGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAAtc                                                                                                                                                 >  4:67329/1‑147 (MQ=255)
   ttGAACCAGCCACGCTCGAACAGGCGAAACTCCTGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATgg                                                                                                                                               <  4:303544/148‑1 (MQ=255)
   ttGAACCAGCCACGCTCGAACAGGCGAAACTCCTGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATgg                                                                                                                                               <  4:116924/148‑1 (MQ=255)
      aaCCAGCCACGCTCGAACAGGCGAAACTCCTGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTg                                                                                                                                            >  6:398653/1‑148 (MQ=255)
       aCCAGCCACGCTCGAACAGGCGAAACTCCTGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGaa                                                                                                                                          <  1:315686/149‑1 (MQ=255)
        ccAGCCACGCTCGAACAGGCGAAACTCCTGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTg                                                                                                                                            >  8:197704/1‑146 (MQ=255)
         cAGCCACGCTCGAACAGGCGAAACTCCTGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAg                                                                                                                                         >  2:312466/1‑148 (MQ=255)
              aCGCTCGAACAGGCGAAACTCCTGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGc                                                                                                                                   >  2:321148/1‑149 (MQ=255)
                gCTCGAACAGGCGAAACTCCTGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCt                                                                                                                                  >  8:248532/1‑148 (MQ=255)
                gCTCGAACAGGCGAAACTCCTGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCt                                                                                                                                  >  8:11589/1‑148 (MQ=255)
                     aaCAGGCGAAACTCCTGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATg                                                                                                                             >  1:347062/1‑148 (MQ=255)
                        aGGCGAAACTCCTGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAg                                                                                                                         >  3:146094/1‑149 (MQ=255)
                          gCGAAACTCCTGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCg                                                                                                                       >  5:93677/1‑149 (MQ=255)
                            gAAACTCCTGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTc                                                                                                                    >  5:377008/1‑149 (MQ=255)
                               aCTCCTGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAgg                                                                                                                 >  5:19389/1‑149 (MQ=255)
                                cTCCTGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGAAGCATAGGGACAGACTCGACACAAGGATGGCGATTCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAgg                                                                                                                 >  3:227012/1‑148 (MQ=255)
                                 tCCTGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGc                                                                                                                        >  3:320901/1‑141 (MQ=255)
                                 tCCTGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGc                                                                                                                        <  4:320901/141‑1 (MQ=255)
                                 tCCTGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGc                                                                                                                >  8:40944/1‑148 (MQ=255)
                                   cTGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGCCGCATAGGGACAGACTCGACTCAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAAt                                                                                                             >  7:55424/1‑149 (MQ=255)
                                   cTGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAAt                                                                                                             >  2:221077/1‑149 (MQ=255)
                                    tGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAAtt                                                                                                            >  4:388164/1‑149 (MQ=255)
                                    tGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAAt                                                                                                             >  2:322827/1‑148 (MQ=255)
                                     gCCCGGCCCAGCCATGCAGCAGCAGAGCAAGGGCTGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTaa                                                                                                                               <  1:373027/130‑1 (MQ=255)
                                     gCCCGGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTaa                                                                                                                               >  2:373027/1‑130 (MQ=255)
                                      cccGGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAAtt                                                                                                            <  5:350829/147‑1 (MQ=255)
                                      cccGGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTc                                                                                                           <  6:377008/148‑1 (MQ=255)
                                       ccGGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCtg                                                                                                         <  3:388164/149‑1 (MQ=255)
                                         ggCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCtgat                                                                                                       <  7:11589/149‑1 (MQ=255)
                                         ggCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCtgat                                                                                                       <  1:221077/149‑1 (MQ=255)
                                           cccAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTaa                                                                                                                               <  5:370331/124‑1 (MQ=255)
                                           cccAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTaa                                                                                                                               >  6:370331/1‑124 (MQ=255)
                                           cccAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCtgatga                                                                                                     >  2:298416/1‑149 (MQ=255)
                                                 cAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGAtgttgt                                                                                               <  7:248532/149‑1 (MQ=255)
                                                 cAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGAtgttgt                                                                                               <  1:362481/149‑1 (MQ=255)
                                                 cAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGAtgttgt                                                                                               <  1:312466/149‑1 (MQ=255)
                                                    gcagcagcaGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGCCGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTAAg                                                                                            >  1:368979/1‑149 (MQ=255)
                                                    gcagcagcaGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAg                                                                                            <  5:398653/149‑1 (MQ=255)
                                                         agcagAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGCCAGACTCGACACAAGGATGGCGAATCTAATCGTGAATGGCTAAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTCtt                                                                                       >  1:330420/1‑149 (MQ=255)
                                                         agcagAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCGTGAATGGCTAAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTCtt                                                                                       >  1:334990/1‑149 (MQ=255)
                                                          gcagAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCtgat                                                                                                       <  3:43072/132‑1 (MQ=255)
                                                          gcagAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCtgat                                                                                                       >  4:43072/1‑132 (MQ=255)
                                                            agagAAATGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGTCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTCttttt                                                                                    <  4:227012/149‑1 (MQ=255)
                                                             gagCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATCGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTCttttt                                                                                    <  1:322827/149‑1 (MQ=255)
                                                             gagCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTCttttt                                                                                    >  6:303218/1‑148 (MQ=255)
                                                                cAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAATGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTCTTTTTGAAt                                                                                <  2:347062/149‑1 (MQ=255)
                                                                 aaGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTCTTTTTGAAt                                                                                <  5:303218/148‑1 (MQ=255)
                                                                 aaGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTCTTTTTGAATc                                                                               >  8:279992/1‑149 (MQ=255)
                                                                   gggCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTCTTTTTGAATCGa                                                                             <  1:298416/149‑1 (MQ=255)
                                                                   gggCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTCTTTTTGAATCGa                                                                             >  1:83631/1‑149 (MQ=255)
                                                                    ggCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTCTTTTTGAATCGAt                                                                            <  6:19389/149‑1 (MQ=255)
                                                                           aGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTCTTTTTGAATCGATAGAAAAg                                                                     >  2:309493/1‑149 (MQ=255)
                                                                             aaCAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTCttttt                                                                                    <  1:204848/132‑1 (MQ=255)
                                                                             aaCAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTCttttt                                                                                    >  2:204848/1‑132 (MQ=255)
                                                                             aaCAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTCTTTTTGAATCGATAGAAAAgcg                                                                   >  5:378857/1‑149 (MQ=255)
                                                                             aaCAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTATGATGATGTTGTCAGTTCTTTTAGAATCGATAGAAAAgcg                                                                   >  5:372645/1‑149 (MQ=255)
                                                                              aCAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTCTTTTTGAATCGATAGAAAAgcgc                                                                  <  2:18763/149‑1 (MQ=255)
                                                                                 aGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTCTTTTTGAATCGATAGAAAAGCGCCgg                                                               <  1:321148/149‑1 (MQ=255)
                                                                                        ggCCAACAAAACGACGCATAGGGACAGACTCGACCCAAGGATCGCGAATCTAATCATCAATGCCTCAAGCAAAGCTTAATCCAAGCGTCAGTCAATTCTGATGATGTTGTCAGTTCTTTTTGAATCGATAGAAAAGCGCCGGCTCGCTg                                                        <  8:55424/149‑1 (MQ=255)
                                                                                         gCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTCTTTTTGAATCGATAGAAACGCGCCGGCTCGCTGa                                                       <  7:279992/149‑1 (MQ=255)
                                                                                            aaCAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTCTTTTTGAATCGATAGAAAAGCGCCGGCTCGCTGACCa                                                    <  6:378857/149‑1 (MQ=255)
                                                                                            aaCAAAACGACGCATAGGGACAGACTAGACACAAGTATGGAGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTCTTTTTGAATCGATAGAAAAGCGCCGGCTCGCTGACCa                                                    <  6:372645/149‑1 (MQ=255)
                                                                                               aaaaCGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTCTTTTTGAATCGa                                                                             <  5:427240/121‑1 (MQ=255)
                                                                                               aaaaCGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTCTTTTTGAATCGa                                                                             >  6:427240/1‑121 (MQ=255)
                                                                                                 aaCGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTCTTTTTGAATCGATAGAAAAGCGCCGGCTCGCTGACCATGTa                                                >  4:347826/1‑148 (MQ=255)
                                                                                                   cgacgCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTCTTTTTGAATCGATAGAAAAGCGCCGGCTCGCTGACCATGTACaa                                             <  3:67329/149‑1 (MQ=255)
                                                                                                     aagCATAGGTACAGAATCTACACAAGGATGGAGAATCTAATCATGAATGTCTTAAGCAAAGCTTAATGCAAGCGTCAGTCAATTCTGATGATTTTGTCAGTTCTTTTTGAATCGATAGAAAAGCGCAGGCTCGCTGACAATGTACAAGt                                           <  7:340063/147‑1 (MQ=255)
                                                                                                      cgCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTCTTTTTGAATCGATAGAAAAGCGCCGGCTCGCTGACCATGTACAAGtt                                          <  7:197704/149‑1 (MQ=255)
                                                                                                       gCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTCTTTTTGAATCGATAGAAAAGCGCCGGCTCGCTGACCATGTACAAGTTg                                         >  1:238095/1‑149 (MQ=255)
                                                                                                         aTAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTCTTTTTGAATCGATAGAAAAGCGCCGGCTCGc                                                          <  2:194085/130‑1 (MQ=255)
                                                                                                         aTAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTCTTTTTGAATCGATAGAAAAGCGCCGGCTCGCTGACCATGTACAAGTTGcc                                       >  1:2280/1‑149 (MQ=255)
                                                                                                              gACAGACTCGACACAAGGATGGCGAATCTAATCGTGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTCTTTTTGAATCGATAGAAAAGCGCCGGCTCGCTGACCATGTACAAGTTGCCCCTgt                                  >  7:330350/1‑149 (MQ=255)
                                                                                                               aCAGACTCGACACAAGGATTGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTc                                                                                         <  7:360583/93‑1 (MQ=255)
                                                                                                               aCAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTc                                                                                         >  8:360583/1‑93 (MQ=255)
                                                                                                               aCAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTCTTTTTGAATCGATAGAAAAGCGCCGGCTCGCTGACCATGTACAAGTTGCCCCTgtc                                 <  4:146094/149‑1 (MQ=255)
                                                                                                                       gACACAAGGATGGCGAATCTAATCGTGAATGGCTAAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTCTTTTTGAATCGATAGAAAAGCGCCGGCTCGCTGACCATGTACAAGTTGCCCCTGTCGTCCATg                          <  2:334990/148‑1 (MQ=255)
                                                                                                                       gACACAAGGATGGCGAATCTAATCGTGAATGGCTAAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTCTTTTTGAATCGATAGAAAAGCGCCGGCTCGCTGACCATGTACAAGTTGCCCCTGTCGTCCATg                          <  2:330420/148‑1 (MQ=255)
                                                                                                                       gACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTCTTTTTGAATCGATAGAAAAGCGCCGGCTCGCTGACCATGTACAAGTTGCCCCTGTCGTCCATg                          <  3:347826/148‑1 (MQ=255)
                                                                                                                        acacAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTCTTTTTGAATCGATAGAAAAGCGCCGGCTCGCTGACCATGTACAAGTTGCCCCTGTCGTCCATGGc                        >  2:376416/1‑149 (MQ=255)
                                                                                                                          acaAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTCTTTTTGAATCGATAGAAAAGCGCCGGCTCGCTGACCATGTACAAGTTGCCCCTGTCGTCCATGGCCa                      >  1:133527/1‑149 (MQ=255)
                                                                                                                                  ggCGAATCTAAACATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTGCAGGCAATTCTGATGATGTTGTCAGTTCTTTTTGAATCGATAGAAAAGCGCCGGCTCGCTGACCATGTACAAGTTGCCCCTGTCGTCCATGGCCACCCCtt                >  1:90425/1‑148 (MQ=255)
                                                                                                                                   gCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTCTTTTTGAATCGATAGAAAAGCGCCGGCTCGCTGACCATGTACAAGTTGCCCCTGTCGTCCATGGCCACCCCTTCAg             >  6:203907/1‑149 (MQ=255)
                                                                                                                                     gAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTCTTTTTGAATCGATAGAAAAGCGCCGGCTCGCTGACCATGTACAAGTTGCCCCTGTCGTCCATGGCCACCCCTTCAGcc           >  1:64158/1‑149 (MQ=255)
                                                                                                                                               aTGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTCTTTTTGAATCGATAGAAAAGCGCCGGCTCGCTGACCATGTACAAGTTGCCCCTGTCGTCCATGGCCACCCCTTCAGCCCGCGGAATg  <  1:376416/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
AGATTGAACCAGCCACGCTCGAACAGGCGAAACTCCTGCCCGGCCCAGCCAAGCAGCAGCAGAGCAAGGGCAGCCAGAACAAGGAACAGGCGAACAAAACGACGCATAGGGACAGACTCGACACAAGGATGGCGAATCTAATCATGAATGGCTGAAGCAAAGCTTAATGCAAGCGTCAGGCAATTCTGATGATGTTGTCAGTTCTTTTTGAATCGATAGAAAAGCGCCGGCTCGCTGACCATGTACAAGTTGCCCCTGTCGTCCATGGCCACCCCTTCAGCCCGCGGAATG  >  NC_002947/5877879‑5878169

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