Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 142,885 A→G 7.2% intergenic (+7/‑168) PP_0134 → / → PP_5428 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947142,8850AG 7.2% 238.2 / 4.3 82intergenic (+7/‑168)PP_0134/PP_5428hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (35/41);  new base G (4/2);  total (39/43)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCAGGCATGGTGCGGCTCTCGGGCGAGGTGGTCAGCGCGGAGCAGAAAACCATCGCCGTGGAGATCGCCCGGCAAATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTGGCATGCAGGCTGCAATTACCTGTTCAGGTTTGAACAAAAAA  >  NC_002947/142739‑143033
                                                                                                                                                    |                                                                                                                                                    
aaGCAGGCATGGTGCGGCTCTCGGGCGAGGTGGTCAGCGCGGAGCAGAAAACCATCGCCGTGGAGATCGCCCGGCAAATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAAt                                                                                                                                                    >  7:374226/1‑149 (MQ=255)
        aTGGTGCGGCTCTCGGGCGAGGTGGTCAGCGCGGAGCAGAAAACCATCGCCGTGGAGATCGCCCGGCAAATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGa                                                                                                                                            <  2:271170/149‑1 (MQ=255)
               ggCTCTCGGGCGAGGTGGTCAGCGCGGAGCGAGAAAACCATCGCCGTGGAGATCGCCCGGCAAATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACcaca                                                                                                                                       >  1:244931/1‑148 (MQ=255)
               ggCTCTCGGGCGAGGTGGTCAGCGCGGAGCAGAAAACCATCGCCGTGGAGATCGCCCGGCAAATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACcaca                                                                                                                                       >  3:163047/1‑147 (MQ=255)
                gCTCTCGGGCGAGGTGGTCAGCGCGGAGCAGAAAACCATCGCCGTGGAGATCGCCCGGCAAATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAAt                                                                                                                                                    <  1:252071/133‑1 (MQ=255)
                gCTCTCGGGCGAGGTGGTCAGCGCGGAGCAGAAAACCATCGCCGTGGAGATCGCCCGGCAAATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAAt                                                                                                                                                    >  2:252071/1‑133 (MQ=255)
                gCTCTCGGGCGAGGTGGTCAGCGCGGAGCAGAAAACCATCGCCGTGGAGATCGCCCGGCAAATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAg                                                                                                                                    >  8:146135/1‑149 (MQ=255)
                   ctcGGGCGAGGTGGTCAGTGCGGAGCAGAAAACCATCGCCGTGGAGATCGCCCGGCATTATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGc                                                                                                                                  >  1:116972/1‑149 (MQ=255)
                    tcGGGCGAGGTGGTCAGCGCGGAGCAGAAAACCATCGCCGTGGAGATCGCCCGGCAAATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGc                                                                                                                                  <  2:244931/147‑1 (MQ=255)
                        gCGAGGTGGTCAGCGCGGAGCAGAAAACCATCGCCGTGGAGATCGCCCGGCAAATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATc                                                                                                                            >  6:31456/1‑149 (MQ=255)
                           aGGTGGTCAGCGCGGAGCAGAAAACCATCGCCGTGGAGATCGCCCGGCAAATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTg                                                                                                                         <  4:380547/149‑1 (MQ=255)
                             gtggtCAGCGCGGAGCAGAAAACCATCGCCGTGGAGATCGCCCGGCAAATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTg                                                                                                                         >  6:395392/1‑147 (MQ=255)
                              tggtCAGCGCGGAGCAGAAAACCATCGCCGTGGAGATCGCCCGGCAAATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCGAATTGCACGATcag                                                                                                                                        >  7:201678/1‑130 (MQ=255)
                              tggtCAGCGCGGAGCAGAAAACCATCGCCGTGGAGATCGCCCGGCAAATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCGAATTGCACGATcag                                                                                                                                        <  8:201678/131‑2 (MQ=255)
                                gtCAGCGCGGAGCAGAAAACCATCGCCGTGGAGATCGCCCGGCAAATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAAt                                                                                                                                                    >  8:226047/1‑117 (MQ=255)
                                gtCAGCGCGGAGCAGAAAACCATCGCCGTGGAGATCGCCCGGCAAATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAAt                                                                                                                                                    <  7:226047/117‑1 (MQ=255)
                                gtCAGCGCGGAGCAGAAAACCATCGCCGTGGAGATCGCCCGGCAAATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTg                                                                                                                         <  7:43698/144‑1 (MQ=255)
                                gtCAGCGCGGAGCAGAAAACCATCGCCGTGGAGATCGCCCGGCAAATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAgg                                                                                                                     <  5:40469/148‑1 (MQ=255)
                                gtCAGCGCGGAGCAGAAAACCATCGCCGTGGAGATCGCCCGGCAAATCATCGGCGTACGGGCCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTg                                                                                                                         >  8:43698/1‑144 (MQ=255)
                                 tCAGTGCGGAGCAGAAAACCATCGCCGTGGAGATCGCCCGGCAAATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGAt                                                                                                                   <  2:116972/149‑1 (MQ=255)
                                       cgGAGCAGAAAACCATCGCCGTGGAGATCGCCCGGCAAATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCAc                                                                                                                                              <  4:179033/116‑1 (MQ=255)
                                       cgGAGCAGAAAACCATCGCCGTGGAGATCGCCCGGCAAATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCAc                                                                                                                                              >  3:179033/1‑116 (MQ=255)
                                            cAGAAAACCATCGCCGTGGAGATCGCCCGGCAAATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTcaca                                                                                                        <  5:31456/149‑1 (MQ=255)
                                            cAGAAAACCATCGCCGTGGAGATCGCCCGGCAAATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTcaca                                                                                                        <  7:146135/149‑1 (MQ=255)
                                                    ccATCGCCGTGGAGATCGCCCGGCAAATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGtttt                                                                                                 >  1:104749/1‑149 (MQ=255)
                                                       tCGCCGTGGAGATCGCCCGGCAAATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTta                                                                                              <  1:249146/149‑1 (MQ=255)
                                                          ccGTGGAGATCGCCCGGCAAATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTtaata                                                                                           <  1:205356/149‑1 (MQ=255)
                                                           cGTGGAGATCGCCCGGCAAATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTtaataa                                                                                          <  8:374226/149‑1 (MQ=255)
                                                              ggAGATCGCCCGGCAAATCATCGGCGTACGCGGCGTGGAAGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCCTCCAGGATACGTCCTCAAAGAGTTTTTTAATAATCa                                                                                       >  2:363528/1‑149 (MQ=255)
                                                               gagaTCGCCCGGCAAATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCac                                                                                      <  4:371977/149‑1 (MQ=255)
                                                                  aTCGCCCGGCAAATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGt                                                                                                                                >  3:305888/1‑104 (MQ=255)
                                                                  aTCGCCCGGCAAATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGt                                                                                                                                <  4:305888/104‑1 (MQ=255)
                                                                      cccGGCAAATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTGTTTAATAATCACTAACTAg                                                                               >  6:242677/1‑149 (MQ=255)
                                                                         ggCAAATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGtt                                                                             >  6:367893/1‑148 (MQ=255)
                                                                         ggCAAATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTg                                                                            <  8:60914/149‑1 (MQ=255)
                                                                           cAAATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTGATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCAGTAACTAGTTGAt                                                                          <  4:117545/149‑1 (MQ=255)
                                                                            aaaTCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTcaca                                                                                                        >  7:285880/1‑118 (MQ=255)
                                                                            aaaTCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTcaca                                                                                                        <  8:285880/118‑1 (MQ=255)
                                                                            aaaTCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTc                                                                                                           <  1:28796/115‑1 (MQ=255)
                                                                            aaaTCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTc                                                                                                           >  2:28796/1‑115 (MQ=255)
                                                                            aaaTCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGAtt                                                                         >  6:174690/1‑149 (MQ=255)
                                                                            aaaTCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGAtt                                                                         <  1:113765/149‑1 (MQ=255)
                                                                                 catcGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTtat                                                                      >  1:156690/1‑148 (MQ=255)
                                                                                  atcGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCAGACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATaa                                                                    <  2:278295/149‑1 (MQ=255)
                                                                                      gCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGAt                                                                <  5:367893/149‑1 (MQ=255)
                                                                                           cGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGTATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATtttat                                                           <  1:245367/149‑1 (MQ=255)
                                                                                            ggggCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATtttatt                                                          <  5:174690/149‑1 (MQ=255)
                                                                                              ggCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATtttattt                                                         >  6:144062/1‑148 (MQ=255)
                                                                                                 gTGGATGCTGACCTGCTGAAAGTCGCGACCATGGTGGAAGGCTGATGGTGCAATTGGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCAGAgtgt                                                                                                    <  5:330582/102‑1 (MQ=255)
                                                                                                 gTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACACGCGTCATCGTGCAGGATACGTCCTCACAgtgt                                                                                                    >  6:330582/1‑102 (MQ=255)
                                                                                                  tGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATc                                                    <  5:395392/149‑1 (MQ=255)
                                                                                                      tGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCaaaa                                                <  4:197290/149‑1 (MQ=255)
                                                                                                       gCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGttt                                                                                                  <  5:178025/98‑1 (MQ=255)
                                                                                                       gCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGttt                                                                                                  >  6:178025/1‑98 (MQ=255)
                                                                                                           aactgcttAAAGTCGCGACCATGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGTTTCGTCCTCACAGTGTTTTTTAATAATGACTAACTAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTg                                           <  1:363528/147‑1 (MQ=255)
                                                                                                           aCCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTg                                           >  5:401867/1‑149 (MQ=255)
                                                                                                            ccTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCaaa                                                 <  1:309139/142‑1 (MQ=255)
                                                                                                            ccTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCaaa                                                 >  2:309139/1‑142 (MQ=255)
                                                                                                             ctgctgAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTgg                                          >  5:14628/1‑148 (MQ=255)
                                                                                                             ctgctgAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTGGc                                         <  5:144062/149‑1 (MQ=255)
                                                                                                               gctgAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTGGCAt                                       >  3:99584/1‑149 (MQ=255)
                                                                                                                  gAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTGGCATga                                     >  3:243669/1‑147 (MQ=255)
                                                                                                                  gAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTGGCATGc                                     <  4:243669/148‑1 (MQ=255)
                                                                                                                      gTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGa                                                                           >  4:246463/1‑106 (MQ=255)
                                                                                                                      gTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTGGCATGCAGGCt                                <  6:64949/149‑1 (MQ=255)
                                                                                                                           gACCAAGGTTGAAGGCTGATCGTGCAATTCGCACGAACACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGAt                                                                          <  7:142947/102‑1 (MQ=255)
                                                                                                                           gACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATAGTGCAGGATACGTCCTCACAGTCTTTTTTAATAAACACTAACTAGTAGAt                                                                          >  8:142947/1‑102 (MQ=255)
                                                                                                                             ccAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTGGCATGCAGGCTGCAATTa                         >  5:45986/1‑149 (MQ=255)
                                                                                                                              cAAGGTGGAAGGCTGATCGTGCGATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAAt                                                     >  7:261051/1‑120 (MQ=255)
                                                                                                                              cAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGAtttt                                                             <  8:116844/112‑1 (MQ=255)
                                                                                                                              cAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGAtttt                                                             >  7:116844/1‑112 (MQ=255)
                                                                                                                              cAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAAt                                                     <  8:261051/120‑1 (MQ=255)
                                                                                                                                aGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCaaa                                                 >  3:371552/1‑122 (MQ=255)
                                                                                                                                aGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCAACGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCaaa                                                 <  4:371552/122‑1 (MQ=255)
                                                                                                                                  gTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTGGCATGCAGGCTGCAATTACCTGt                    <  7:8481/149‑1 (MQ=255)
                                                                                                                                            gATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGGTACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTGGCATGCAGGCTGCAATTACCTGTTCAGGTTTg           >  3:183304/1‑148 (MQ=255)
                                                                                                                                             gtCGTGCGAATTGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTGGCATGCAGGCTGCAATTACCTGTTCAGGTTTGa          >  8:38679/2‑148 (MQ=255)
                                                                                                                                             aTCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTGGCATGCAGGCTGCAATTACCTGTTCAGGTTTGa          <  2:104749/148‑1 (MQ=255)
                                                                                                                                              tCGTGCGAATTGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTGGCATGCAGGCTg                               >  8:123297/1‑126 (MQ=255)
                                                                                                                                              tCGTGCGAATTGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTGGCATGCAGGCTg                               <  7:123297/126‑1 (MQ=255)
                                                                                                                                                   cAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTgg                                          <  1:340022/110‑1 (MQ=255)
                                                                                                                                                   cAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAAATAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTgg                                          >  2:340022/1‑110 (MQ=255)
                                                                                                                                                    aaTTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTGGCATGCAGGCTGCAATTACCTGTTCAGGTTTGAACaaaaaa  <  6:401867/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AAGCAGGCATGGTGCGGCTCTCGGGCGAGGTGGTCAGCGCGGAGCAGAAAACCATCGCCGTGGAGATCGCCCGGCAAATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTGGCATGCAGGCTGCAATTACCTGTTCAGGTTTGAACAAAAAA  >  NC_002947/142739‑143033

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: