Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,153,420 G→A 5.6% intergenic (‑322/+276) ppc ← / ← argE phosphoenolpyruvate carboxylase/acetylornithine deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,153,4200GA 5.6% 579.2 / 16.4 177intergenic (‑322/+276)ppc/argEphosphoenolpyruvate carboxylase/acetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base G (68/99);  new base A (2/8);  total (70/107)
Fisher's exact test for biased strand distribution p-value = 3.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.21e-01

Only 100 of 186 total aligned reads displayed.

ATTAAATTCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATC  >  NC_000913/4153273‑4153561
                                                                                                                                                   |                                                                                                                                             
aTTAAATTCACGACGCTTTAAATAAGCGTAACTTCTGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATgcg                                                                                                                                             <  2:2696212/150‑1 (MQ=255)
aTTAAATTCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATgcg                                                                                                                                             >  1:1529177/1‑150 (MQ=255)
 ttAAATTCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATgcgg                                                                                                                                            >  1:735628/1‑150 (MQ=255)
 ttAAATTCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATgcgg                                                                                                                                            <  2:1260977/150‑1 (MQ=255)
   aaaTTCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGAGCCGGATgcggcg                                                                                                                                          >  1:284517/1‑150 (MQ=255)
    aaTTCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGgcgc                                                                                                                                         <  2:1061711/150‑1 (MQ=255)
      ttCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGa                                                                                                                                       >  2:1447874/1‑150 (MQ=255)
      ttCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGa                                                                                                                                       <  1:1089924/150‑1 (MQ=255)
       tCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAg                                                                                                                                      >  1:1655072/1‑150 (MQ=255)
       tCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAg                                                                                                                                      >  1:40225/1‑150 (MQ=255)
             cgCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCtt                                                                                                                                <  2:735871/150‑1 (MQ=255)
              gCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTa                                                                                                                               >  1:2084958/1‑150 (MQ=255)
              gCTTAAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACCCCAAAGGTGTAGGTCGGATAAGATGCGCACGTATCGCAGCCGACATTATTGAGGCACTGGAGTTTGGCAACAGAGCCGGATGCGGCGCGAGCGCCTTa                                                                                                                               >  1:2268989/1‑150 (MQ=21)
               cTTTAAATAAGCGTAACTTATTGAAATTTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATAAGGCGCTCGCGccgcat                                                                                                                              <  2:1903677/150‑5 (MQ=21)
               cTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTAt                                                                                                                              >  1:134867/1‑150 (MQ=255)
               cTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATAAGGCGCTCGCGccgcat                                                                                                                              <  1:715229/150‑5 (MQ=25)
                tttAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATc                                                                                                                             >  1:236659/1‑150 (MQ=255)
                tttAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATc                                                                                                                             >  1:2729172/1‑150 (MQ=255)
                tttAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATc                                                                                                                             >  2:2330562/1‑150 (MQ=255)
                 ttAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTACCACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATcc                                                                                                                            <  1:759168/150‑1 (MQ=255)
                 ttAAATAAGCGTAACGTATGGAACTGTTAAACCATCGCCCCAAGTAACACCACAGGTGTAGGGCGGATACGATGCGCAAGTATCGCATCCGACATTCTTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATcc                                                                                                                            <  2:1305582/150‑1 (MQ=18)
                  tAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATAAGGCGCTCGCGCCGCATCCg                                                                                                                           <  1:977170/150‑1 (MQ=25)
                  tAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATAAGGCGCTCGCGCCGCATCCg                                                                                                                           <  1:835872/150‑1 (MQ=25)
                    aaTAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGc                                                                                                                         <  1:1417195/150‑1 (MQ=255)
                    aaTAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGc                                                                                                                         <  1:1417190/150‑1 (MQ=255)
                    aaTAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATAAGGCGCTCGCGCCGCATCCGGc                                                                                                                         <  2:1474809/150‑1 (MQ=21)
                    aaTAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGc                                                                                                                         >  2:2348967/1‑150 (MQ=255)
                    aaTAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGAGCCGGATGCGGCGCGAGCGCCTTATCCGGc                                                                                                                         >  2:1141988/1‑150 (MQ=255)
                     aTAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGcc                                                                                                                        >  2:929296/1‑150 (MQ=255)
                     aTAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGcc                                                                                                                        >  1:2240269/1‑150 (MQ=255)
                      tAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCt                                                                                                                       >  1:2163513/1‑150 (MQ=255)
                      tAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCt                                                                                                                       >  1:2553006/1‑150 (MQ=255)
                       aaGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTa                                                                                                                      >  1:2312658/1‑150 (MQ=255)
                          cGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGAGCCGGATGCGGAGCTAGCGCCATAACAGGACTaaac                                                                                                                   >  1:1450170/1‑147 (MQ=21)
                              aCTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCACAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTgg                                                                                                               <  1:2043757/150‑1 (MQ=255)
                              aCTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTgg                                                                                                               <  1:862305/150‑1 (MQ=255)
                                ttATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGc                                                                                                             >  1:577669/1‑150 (MQ=255)
                                 tATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCgg                                                                                                                          <  2:177970/136‑1 (MQ=255)
                                    ggAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACCGTTGGGCATCg                                                                                                         >  2:985931/1‑150 (MQ=255)
                                       aaTGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGttt                                                                                                      <  1:1733068/150‑1 (MQ=255)
                                        aTGTTAACAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTg                                                                                                     <  2:875846/150‑1 (MQ=255)
                                         tGTTAAAAAATCGCCCCAAGTAACACCAACGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGa                                                                                                    <  1:1247398/150‑1 (MQ=255)
                                          gTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAg                                                                                                   <  2:974103/150‑1 (MQ=255)
                                            tAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACCGTTTTGCAACGGATGAGTc                                                                                                 >  1:2420000/1‑150 (MQ=17)
                                             aaaaaaTCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCa                                                                                                >  1:233143/1‑150 (MQ=255)
                                              aaaaaTCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTGTGAGTCAc                                                                                               >  2:2299757/1‑150 (MQ=255)
                                                   tCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTgtcg                                                                                          <  2:1674607/150‑1 (MQ=255)
                                                     gCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTgtcggt                                                                                        <  1:1008634/150‑1 (MQ=255)
                                                         cAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGa                                                                                    >  1:2183935/1‑150 (MQ=255)
                                                           aGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATa                                                                                  <  2:1251099/150‑1 (MQ=255)
                                                                 acCAACGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGcc                                                                                                                                  <  2:1651988/96‑1 (MQ=255)
                                                                 acCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTTGGCATCTTATGAGTCCCTGTCGGTCGTATAAGATgc                                                                            >  2:2618597/1‑150 (MQ=14)
                                                                 acCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATgc                                                                            <  1:2595229/150‑1 (MQ=255)
                                                                 acCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATgc                                                                            <  1:929296/150‑1 (MQ=255)
                                                                 acCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATgc                                                                            <  1:1235858/150‑1 (MQ=255)
                                                                 acCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTAATACGGCCTACAGATGGGCATCGGTTGAGTCACTGTCGGTCGGATAAGATgc                                                                            >  1:1652017/1‑150 (MQ=255)
                                                                 acCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACTGTTGGGGGTCGGCTGAGTGGGTGTCGTTCGGGTAAGAtcc                                                                            >  2:456085/1‑148 (MQ=14)
                                                                  ccAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGATgcg                                                                           >  2:2664398/1‑150 (MQ=17)
                                                                  ccAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATgcg                                                                           <  1:681626/150‑1 (MQ=255)
                                                                   cAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATgcgc                                                                          <  2:659014/150‑1 (MQ=255)
                                                                   cAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATgcgc                                                                          >  1:1803422/1‑150 (MQ=255)
                                                                     aaGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCaa                                                                        <  1:1931590/150‑1 (MQ=40)
                                                                     aaGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCaa                                                                        <  1:1542186/150‑1 (MQ=40)
                                                                     aaGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCaa                                                                        <  1:2449987/150‑1 (MQ=40)
                                                                         tgtAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTAt                                                                    <  1:996938/150‑1 (MQ=39)
                                                                         tgtAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTAt                                                                    <  1:243566/150‑1 (MQ=39)
                                                                           tAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCg                                                                  <  1:1895295/150‑1 (MQ=38)
                                                                            aGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGc                                                                 >  2:2508580/1‑150 (MQ=38)
                                                                               tCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACATTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATc                                                              >  1:718879/1‑150 (MQ=37)
                                                                                  gATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGttt                                                                                                      >  2:557381/1‑107 (MQ=32)
                                                                                  gATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGttt                                                                                                      <  1:557381/107‑1 (MQ=32)
                                                                                  gATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGa                                                           <  2:577669/150‑1 (MQ=37)
                                                                                  gATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGa                                                           <  2:630272/150‑1 (MQ=37)
                                                                                  gATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGa                                                           <  2:2112678/150‑1 (MQ=37)
                                                                                  gATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGa                                                           <  1:2089970/150‑1 (MQ=37)
                                                                                   aTAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGAc                                                          <  1:1116797/150‑1 (MQ=37)
                                                                                            gcAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCTAGTATCGCATCCGACATTATTGCg                                                 >  2:627606/1‑150 (MQ=37)
                                                                                            gcAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCg                                                 <  1:1656626/150‑1 (MQ=37)
                                                                                             cAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCgg                                                <  2:1355916/150‑1 (MQ=37)
                                                                                                    cGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGa                                         <  1:538369/150‑1 (MQ=37)
                                                                                                       aTCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGtt                                      >  1:213532/1‑150 (MQ=37)
                                                                                                        tCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGttt                                     >  1:164568/1‑150 (MQ=37)
                                                                                                          cGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTgg                                   <  1:490866/150‑1 (MQ=37)
                                                                                                              attattGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAAc                               <  1:848245/150‑1 (MQ=37)
                                                                                                              attattGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAAc                               <  1:848207/150‑1 (MQ=37)
                                                                                                              attattGCGGCACTGGAGTTTGGCAACAGTGCAGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAAc                               <  1:848231/150‑1 (MQ=37)
                                                                                                                     cGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCg                        <  2:1548531/150‑1 (MQ=39)
                                                                                                                      tgCACTTGAGTTTGGCTACAGTTCATTATGCGGCTCGAGCGCGTTATCATGCATAATGTTGGGAATCGTTTGAGTCCCTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCgc                       <  2:331564/149‑2 (MQ=11)
                                                                                                                      ggCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCTCTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCgg                       <  2:1655072/150‑1 (MQ=255)
                                                                                                                      ggCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCgg                       <  2:483925/150‑1 (MQ=39)
                                                                                                                      ggCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCATTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCgg                       <  2:430764/150‑1 (MQ=255)
                                                                                                                      ggCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCAGGACTACAGTTGGGCATCTTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCgg                       <  2:1654750/150‑1 (MQ=255)
                                                                                                                         aCTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCTCTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATg                    <  2:1373799/150‑1 (MQ=255)
                                                                                                                            ggAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCgcaataat                 <  2:130095/150‑8 (MQ=255)
                                                                                                                             gAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATgcggc                >  1:1007051/1‑150 (MQ=39)
                                                                                                                                   ggCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCAGCGTGTGAGTCACTGTCGGTCGGATAAGATGCGCAAGGATCGCATCAGACATTATGGCGGCACGGGAGATTGGCAACAGTGCCGGATGCGGCGCGAgc          >  2:1318428/1‑150 (MQ=11)
                                                                                                                                      aaCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGcc       <  2:421401/150‑1 (MQ=39)
                                                                                                                                       aCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCt      <  1:2327377/150‑1 (MQ=39)
                                                                                                                                         aGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTa    >  2:621225/1‑150 (MQ=21)
                                                                                                                                           ggCCGGATAAGGCGCTCGCGCCGCATACCGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACAATATTTCGGCACGGGAGTTTGGCAACATTGCCGGATGCGGCGCGAGCGCCGTATc  >  2:1269785/2‑150 (MQ=14)
                                                                                                                                                   |                                                                                                                                             
ATTAAATTCACGACGCTTTAAATAAGCGTAACTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATC  >  NC_000913/4153273‑4153561

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