Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,670,458 T→A 5.8% intergenic (+47/‑38) metF → / → PP_4978 5,10‑methylenetetrahydrofolate reductase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,670,4580TA 5.8% 258.8 / 3.2 84intergenic (+47/‑38)metF/PP_49785,10‑methylenetetrahydrofolate reductase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (42/37);  new base A (2/3);  total (44/40)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGCAAGGCGGCGCGCCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCAGCATCGTCAGCGCCCTCGC  >  NC_002947/5670319‑5670605
                                                                                                                                           |                                                                                                                                                   
ttGCAAGGCGGCGCGCCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTaaa                                                                                                                                             <  1:321327/148‑1 (MQ=255)
ttGCAAGGCGGCGCGCCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAg                                                                                                                                            >  4:245865/1‑149 (MQ=255)
   cAAGGCGGCGCGCCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAATTGATTAAGGCTTTGGTCATAGACTAAAGCCt                                                                                                                                         >  4:289420/1‑149 (MQ=255)
   cAAGGCGGCGCGCCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCt                                                                                                                                         >  4:289419/1‑149 (MQ=255)
     aGGCGGCGCGCCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTACGGCTTTGGTCATAGACTAAAGCCTTa                                                                                                                                       <  5:72817/149‑1 (MQ=255)
      ggcggcGCGCCAGGCCTGTACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTAt                                                                                                                                      <  3:179152/149‑1 (MQ=255)
       gcggcgCGCCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAGATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTAGTCTATGACCAAAGCCTTaat                                                                                                                                     <  2:206022/149‑3 (MQ=39)
         ggcgCGCCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCa                                                                                                                                   >  5:290562/1‑149 (MQ=255)
          gcgcgcCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCAt                                                                                                                                  <  5:202410/149‑1 (MQ=255)
              gcCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTc                                                                                                                              <  1:153227/149‑1 (MQ=255)
               ccAGGCCTGCACTTCTATACCTTGAACCATGCTGAACCGAGCCTGGCGATCTTGAACAACCTGAAGCTGCCGCTCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCg                                                                                                                             <  5:201097/149‑1 (MQ=255)
                cAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCgg                                                                                                                            <  3:289419/149‑1 (MQ=255)
                cAGGCCTGCACTTCTATACCTTGAACCAGGATGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCgg                                                                                                                            <  3:289420/149‑1 (MQ=255)
                        cACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTaaa                                                                                                                      >  1:19571/1‑147 (MQ=255)
                         aCTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAg                                                                                                                   >  6:236577/1‑149 (MQ=255)
                           ttCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGt                                                                                                                 >  8:106296/1‑149 (MQ=255)
                             cTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCt                                                                                                               >  5:349006/1‑149 (MQ=255)
                             cTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCt                                                                                                               >  7:293218/1‑149 (MQ=255)
                                taCCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTtgc                                                                                                            <  7:106296/149‑1 (MQ=255)
                                        aCCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTAGTCTATGACCAAAGCCTTAATCAt                                                                                                                                  <  3:117707/119‑1 (MQ=38)
                                        aCCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTAGTCTATGACCAAAGCCTTAATCAt                                                                                                                                  >  4:117707/1‑119 (MQ=38)
                                          cAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATgca                                                                                                  <  3:245865/149‑1 (MQ=255)
                                          cAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATg                                                                                                    >  4:343320/1‑147 (MQ=255)
                                             gCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATgcagca                                                                                               >  4:13109/1‑149 (MQ=255)
                                                    cGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGAATAAGGCTTTAGTCTATGACCAAAGCCTTAATCAtttcttaa                                                                                                                           >  2:299092/1‑109 (MQ=21)
                                                    cGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGAATAAGGCTTTAGTCTATGACCAAAGCCTTAATCAtttcttaa                                                                                                                           <  1:299092/114‑6 (MQ=21)
                                                     gAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAAc                                                                                       <  4:122098/149‑1 (MQ=255)
                                                       gCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCt                                                                                     <  8:293218/149‑1 (MQ=255)
                                                        ccTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCtt                                                                                    <  3:13109/149‑1 (MQ=255)
                                                            gCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTcgcg                                                                                <  2:19571/149‑1 (MQ=255)
                                                            gCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTcgcg                                                                                >  7:150647/1‑149 (MQ=255)
                                                                tCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCca                                                                            >  3:143743/1‑149 (MQ=255)
                                                                   ggAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAg                                                                         >  7:68702/1‑149 (MQ=255)
                                                                   ggAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAg                                                                         >  1:71666/1‑149 (MQ=255)
                                                                    gAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGc                                                                        >  3:78052/1‑149 (MQ=255)
                                                                     aacaacCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCACCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATGTTCGGGTTAAACAGGTCTAGCCTCCCATGAAGCTCACCCCACCTTCGCGCCCACGGCt                                                                       >  8:129554/1‑149 (MQ=255)
                                                                      acaacCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCtt                                                                      <  6:349006/149‑1 (MQ=255)
                                                                          ccTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCAttt                                                                                                                                >  1:122954/1‑87 (MQ=255)
                                                                          ccTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCAttt                                                                                                                                <  2:122954/87‑1 (MQ=255)
                                                                          ccTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGta                                                                   >  8:36069/1‑148 (MQ=255)
                                                                           cTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGtata                                                                 >  8:101633/1‑149 (MQ=255)
                                                                             gAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGTAGCACACCCAACCTTCGCGCCCACAGCTTGTATAcc                                                               >  2:25001/1‑149 (MQ=255)
                                                                             gAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATAcc                                                               >  2:280402/1‑149 (MQ=255)
                                                                              aaGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATAc                                                                <  7:62258/147‑1 (MQ=255)
                                                                                    ccGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCt                                                          <  6:290562/147‑1 (MQ=255)
                                                                                      gcgcTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATAAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTg                                                      >  4:264479/1‑149 (MQ=255)
                                                                                         cTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGcc                                                   <  4:78052/149‑1 (MQ=255)
                                                                                          tGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCATCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGccc                                                  >  5:352776/1‑149 (MQ=255)
                                                                                          tGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGccc                                                  >  6:176527/1‑149 (MQ=255)
                                                                                            aaaTCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCgg                                                >  7:233927/1‑149 (MQ=255)
                                                                                            aaaTCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCg                                                 <  8:68702/148‑1 (MQ=255)
                                                                                            aaaTCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCg                                                 <  3:343320/148‑1 (MQ=255)
                                                                                               tCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCaa                                             >  1:313039/1‑149 (MQ=255)
                                                                                                   ttGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAAc                                            <  1:280402/146‑1 (MQ=255)
                                                                                                   ttGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAAc                                            <  5:176527/146‑1 (MQ=255)
                                                                                                      gCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAAc                                                                                                                     >  3:303268/1‑70 (MQ=255)
                                                                                                      gCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAAc                                                                                                                     <  4:303268/70‑1 (MQ=255)
                                                                                                      gCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCa                                      >  7:188173/1‑149 (MQ=255)
                                                                                                        aaaGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCAt                                     >  1:61290/1‑148 (MQ=255)
                                                                                                        aaaGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTAATTTTCGGGTTAATCAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACACCTTGTATACCTGGCTTTTGGCCCGGCAACCTACAAt                                     >  1:325987/1‑148 (MQ=255)
                                                                                                         aaGTCTGCTTAAGAAATGATAAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCa                                   <  3:264479/149‑1 (MQ=255)
                                                                                                          aGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGTCCCGGCAACCTACCATCAg                                  <  8:233927/149‑1 (MQ=255)
                                                                                                          aGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAg                                  >  6:232069/1‑149 (MQ=255)
                                                                                                           gTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCtt                                                                      <  1:103816/112‑1 (MQ=255)
                                                                                                           gTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCtt                                                                      >  2:103816/1‑112 (MQ=255)
                                                                                                                cTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGaa                               <  4:143743/146‑1 (MQ=255)
                                                                                                                cTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGcc                            >  5:117839/1‑149 (MQ=255)
                                                                                                                cTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGcc                            >  5:92248/1‑149 (MQ=255)
                                                                                                                cTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGc                             <  5:232069/148‑1 (MQ=255)
                                                                                                                cTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGc                             <  8:150647/148‑1 (MQ=255)
                                                                                                                   aaGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCACGGCAACCTACCATCAGGAAGCCTGc                         >  4:348219/1‑149 (MQ=255)
                                                                                                                    aGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGc                         <  7:101633/148‑1 (MQ=255)
                                                                                                                       aaTGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCa                     <  6:92248/149‑1 (MQ=255)
                                                                                                                       aaTGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCa                     <  7:36069/149‑1 (MQ=255)
                                                                                                                        aTGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTcgcg                                                                                >  8:128797/1‑89 (MQ=255)
                                                                                                                        aTGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTcgcg                                                                                <  7:128797/89‑1 (MQ=255)
                                                                                                                        aTGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCAg                    <  3:348219/149‑1 (MQ=255)
                                                                                                                            ttAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCAGCATc                >  5:129176/1‑149 (MQ=255)
                                                                                                                                   tttGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCa                     >  6:321134/1‑137 (MQ=255)
                                                                                                                                   tttGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCa                     <  5:321134/137‑1 (MQ=255)
                                                                                                                                      ggTCTTAGACTAAAGCCTTATTCATTTTCGGGTTAATCAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGtata                                                                 <  6:179839/90‑1 (MQ=255)
                                                                                                                                      ggTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGACTTGCCTGCCAAGCAGCACACCCAACCTTCGCGCCCACAGCTTGtata                                                                 >  5:179839/1‑90 (MQ=255)
                                                                                                                                         cATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCAGCATCGTCAGCGCCCTCg   >  1:57596/1‑149 (MQ=255)
                                                                                                                                          aTAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCAGCATCGTCAGCGCCCTCGc  >  1:209127/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                   
TTGCAAGGCGGCGCGCCAGGCCTGCACTTCTATACCTTGAACCAGGCTGAACCGAGCCTGGCGATCTGGAACAACCTGAAGCTGCCGCGCTGAAATCTTTTGGCAAAGTCTGCTTAAGAAATGATTAAGGCTTTGGTCATAGACTAAAGCCTTATTCATTTTCGGGTTAAACAGGTCTTGCCTGCCATGCAGCACACCCAACCTTCGCGCCCACAGCTTGTATACCTGGCTTTTGGCCCGGCAACCTACCATCAGGAAGCCTGCTTCAGCATCGTCAGCGCCCTCGC  >  NC_002947/5670319‑5670605

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: