Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,537,076 C→G 6.7% L54L (CTG→CTC)  pdxK ← pyridoxal kinase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,537,0760CG 6.7% 483.1 / 17.3 144L54L (CTG→CTC) pdxKpyridoxal kinase I
Reads supporting (aligned to +/- strand):  ref base C (72/61);  new base G (6/4);  total (78/66)
Fisher's exact test for biased strand distribution p-value = 7.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 144 total aligned reads displayed.

GATTTGCGATGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATCGTTAAACAACA  >  NC_000913/2536929‑2537224
                                                                                                                                                   |                                                                                                                                                    
gATTTGCGATGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAg                                                                                                                                                    >  1:2290787/1‑150 (MQ=255)
gATTTGCGATGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAg                                                                                                                                                    >  1:124240/1‑150 (MQ=255)
  tttGCGATGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAACTGTCATAATGCGGCGTATTGCTCAGCa                                                                                                                                                  <  1:599940/150‑1 (MQ=255)
        aTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTGAGCAATAcgc                                                                                                                                            <  2:1290512/150‑3 (MQ=255)
           ccGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTcg                                                                                                                                         >  2:534781/1‑150 (MQ=255)
           ccGTTCCCATATAGCCAGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTcg                                                                                                                                         >  1:2300081/1‑150 (MQ=255)
             gTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTcggc                                                                                                                                       <  2:1590796/150‑1 (MQ=255)
                 ccATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTcggcacgg                                                                                                                                   <  2:994171/150‑1 (MQ=255)
                      taGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAg                                                                                                                              <  2:1603095/150‑1 (MQ=255)
                          ccGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACAt                                                                                                                          >  2:61402/1‑150 (MQ=255)
                           cGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACAtt                                                                                                                         <  2:1243573/150‑1 (MQ=255)
                            gTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTc                                                                                                                        <  2:250391/150‑1 (MQ=255)
                              ggTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAg                                                                                                                      <  1:1246013/150‑1 (MQ=255)
                                ttACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGc                                                                                                                    >  1:1377594/1‑150 (MQ=255)
                                    aGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGtt                                                                                                                >  2:2316085/1‑150 (MQ=255)
                                     gCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTc                                                                                                               >  1:159368/1‑150 (MQ=255)
                                     gCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTc                                                                                                               >  2:2471392/1‑150 (MQ=255)
                                     gCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTc                                                                                                               >  1:1780538/1‑150 (MQ=255)
                                     gCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTc                                                                                                               >  2:1546597/1‑150 (MQ=255)
                                     gCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTc                                                                                                               <  2:2376162/150‑1 (MQ=255)
                                      cACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAACCCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCt                                                                                                              >  1:761136/1‑150 (MQ=255)
                                       aCGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTg                                                                                                             >  2:1455874/1‑150 (MQ=255)
                                        cGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGt                                                                                                            >  2:1068546/1‑150 (MQ=255)
                                         gAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGtt                                                                                                           >  2:2475231/1‑150 (MQ=255)
                                           aGTTGGCGCAGCGCATCACGCTCCTGAATCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTg                                                                                                         <  2:1317874/150‑1 (MQ=255)
                                           aGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTg                                                                                                         <  1:2305344/150‑1 (MQ=255)
                                           aGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTg                                                                                                         <  1:573558/150‑1 (MQ=255)
                                                gcgcagcgcaTCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGc                                                                                                    >  1:1698529/1‑150 (MQ=255)
                                                    agcgcaTCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGc                                                                                                <  1:1168026/150‑1 (MQ=255)
                                                         aTCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGc                                                                                           <  1:1385152/150‑1 (MQ=255)
                                                          tCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCa                                                                                          >  1:287390/1‑150 (MQ=255)
                                                            aCGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAAt                                                                                        <  2:1266801/150‑1 (MQ=255)
                                                             cGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATg                                                                                       >  1:1239171/1‑150 (MQ=255)
                                                              gCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGc                                                                                      >  2:755344/1‑150 (MQ=255)
                                                                tCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCaacagca                                                                                    <  1:1315448/150‑6 (MQ=255)
                                                                  cTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCAGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGtt                                                                                  >  2:171463/1‑150 (MQ=255)
                                                                   tGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGc                                                                                      <  1:805049/145‑1 (MQ=255)
                                                                    gatgcgcACTCAAATATCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGc                                                                                <  2:1409899/147‑1 (MQ=255)
                                                                      aGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGccc                                                                              >  2:2198951/1‑150 (MQ=255)
                                                                         gcacgcaAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACg                                                                           >  2:2124992/1‑150 (MQ=255)
                                                                          cacgcaAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCAc                                                                            <  1:2389366/148‑1 (MQ=255)
                                                                          cacgcaAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCAc                                                                            >  2:2389366/1‑148 (MQ=255)
                                                                            cgcaAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTg                                                                        >  1:2108806/1‑150 (MQ=255)
                                                                             gcaAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGAGATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGc                                                                       <  2:284113/150‑1 (MQ=255)
                                                                             gcaAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGc                                                                       <  2:2249721/150‑1 (MQ=255)
                                                                              caAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGcc                                                                      >  2:542676/1‑150 (MQ=255)
                                                                                aaTAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGt                                                                    <  2:1415085/150‑1 (MQ=255)
                                                                                  tAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGCCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTaa                                                                  >  2:1797278/1‑150 (MQ=255)
                                                                                   aGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGTCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTaaa                                                                 <  2:1590289/150‑1 (MQ=255)
                                                                                   aGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGGCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTaaa                                                                 <  1:727485/150‑1 (MQ=255)
                                                                                    gCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCAACGCTGCCGTAAac                                                                >  2:1792675/1‑150 (MQ=255)
                                                                                    gCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTCATAGCAGGCACGGCAATGCTGTTGCCCACCCTGCCGTAAac                                                                >  2:2380607/1‑150 (MQ=255)
                                                                                         tAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCt                                                           <  1:1970857/150‑1 (MQ=255)
                                                                                             ccATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGa                                                       >  2:103877/1‑150 (MQ=255)
                                                                                               aTTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACt                                                     <  2:2162067/150‑1 (MQ=255)
                                                                                                ttCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTg                                                    <  1:498291/150‑1 (MQ=255)
                                                                                                  cGTCCTGAATCGAACCAGCGTAGAAAGTGTCATAATGCGGCGTATTTCTCAGCAATACCGTCGGCAGGGAAAAGACATTCAGGCCGTTCTGTTTTATAGCATGCACGGCAATGATGTTGCCCACGCTGTCTTAAACGACCTGCGACTGCa                                                  <  2:2277154/150‑1 (MQ=255)
                                                                                                  cGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCa                                                  >  2:688927/1‑150 (MQ=255)
                                                                                                  cGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCa                                                  >  2:803823/1‑150 (MQ=255)
                                                                                                    tCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACg                                                <  1:534781/150‑1 (MQ=255)
                                                                                                    cccGGCAGCGCACCACCCTAGACCGCGGCATAATGCGGCGTATTGCGCAGCAATACCGGCGGCACGGCAAAGACATTCAGGCCGGCCTGCTTGCTAGCCGGCACGGCATTGCTGGGGCCCACGCGGCCGTAAACCCCCTGCGACTGCACg                                                <  2:999340/149‑1 (MQ=255)
                                                                                                     ccGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACt                                                     <  2:413951/144‑1 (MQ=255)
                                                                                                     ccGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACt                                                     >  1:413951/1‑144 (MQ=255)
                                                                                                      cGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGc                                              >  2:446292/1‑150 (MQ=255)
                                                                                                      cGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGc                                              >  2:81239/1‑150 (MQ=255)
                                                                                                      cGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGc                                              >  1:1347585/1‑150 (MQ=255)
                                                                                                       ggAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGcg                                             >  2:2097764/1‑150 (MQ=255)
                                                                                                       ggAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGGCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGcg                                             <  2:2124850/150‑1 (MQ=255)
                                                                                                        gAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGcga                                            >  1:1254454/1‑150 (MQ=255)
                                                                                                         aaTCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGcgac                                           <  2:2108806/150‑1 (MQ=255)
                                                                                                         aaTCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCACACGCTGCCGTAAACCACCTGCGACTGCACGGcgac                                           >  1:1130244/1‑150 (MQ=255)
                                                                                                          aTCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCGCGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGcgacg                                          <  1:643882/150‑1 (MQ=255)
                                                                                                            cGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGat                                        >  1:1988621/1‑150 (MQ=255)
                                                                                                             gCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGata                                       >  1:642726/1‑150 (MQ=255)
                                                                                                              caccacCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGatat                                      >  1:1595591/1‑150 (MQ=255)
                                                                                                                 caccGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCg                                   <  2:159368/150‑1 (MQ=255)
                                                                                                                    cGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCt                                <  1:1468491/150‑1 (MQ=255)
                                                                                                                    cGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCt                                >  1:777870/1‑150 (MQ=255)
                                                                                                                      tAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGc                              <  1:101515/150‑1 (MQ=255)
                                                                                                                            gtgtCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGcc                        >  2:1750594/1‑150 (MQ=255)
                                                                                                                                cATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTAc                    >  2:2274873/1‑150 (MQ=255)
                                                                                                                                cATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTAc                    >  1:345588/1‑150 (MQ=255)
                                                                                                                                 aTAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTAct                   >  1:819502/1‑150 (MQ=255)
                                                                                                                                  tAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTActc                  >  2:1022159/1‑150 (MQ=255)
                                                                                                                                      gcggcgTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTAt              >  1:716116/1‑150 (MQ=255)
                                                                                                                                       gacggTATTGCTGAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATc             >  1:571976/5‑150 (MQ=255)
                                                                                                                                       gacggTATTGCTGAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATc             >  2:1223754/5‑150 (MQ=255)
                                                                                                                                        acggTATTGCTGAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATCg            >  1:964429/4‑150 (MQ=255)
                                                                                                                                         cggTATTGCTGAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATCGt           >  2:837485/3‑150 (MQ=255)
                                                                                                                                          cgTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATCGtt          <  1:853842/150‑1 (MQ=255)
                                                                                                                                           gTATTGCTGAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATCGTTa         >  2:1849607/1‑150 (MQ=255)
                                                                                                                                            tATTGCTCAGCAATAGGGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTTTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACTGCGACGATATCCGCCTGCAGTGCCCTACTCTTATCGTTaa        <  1:2236010/150‑1 (MQ=255)
                                                                                                                                             aTTGCTCAGCAATACCGCCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATCGTTaaa       <  1:1797278/150‑1 (MQ=255)
                                                                                                                                              ttGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATCGTTAaac      >  2:1185236/1‑150 (MQ=255)
                                                                                                                                              ttGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATCGTTAaac      <  2:642726/150‑1 (MQ=255)
                                                                                                                                              ttGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATCGTTAaac      <  1:1455874/150‑1 (MQ=255)
                                                                                                                                              ttGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATCGTTAaac      <  2:2353279/150‑1 (MQ=255)
                                                                                                                                                gCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATCGTTAaacaa    >  1:604454/1‑150 (MQ=255)
                                                                                                                                                 cTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATCGTTAaacaac   >  1:1352774/1‑150 (MQ=255)
                                                                                                                                                  tCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATCGTTAaacaaca  <  1:2359876/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GATTTGCGATGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATTCGTCCGGAATCGCACCACCGTAGAAAGTGTCATAATGCGGCGTATTGCTCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATCGTTAAACAACA  >  NC_000913/2536929‑2537224

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: ‑