Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,151,969 T→A 7.4% R168S (AGA→AGT pphC ← protein‑serine/threonine phosphatase PphC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,151,9690TA 7.4% 242.4 / 7.5 81R168S (AGA→AGTpphCprotein‑serine/threonine phosphatase PphC
Reads supporting (aligned to +/- strand):  ref base T (42/33);  new base A (2/4);  total (44/37)
Fisher's exact test for biased strand distribution p-value = 4.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

TCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCGATG  >  NC_000913/2151820‑2152117
                                                                                                                                                     |                                                                                                                                                    
tCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGt                                                                                                                                                      >  1:1028430/1‑150 (MQ=255)
tCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGt                                                                                                                                                      >  2:1275101/1‑150 (MQ=255)
   tGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTg                                                                                                                                                   >  1:1310907/1‑150 (MQ=255)
     cGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGa                                                                                                                                                 >  2:97945/1‑150 (MQ=255)
        tgcttcCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCGAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAAc                                                                                                                                              <  2:1328157/150‑1 (MQ=255)
         gctgcCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACg                                                                                                                                             <  1:708427/150‑1 (MQ=255)
           tgcCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGc                                                                                                                                           >  1:1375125/1‑150 (MQ=255)
           tgcCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGACTGGAAActtt                                                                                                                                           <  1:1364285/150‑4 (MQ=255)
           tgcCGCCAGTCCATTGAAAAACGGCGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGACTGGAAActtt                                                                                                                                           <  1:1369432/150‑4 (MQ=255)
             ccgccAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGACTGGAAActttca                                                                                                                                         <  1:118/150‑6 (MQ=255)
              cgccAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGACTGGAAActttcac                                                                                                                                        <  1:727355/150‑7 (MQ=255)
               gccAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCt                                                                                                                                       <  2:252257/150‑1 (MQ=255)
               gccAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCt                                                                                                                                       >  1:549487/1‑150 (MQ=255)
                  aGTCCATTGAAAAACGGGGTAAAGAAAGGCCCATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTtcg                                                                                                                                    >  2:1016695/1‑150 (MQ=255)
                     ccATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTtcgtcg                                                                                                                                 <  1:126227/150‑1 (MQ=255)
                      cATTGAAAAACGGGGTAAAGAAAGGCACATGGGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCgg                                                                                                                                >  1:1128277/1‑150 (MQ=255)
                      cATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCgg                                                                                                                                >  1:955656/1‑150 (MQ=255)
                       aTTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGt                                                                                                                               >  2:1408812/1‑150 (MQ=255)
                       aTTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGt                                                                                                                               <  2:1071730/150‑1 (MQ=255)
                       aTTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGt                                                                                                                               >  2:1475500/1‑150 (MQ=255)
                       aTTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGt                                                                                                                               <  1:910283/150‑1 (MQ=255)
                            aaaaCGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGa                                                                                                                          <  1:1436673/150‑1 (MQ=255)
                                     tAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCa                                                                                                                 <  1:1286647/150‑1 (MQ=255)
                                     tAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCa                                                                                                                 <  2:12950/150‑1 (MQ=255)
                                         gaaagGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGtt                                                                                                             <  2:882923/150‑1 (MQ=255)
                                          aaagGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTg                                                                                                            <  2:693342/150‑1 (MQ=255)
                                            agGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGGGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGc                                                                                                          <  2:1272039/150‑1 (MQ=255)
                                              gCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCat                                                                                                        >  2:1564228/1‑150 (MQ=255)
                                               cacaTGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCata                                                                                                       >  1:166722/1‑150 (MQ=255)
                                                     aggagAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCAcc                                                                                                 >  2:1068627/1‑150 (MQ=255)
                                                      ggagAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCAccc                                                                                                >  2:72557/1‑150 (MQ=255)
                                                      ggagAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCAccc                                                                                                >  1:284675/1‑150 (MQ=255)
                                                       gagaATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCg                                                                                               <  1:1564228/150‑1 (MQ=255)
                                                        agaATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCCCGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACCGCATCATCGACGGTAATGAAGTGCGTCATGTTGGCATATACCCCCGc                                                                                              >  1:1577577/1‑150 (MQ=255)
                                                               tCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGc                                                                                       <  1:1507506/150‑1 (MQ=255)
                                                                ccAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCg                                                                                      >  1:951769/1‑150 (MQ=255)
                                                                      tatTTAACGCCAGTCGCTGGATACGATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGagag                                                                                >  2:1052306/1‑150 (MQ=255)
                                                                      tatTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGagag                                                                                >  1:1240604/1‑150 (MQ=255)
                                                                       atTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAgg                                                                               >  2:732920/1‑150 (MQ=255)
                                                                        tttAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGc                                                                              >  2:484052/1‑150 (MQ=255)
                                                                          tAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCaa                                                                            <  2:1128277/150‑1 (MQ=255)
                                                                           aaCGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAAt                                                                           >  2:435749/1‑150 (MQ=255)
                                                                              gCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGc                                                                        <  1:97945/150‑1 (MQ=255)
                                                                               ccAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCTTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCa                                                                       <  2:1310907/150‑1 (MQ=255)
                                                                                      cTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTg                                                                <  2:1567496/150‑1 (MQ=255)
                                                                                      cTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTg                                                                <  2:348259/150‑1 (MQ=255)
                                                                                        ggATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGAc                                                              <  2:967582/150‑1 (MQ=255)
                                                                                          aTACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGAccc                                                            <  1:1408812/150‑1 (MQ=255)
                                                                                             ccATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTgaga                                                                                 <  1:209535/126‑1 (MQ=255)
                                                                                             ccATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTgaga                                                                                 >  2:209535/1‑126 (MQ=255)
                                                                                             ccATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCaaa                                                         >  1:595949/1‑150 (MQ=255)
                                                                                             ccATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCaaa                                                         >  1:1221133/1‑150 (MQ=255)
                                                                                              cATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAAt                                                        >  1:384064/1‑150 (MQ=255)
                                                                                              cATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAAt                                                        >  2:508513/1‑150 (MQ=255)
                                                                                                   gTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTacc                                                   <  1:724635/150‑1 (MQ=255)
                                                                                                     aaaTGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTaccac                                                 <  2:827303/150‑1 (MQ=255)
                                                                                                      aaTGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAgg                                                                               <  1:660943/119‑1 (MQ=255)
                                                                                                      aaTGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAgg                                                                               >  2:660943/1‑119 (MQ=255)
                                                                                                      aaTGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTaccacc                                                <  1:806709/150‑1 (MQ=255)
                                                                                                      aaTGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTaccacc                                                <  1:1423174/150‑1 (MQ=255)
                                                                                                       aTGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTAccacca                                               <  1:160795/150‑1 (MQ=255)
                                                                                                           tgcAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCAcgcc                                           >  2:1206915/1‑150 (MQ=255)
                                                                                                            gcAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCAcgccg                                          >  1:784968/1‑150 (MQ=255)
                                                                                                             cAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCAcgccgc                                         >  2:327567/1‑150 (MQ=255)
                                                                                                             cAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCAcgccgc                                         >  2:1567096/1‑150 (MQ=255)
                                                                                                                  ttGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGGCATGTTGGCATATGCACCCGCCATCGCCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCg                                    <  2:983649/150‑1 (MQ=255)
                                                                                                                           cGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGTAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGc                           <  1:732920/150‑1 (MQ=255)
                                                                                                                            gCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCa                          <  1:1052306/150‑1 (MQ=255)
                                                                                                                             cccAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCAt                         >  1:1586562/1‑150 (MQ=255)
                                                                                                                                  tGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGTAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCa                    >  1:174856/1‑150 (MQ=255)
                                                                                                                                      ggTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGt                >  1:26209/1‑150 (MQ=255)
                                                                                                                                      ggTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTCCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGt                >  2:907511/1‑150 (MQ=255)
                                                                                                                                         gAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGcc             >  2:428359/1‑150 (MQ=255)
                                                                                                                                          gccgTTTCCAGACTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTTGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCg            >  2:245456/4‑150 (MQ=255)
                                                                                                                                          aatgtgtCCAGTCTGGAAACGGCATCTTGGGCTGTTTTGAAGGGCGTCTTGTTGGCATATTCACCGGCCATCGGCGTGAGAGGCAATTGCAGACCGTGAGCCAAATCTACCACCAGGCCGCCATCGCCTATCTGCATAATCAAGGTGCCg            <  2:1577577/144‑1 (MQ=255)
                                                                                                                                           ccgTTTCCAGACTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGt           >  1:1003317/3‑150 (MQ=255)
                                                                                                                                            aGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCCCGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTc          >  1:688883/1‑150 (MQ=255)
                                                                                                                                            aGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTc          >  2:13875/1‑150 (MQ=255)
                                                                                                                                            aGTTTCCAGTATGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTAACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCAAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTc          >  2:1505147/1‑150 (MQ=255)
                                                                                                                                               ttCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCAgg       >  1:572094/1‑150 (MQ=255)
                                                                                                                                                 ccAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGCCCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCg     >  2:1306537/1‑150 (MQ=255)
                                                                                                                                                 ccAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCg     >  1:843855/1‑150 (MQ=255)
                                                                                                                                                    gTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCGATg  <  1:1405237/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAATGCTGCAACTTTGTGCGCACGCCCAGTGCTGGTGAAAGTTTCCAGTCTGGAAACGGCATCTTCGTCGGTAATGAAGTGCGTCATGTTGGCATATTCACCCGCCATCGGCGTGAGAGGCAATTGCAGACCGTGACCCAAATCTACCACCACGCCGCCATCGCCAATCTGCATAATCAAGGTGCCGTCAGGCGATG  >  NC_000913/2151820‑2152117

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