Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,135,575 G→A 6.3% A266V (GCC→GTC)  pitB ← metal phosphate:H(+) symporter PitB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,135,5750GA 6.3% 239.6 / 3.0 77A266V (GCC→GTC) pitBmetal phosphate:H(+) symporter PitB
Reads supporting (aligned to +/- strand):  ref base G (45/27);  new base A (3/2);  total (48/29)
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

TAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAA  >  NC_000913/3135429‑3135724
                                                                                                                                                  |                                                                                                                                                     
tAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCAt                                                                                                                                                    >  1:1348865/1‑150 (MQ=255)
  gCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATcg                                                                                                                                                  >  1:374623/1‑150 (MQ=255)
     tCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGcgg                                                                                                                                               >  1:1126447/1‑150 (MQ=255)
       aaaGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGTATCGcgggt                                                                                                                                             >  1:709970/1‑150 (MQ=255)
            tATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGcgggtacggg                                                                                                                                        >  1:122582/1‑150 (MQ=255)
                 gCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAAtt                                                                                                                                   >  1:1211285/1‑150 (MQ=255)
                     gATGACACTGAAACTCTGGTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCACTCAACTTCTGCGGCAGTGCAGGATGCTGTTGCAGGTAGTGGTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTcct                                                                                                                               <  2:1344487/150‑3 (MQ=255)
                          cAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCagtgccc                                                                                                                          >  1:608957/1‑144 (MQ=255)
                          cAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCg                                                                                                                          >  1:11972/1‑150 (MQ=255)
                           aGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGTGTACAGGTAATTTCATAGCCgg                                                                                                                         <  2:1441498/150‑1 (MQ=255)
                                  ctctGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATcgcg                                                                                                                                                >  1:390318/1‑120 (MQ=255)
                                  ctctGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATcgcg                                                                                                                                                <  2:390318/120‑1 (MQ=255)
                                    ctGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCa                                                                                                                >  1:671729/1‑150 (MQ=255)
                                    ctGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCa                                                                                                                <  1:242076/150‑1 (MQ=255)
                                    ctGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCa                                                                                                                >  1:1411586/1‑150 (MQ=255)
                                    ctGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCa                                                                                                                >  1:1405134/1‑150 (MQ=255)
                                         aCTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGACATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTg                                                                                                           >  2:1361760/1‑150 (MQ=255)
                                         aCTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGACATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTg                                                                                                           >  2:947848/1‑150 (MQ=255)
                                          cTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTga                                                                                                          <  1:378900/150‑1 (MQ=255)
                                            tGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTgacg                                                                                                        <  2:831418/150‑1 (MQ=255)
                                            tGCGTCCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCAACGCAATCAACTTATCCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTACATAGCCGGACGCATTCATATTgacg                                                                                                        >  2:694416/1‑150 (MQ=255)
                                               gTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGacgacg                                                                                                     <  2:1354706/150‑1 (MQ=255)
                                               gTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATCCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGAATCGCGGGTAAGGGTAATTTCATAGCAGGACGCATTCATATTGacgacg                                                                                                     >  2:1275682/1‑150 (MQ=255)
                                               gTGCCATCACTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGacgacg                                                                                                     >  2:1082104/1‑150 (MQ=255)
                                                   cATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGc                                                                                                 >  2:1186808/1‑150 (MQ=255)
                                                    aTCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGcc                                                                                                <  2:4659/150‑1 (MQ=255)
                                                       aGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCcagc                                                                                             >  2:190606/1‑150 (MQ=255)
                                                        gTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCcagca                                                                                            <  2:1514225/150‑1 (MQ=255)
                                                             tgctgcAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGc                                                                                       <  1:1274802/150‑1 (MQ=255)
                                                             tgctgcAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGc                                                                                       >  1:605690/1‑150 (MQ=255)
                                                             tgctgcAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGc                                                                                       <  2:796372/150‑1 (MQ=255)
                                                                  cAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATcc                                                                                  >  2:238679/1‑150 (MQ=255)
                                                                     gCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCcca                                                                               >  1:514010/1‑150 (MQ=255)
                                                                     gCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCcca                                                                               <  1:1160568/150‑1 (MQ=255)
                                                                      cAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCccac                                                                              <  1:887086/150‑1 (MQ=255)
                                                                         tGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAg                                                                           >  1:194496/1‑150 (MQ=255)
                                                                         tGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAg                                                                           >  2:211681/1‑150 (MQ=255)
                                                                         tGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAg                                                                           >  1:1205141/1‑150 (MQ=255)
                                                                          ggAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGt                                                                          >  1:192285/1‑150 (MQ=255)
                                                                            aGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTAc                                                                        <  2:11972/150‑1 (MQ=255)
                                                                             ggTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTAcc                                                                       <  2:122582/150‑1 (MQ=255)
                                                                                 ccATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGGAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAg                                                                                                  >  2:183954/1‑119 (MQ=255)
                                                                                 ccATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGGAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAg                                                                                                  <  1:183954/119‑1 (MQ=255)
                                                                                 ccATCGCAATCAACTTCTGCGGAAGTTCAGTATACTGATGCAGGTAGTGTTAGAAGTTAATAACTACATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCa                                                                   >  2:1441469/1‑150 (MQ=255)
                                                                                     cGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTAc                                                               >  2:695150/1‑150 (MQ=255)
                                                                                       cAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCa                                                             >  1:449947/1‑150 (MQ=255)
                                                                                                 cTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATccc                                                   >  2:241742/1‑150 (MQ=255)
                                                                                                       cAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTg                                             <  1:1082104/150‑1 (MQ=255)
                                                                                                       cAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTg                                             <  1:888424/150‑1 (MQ=255)
                                                                                                         gTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTc                                           >  1:525724/1‑150 (MQ=255)
                                                                                                              ggATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCcgtcg                                      >  2:780264/1‑150 (MQ=255)
                                                                                                                aTGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGtt                                    <  2:211599/150‑1 (MQ=255)
                                                                                                                 tGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGCTCCCTTTTTGTCCGTCGTTc                                   >  1:619524/1‑150 (MQ=255)
                                                                                                                 tGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTc                                   <  2:1405134/150‑1 (MQ=255)
                                                                                                                     gTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCgcgc                               >  2:340985/1‑150 (MQ=255)
                                                                                                                       tGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCg                             >  1:825653/1‑150 (MQ=255)
                                                                                                                        gCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGt                            >  2:328234/1‑150 (MQ=255)
                                                                                                                        gCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGt                            >  1:185807/1‑150 (MQ=255)
                                                                                                                        gCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTAATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACAAGCATTACCAGGCCGATCCCTTTTTCTCCGTCGTTcgctacgt                            >  2:1337588/1‑145 (MQ=255)
                                                                                                                         cAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTg                           >  1:456693/1‑150 (MQ=255)
                                                                                                                          aGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGc                          >  1:110261/1‑150 (MQ=255)
                                                                                                                          aGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGc                          <  1:1232084/150‑1 (MQ=255)
                                                                                                                               gtgtTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGaaaa                     >  2:1206807/1‑150 (MQ=255)
                                                                                                                                 gtTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACTCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCAGAGCGaaaaag                   >  1:786955/1‑148 (MQ=255)
                                                                                                                                 gtTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAacg                   >  1:325211/1‑150 (MQ=255)
                                                                                                                                  ttCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAacgc                  <  2:1411586/150‑1 (MQ=255)
                                                                                                                                  ttCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAacgc                  <  1:1450715/150‑1 (MQ=255)
                                                                                                                                  ttCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAacgc                  <  2:671729/150‑1 (MQ=255)
                                                                                                                                    cGAAGTTGGTAACGGCATCGCGGGTATGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAacgcca                >  2:1234190/1‑150 (MQ=255)
                                                                                                                                    cGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAacgcca                >  2:219508/1‑150 (MQ=255)
                                                                                                                                          tGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCg          <  1:780264/150‑1 (MQ=255)
                                                                                                                                          tGGTAACGACATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCg          <  1:947848/150‑1 (MQ=255)
                                                                                                                                          tGGTAACGACATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCg          <  1:1361760/150‑1 (MQ=255)
                                                                                                                                           ggTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGCGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGc         >  2:809747/1‑150 (MQ=255)
                                                                                                                                             tAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGAATTCATATTGACGACGAAGCCAGCAGGGTCAATCCACACCAGTACCAGCATTACCACCCCGAACCCTTTATGTCCGTCGTTCGCGCCCTGCGAAAACGCCACGCcagcag       >  2:383288/1‑150 (MQ=255)
                                                                                                                                             tAACGGCAGCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCTGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGGGCGAAAACGCCACGCCCGCAg       <  2:980253/150‑1 (MQ=255)
                                                                                                                                                cGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACGAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCgg    <  2:374623/150‑1 (MQ=255)
                                                                                                                                                 ggCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGa   >  2:112068/1‑150 (MQ=255)
                                                                                                                                                  gCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGaa  <  2:605690/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAA  >  NC_000913/3135429‑3135724

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