Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,896,474 C→A 5.3% P102Q (CCG→CAG)  nudL → putative CoA pyrophosphohydrolase, weak 3‑phosphohydroxypyruvate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,896,4740CA 5.3% 296.6 / 2.9 93P102Q (CCG→CAG) nudLputative CoA pyrophosphohydrolase, weak 3‑phosphohydroxypyruvate phosphatase
Reads supporting (aligned to +/- strand):  ref base C (53/34);  new base A (4/2);  total (57/36)
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

GATTCATCTGCGTAAACACGCTGGACAAGTGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACATCTGGGTCGTTATCACC  >  NC_000913/1896325‑1896620
                                                                                                                                                     |                                                                                                                                                  
gATTCATCTGCGTAAACACGCTGGACAAGTGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTgcc                                                                                                                                                    >  2:989602/1‑150 (MQ=255)
gATTCATCTGCGTAAACACGCTGGACAAGTGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTgcc                                                                                                                                                    <  2:301491/150‑1 (MQ=255)
  ttCATCTGCGTAAACACGCTGGACAAGTGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTgccgc                                                                                                                                                  >  1:1509023/1‑150 (MQ=255)
     ctCTGCGTCACCCCGCTGGACCAGTGGCATTCCCTGGAGGGGCAGTCGAGGACACGGCCGCAGCAGCTCGCGCCGCCGCGCGGCGCGAAGCGGAAGAAGCGGGCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCg                                                                                                                                               <  1:704728/149‑1 (MQ=255)
          cGTAAACACGCTGGACAAGTGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGAt                                                                                                                                          >  1:254443/1‑150 (MQ=255)
          cGTAAACACGCTGGACAAGTGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGAt                                                                                                                                          <  2:359701/150‑1 (MQ=255)
             aaaCACGCTGGACAAGTGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGc                                                                                                                                       <  1:214562/150‑1 (MQ=255)
             aaaCACGCTGGACAAGTGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGc                                                                                                                                       >  2:5182/1‑150 (MQ=255)
              aaCACGCTGGACAAGTGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGACGGGCGGCAGCACGCCGATaact                                                                                                                                      <  2:513352/150‑4 (MQ=255)
              aaCACGCTGGACAAGTGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGACGGGCGGCAGCACGCCGATaact                                                                                                                                      <  2:1004815/150‑4 (MQ=255)
              aaCACGCTGGAAAATTGTCATTCCCTGCACGTGCAGTCTATCACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAATCTGAAGAAGAGGTCGCTAAACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCAGCCCGTAGaaaccg                                                                                                                                      >  2:763030/1‑145 (MQ=255)
                  cGCTGGACAAGTGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCAc                                                                                                                                  >  2:139109/1‑150 (MQ=255)
                     tGGACAAGTGGCATTCCCTGGAGGTGCAGTCGATGCCACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTgg                                                                                                                               >  2:731875/1‑150 (MQ=255)
                     tGGACAAGTGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTgg                                                                                                                               >  2:972106/1‑150 (MQ=255)
                      ggACAAGTGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGc                                                                                                                              >  2:739200/1‑150 (MQ=255)
                         cAAGTGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTAc                                                                                                                           >  1:1053195/1‑150 (MQ=255)
                           aGTGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCAGCCCGTCGATAGCGTCACTGGCTACCa                                                                                                                         >  2:284397/1‑150 (MQ=255)
                             tGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAgg                                                                                                                       <  2:1481964/150‑1 (MQ=255)
                              ggaaTTCCCTGGAGGTTCATTCGATGACACGTACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGt                                                                                                                      <  2:895964/147‑1 (MQ=255)
                               gCATTCCCTGGAGGTGCAGTCGATGACACCGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGATGTTTTCGGCGAGCTGCCGCCCGTCGATAGCGTCTCTGGCTACCAGGTa                                                                                                                     >  2:57985/1‑150 (MQ=255)
                                 aTTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAAc                                                                                                                   >  2:999132/1‑150 (MQ=255)
                                  ttCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAAcc                                                                                                                  <  2:1053195/150‑1 (MQ=255)
                                      cTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAg                                                                                                              >  2:98126/1‑150 (MQ=255)
                                       tGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGt                                                                                                             <  2:1274876/150‑1 (MQ=255)
                                        ggAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCATCGATAGCGTCACTGGCTACCAGGTAACCCCAGTg                                                                                                            >  1:1581944/1‑150 (MQ=255)
                                            gTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCg                                                                                                        <  2:703165/150‑1 (MQ=255)
                                            gTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCg                                                                                                        >  1:1431234/1‑150 (MQ=255)
                                             ggCAGTCGATGACACGGACGCATCAGATATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCgg                                                                                                       <  1:751393/149‑1 (MQ=255)
                                              gCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGc                                                                                                      >  2:1030166/1‑150 (MQ=255)
                                                aGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCAt                                                                                                    <  1:556553/150‑1 (MQ=255)
                                                   cGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCTGTGGTCGGCATTAt                                                                                                 >  1:1626912/1‑150 (MQ=255)
                                                     aggACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATcc                                                                                               <  2:1442335/148‑1 (MQ=255)
                                                      tGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATccc                                                                                              >  2:148233/1‑150 (MQ=255)
                                                       gACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGACGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATcccg                                                                                             <  2:1415342/150‑1 (MQ=255)
                                                        acacGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTACCCCCAGTGGTCGGCATTATcccgc                                                                                            >  1:55009/1‑150 (MQ=255)
                                                           cGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATcccgcccg                                                                                         >  1:149306/1‑150 (MQ=255)
                                                                gCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTg                                                                                    >  1:480143/1‑150 (MQ=255)
                                                                 cATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGATGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTTTCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTAGCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGCTCTGc                                                                                   >  1:1169207/1‑150 (MQ=255)
                                                                     aGCTATCGCCGCCGCGCTCCTCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGATGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTa                                                                               >  2:852042/1‑150 (MQ=255)
                                                                      gCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTAt                                                                              <  2:480143/150‑1 (MQ=255)
                                                                        tagcgccgccgcGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTATGCAGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATcg                                                                            <  2:275809/147‑1 (MQ=255)
                                                                          tCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATcgcg                                                                          >  1:1667691/1‑150 (MQ=255)
                                                                            gccgccgcGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGcc                                                                        >  1:609430/1‑150 (MQ=255)
                                                                            gccgccgcGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGcc                                                                        >  2:1339650/1‑150 (MQ=255)
                                                                                ccgcGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTg                                                                    <  2:1254811/150‑1 (MQ=255)
                                                                                      tGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTAACCCGCCCGATCTGCCGAATCGCGCCAGTGAAGATg                                                              >  2:757637/1‑150 (MQ=255)
                                                                                         gcgAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAg                                                           <  1:450198/150‑1 (MQ=255)
                                                                                         gcgAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAg                                                           <  1:5182/150‑1 (MQ=255)
                                                                                          cgAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGt                                                          <  1:148233/150‑1 (MQ=255)
                                                                                              gCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTcg                                                      >  1:1341893/1‑150 (MQ=255)
                                                                                                tGAAGAAGAGGTCGCTATACCGCCTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTcggcg                                                   >  2:531839/1‑150 (MQ=255)
                                                                                                tGAAGAAGAGGTCGCTATACAGCCTTCCGCCGTGGAAGTTATCTGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCATTGGTCGGCATGATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAATtcgaggc                                                    >  2:229912/1‑145 (MQ=255)
                                                                                                 gaagaaGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGGGAAGATGAAGTCTcggcg                                                   >  1:202253/1‑150 (MQ=255)
                                                                                                    gaagaGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGgtg                                                >  2:1450570/1‑150 (MQ=255)
                                                                                                     aagaGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGgtgt                                               >  1:1338784/1‑150 (MQ=255)
                                                                                                     aagaGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGgtgt                                               >  1:809629/1‑150 (MQ=255)
                                                                                                      agagGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGtt                                              <  2:1621568/150‑1 (MQ=255)
                                                                                                       gagGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGttt                                             <  1:98126/150‑1 (MQ=255)
                                                                                                        agGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTg                                            <  2:1509023/150‑1 (MQ=255)
                                                                                                         ggTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTACCCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGGATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGa                                           >  1:225712/1‑150 (MQ=255)
                                                                                                         ggTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGa                                           >  2:1357827/1‑150 (MQ=255)
                                                                                                         ggTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGa                                           <  2:1667691/150‑1 (MQ=255)
                                                                                                           tCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGaaa                                         <  2:1431234/150‑1 (MQ=255)
                                                                                                            cGCTATACCGCCTTCCGACGTTGAAGTTATCGGCGTGCTGCCGCCCTTCGATAGCGTCACTGGCTACCAGGTTACCCCAGTTGTCGGCATTATCCCGCCCGATCTGCCTTATCGCGCAAGTGAAGATGAAGTCTCGGCGGTGTTTGAAAt                                        <  1:999255/150‑1 (MQ=255)
                                                                                                              cTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGc                                      >  1:1186128/1‑150 (MQ=255)
                                                                                                              cTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGAAAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCCCCCGATCTGCCGAATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGc                                      >  1:1440843/1‑150 (MQ=255)
                                                                                                                 taCCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGc                                   >  1:1146587/1‑150 (MQ=255)
                                                                                                                  aCCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCt                                  >  2:843722/1‑150 (MQ=255)
                                                                                                                   ccGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCACGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGCTGAAGTCTCGGGGGTGTTTGAAGTGCCGCTc                                 >  1:101061/1‑150 (MQ=255)
                                                                                                                     gCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGc                               >  2:450820/1‑150 (MQ=255)
                                                                                                                     gCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGACGGTGTTTGAATGCCGCTCGcc                              <  2:1253191/150‑1 (MQ=255)
                                                                                                                      ccTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGcc                              >  1:1685210/1‑150 (MQ=255)
                                                                                                                         tCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCTGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAg                           <  2:1626912/150‑1 (MQ=255)
                                                                                                                           cgccgTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGc                         >  2:1691769/1‑150 (MQ=255)
                                                                                                                            gccgTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCa                        <  1:1030166/150‑1 (MQ=255)
                                                                                                                              cgTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCAtt                      >  2:676298/1‑150 (MQ=255)
                                                                                                                                 tGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTTGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACa                   <  1:55973/150‑1 (MQ=255)
                                                                                                                                 tGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACa                   >  2:1256868/1‑150 (MQ=255)
                                                                                                                                  gAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACAt                  >  2:360847/1‑150 (MQ=255)
                                                                                                                                  gAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACAt                  >  2:1606763/1‑150 (MQ=255)
                                                                                                                                  gAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACAt                  >  2:1540853/1‑150 (MQ=255)
                                                                                                                                  gAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACAt                  >  2:1007257/1‑150 (MQ=255)
                                                                                                                                   aaGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACATc                 <  2:1387140/150‑1 (MQ=255)
                                                                                                                                    cgctATCGACGGGCGGCAGCACGCCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACATCt                >  2:1393621/4‑150 (MQ=255)
                                                                                                                                     ggtATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACATCTg               <  2:706691/148‑1 (MQ=255)
                                                                                                                                      ttATCGGCGTGCTGCCGCCAGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGTCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACATCTgg              >  2:226863/1‑150 (MQ=255)
                                                                                                                                      ctATCGACGGGCGGCAGCACGCCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACATCTgg              >  2:1173844/2‑150 (MQ=255)
                                                                                                                                          cGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACATCTGGGTCg          >  1:391773/1‑150 (MQ=255)
                                                                                                                                          cGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACATCTGGGTCg          <  2:606685/150‑1 (MQ=255)
                                                                                                                                            gCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACATCTGGGTCGtt        >  1:159791/1‑150 (MQ=255)
                                                                                                                                            gCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACATCTGGGTCGtt        <  2:1169207/150‑1 (MQ=255)
                                                                                                                                             cGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAATTGCCGCTCGCCCAGGCATTACATCTGGGTCGtga       >  2:1055588/1‑148 (MQ=255)
                                                                                                                                             cGTGCTGCCGCCCGTCGAAAGCGTCACTGGCTTCCAGGTAACCCCAGTGGGCGGCATTATCCCGCACGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATTCCGCTCGCCCTGGCATTACATCTGGGTCGTTa       <  1:284397/150‑1 (MQ=255)
                                                                                                                                               ggctgcCGCCCGTCGATAGCGTCACTGGCTTCCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTAGCGCGCCAGTGAAGATGAAGTCTCGGCGTTGTTTGAAATGCCGCTCGCCCAGGCATTACATCTGTGTCGTTATc     <  2:1440843/149‑1 (MQ=255)
                                                                                                                                                 ctgcCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACATCTGGGTCGTTATCAc   >  1:1117191/1‑150 (MQ=255)
                                                                                                                                                  tgcCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGATGATGCAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACATCTGGGTCGTTATCGcc  >  1:256740/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GATTCATCTGCGTAAACACGCTGGACAAGTGGCATTCCCTGGAGGTGCAGTCGATGACACGGACGCATCAGCTATCGCCGCCGCGCTGCGCGAAGCTGAAGAAGAGGTCGCTATACCGCCTTCCGCCGTTGAAGTTATCGGCGTGCTGCCGCCCGTCGATAGCGTCACTGGCTACCAGGTAACCCCAGTGGTCGGCATTATCCCGCCCGATCTGCCGTATCGCGCCAGTGAAGATGAAGTCTCGGCGGTGTTTGAAATGCCGCTCGCCCAGGCATTACATCTGGGTCGTTATCACC  >  NC_000913/1896325‑1896620

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