Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,908,248 G→T A40E (GCG→GAG)  pstB ← phosphate ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,908,2480GT98.8% 312.7 / NA 82A40E (GCG→GAG) pstBphosphate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (25/56);  total (25/57)
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

CAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTACTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGAT  >  NC_000913/3908100‑3908393
                                                                                                                                                    |                                                                                                                                                 
cagcagTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTcc                                                                                                                                                 <  1:712839/151‑1 (MQ=255)
   cagTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGtt                                                                                                                                              <  2:454338/151‑1 (MQ=255)
    agTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTa                                                                                                                                             <  2:798650/151‑1 (MQ=255)
      tGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTAcc                                                                                                                                           >  2:1150365/1‑151 (MQ=255)
         gATATCCTGAGAGTTGGTCAGGATGTTGTCTCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTgg                                                                                                                                        <  2:754168/151‑1 (MQ=255)
            atCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGttt                                                                                                                                     <  2:1119520/151‑1 (MQ=255)
            atCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGttt                                                                                                                                     <  1:928185/151‑1 (MQ=255)
            atCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGtt                                                                                                                                      <  1:470177/150‑1 (MQ=255)
              ccTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGttttt                                                                                                                                   >  2:1068264/1‑151 (MQ=255)
                 gagaGTTGGTCAGGATGTTGTCGCCCTCAAGCAGAATTTCACCTCCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTa                                                                                                                                  <  1:93865/149‑1 (MQ=255)
                 gagaGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAg                                                                                                                                 <  1:787108/150‑1 (MQ=255)
                  agagTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCg                                                                                                                               <  2:457293/151‑1 (MQ=255)
                   gagTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGa                                                                                                                              <  1:1026756/151‑1 (MQ=255)
                     gTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGat                                                                                                                             >  2:736128/1‑150 (MQ=255)
                          tCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATcc                                                                                                                         >  2:985081/1‑149 (MQ=255)
                           cAGGATGTTGTCGCCATCAAGCAGAATTTCACCTCCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGCTTTCCCGCAGCCGGACGGCCCGATAACCTCCGTTACCTGGTTTTTAGCGATATCCAgg                                                                                                                      <  1:845062/151‑1 (MQ=255)
                                 ggtgtCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATGTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAgct                                                                                                                     <  1:67875/144‑3 (MQ=255)
                                    gtCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGtt                                                                                                              <  1:257774/150‑1 (MQ=255)
                                    gtCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGtt                                                                                                              <  1:1267158/150‑1 (MQ=255)
                                      cGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGttttt                                                                                                           <  2:878712/151‑1 (MQ=255)
                                       gCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTc                                                                                                          >  2:981231/1‑151 (MQ=255)
                                        ccATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTc                                                                                                          >  2:603533/1‑150 (MQ=255)
                                         cATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAg                                                                                                        <  2:1051954/151‑1 (MQ=255)
                                              aGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCAt                                                                                                   <  2:803985/151‑1 (MQ=255)
                                               gCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATg                                                                                                  <  1:752895/151‑1 (MQ=255)
                                               gCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATg                                                                                                  <  1:1304964/151‑1 (MQ=255)
                                                   aaTTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGa                                                                                                <  2:833472/149‑1 (MQ=255)
                                                   aaTTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAAt                                                                                              <  2:1255751/151‑1 (MQ=255)
                                                        cACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGcc                                                                                         >  1:886577/1‑151 (MQ=255)
                                                          ccTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGcc                                                                                         >  1:1205756/1‑149 (MQ=255)
                                                             tCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCgta                                                                                      <  2:980035/149‑1 (MQ=255)
                                                              ccGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCgtagt                                                                                    <  2:64720/150‑1 (MQ=255)
                                                               cGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGATTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCgtagt                                                                                    >  2:161058/1‑149 (MQ=255)
                                                                gCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGtagtag                                                                                  <  2:681062/150‑1 (MQ=255)
                                                                gCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAAGGCATGGAATTTGCCgtagt                                                                                    >  2:184950/1‑148 (MQ=255)
                                                                  aCGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCGCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTGTTCAGGGCATGGAATTTGCCggagta                                                                                   >  1:146003/1‑147 (MQ=255)
                                                                  aCGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAg                                                                               <  1:543176/151‑1 (MQ=255)
                                                                    gctgctCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGtt                                                                             <  2:301982/151‑1 (MQ=255)
                                                                          ccGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCaaa                                                                         <  2:236926/149‑1 (MQ=255)
                                                                          ccGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCaaa                                                                         <  1:828209/149‑1 (MQ=255)
                                                                          ccGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAAt                                                                        <  2:857190/150‑1 (MQ=255)
                                                                            gggTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTa                                                                      >  2:888602/1‑150 (MQ=255)
                                                                            gggTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTAc                                                                     >  2:456955/1‑151 (MQ=255)
                                                                                aCAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGCCGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGaa                                                                  <  2:286185/150‑1 (MQ=255)
                                                                                  aGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCt                                                               >  2:827572/1‑151 (MQ=255)
                                                                                   gTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTg                                                              >  2:887543/1‑151 (MQ=255)
                                                                                     tCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGa                                                             <  1:795544/150‑1 (MQ=255)
                                                                                       aaaCATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAAtt                                                          >  1:557229/1‑151 (MQ=255)
                                                                                         aCATTTTGTTGAAGGTACGCAGCAGCGTCGCTTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAAtttt                                                        <  2:711232/151‑1 (MQ=255)
                                                                                         aCATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAAtttt                                                        <  2:264923/151‑1 (MQ=255)
                                                                                            ttttGTTGACGGTACGCAGCAGCGTCGATTTACAGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGAAGAAGTTAAAATTACGAACCTGAATTTTACt                                                     >  1:634741/1‑151 (MQ=255)
                                                                                             tttGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTACt                                                     <  2:537556/150‑1 (MQ=255)
                                                                                               tgttgAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTACTCg                                                   <  1:897222/150‑1 (MQ=255)
                                                                                                 ttgAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAACTTGCCGTAGAAGAAGTTCAAATTACGAACCTGAATTTTACTCgggg                                                >  1:735196/1‑151 (MQ=255)
                                                                                                   gAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTACTCGGGGc                                               <  1:8939/150‑1 (MQ=255)
                                                                                                   gAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTACTCGGGGCa                                              <  1:246880/151‑1 (MQ=255)
                                                                                                    aaGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTACTCGGGGc                                               >  1:323573/1‑149 (MQ=255)
                                                                                                     aGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTACTCGGGGCAGt                                            <  2:1025151/151‑1 (MQ=255)
                                                                                                        tACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTACTCGGGGCAGttt                                          <  1:1150397/150‑1 (MQ=255)
                                                                                                         aCGGAGCAGCGTCGATTTACCGCAGCCGGACGGTCCGATAAACGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTACTCGGGGCAGTTTCa                                        <  1:337647/151‑1 (MQ=255)
                                                                                                         aCGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTACTCGGGGCAGTTTc                                         >  2:906725/1‑150 (MQ=255)
                                                                                                         aCGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTACTCGGGGCAGTTTc                                         <  1:509079/150‑1 (MQ=255)
                                                                                                            cagcagCGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTACTCGGGGCAGTTTCAAcc                                     <  2:1062727/151‑1 (MQ=255)
                                                                                                              gcagcGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTACTCGGGGCAGTTTCAACCAt                                   <  1:505389/151‑1 (MQ=255)
                                                                                                                 gcGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTACTCGGGGCAGTTTCAACCATACt                                <  1:1219786/151‑1 (MQ=255)
                                                                                                                  cGTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTACTCGGGGCAGTTTCAACCATACTc                               >  2:1143143/1‑151 (MQ=255)
                                                                                                                   gTCGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTACTCGGGGCAGTTTCAACCATACTCa                              <  1:922694/151‑1 (MQ=255)
                                                                                                                     cGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTACTCGGGGCAGTTTCAACCATACTCAtt                            <  1:549979/151‑1 (MQ=255)
                                                                                                                     cGATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTACTCGGGGCAGTTTCAACCATACTCAtt                            <  2:989042/151‑1 (MQ=255)
                                                                                                                      gATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTACTCGGGGCAGTTTCAACCATACTCATTg                           <  1:1042951/151‑1 (MQ=255)
                                                                                                                      gATTTACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTACTCGGGGCAGTTTCAACCATACTCATTg                           >  2:784575/1‑151 (MQ=255)
                                                                                                                         ttACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTACTCGGGGCAGTTTCAACCATACTCATTGCa                         >  2:1149818/1‑150 (MQ=255)
                                                                                                                          tACCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTACTCGGGGCAGTTTCAACCATACTCATTGCAAt                       >  1:714253/1‑151 (MQ=255)
                                                                                                                           aCCGCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTACTCGGGGCAGTTTCAACCATACTCATTGCAAt                       <  2:1292343/150‑1 (MQ=255)
                                                                                                                              gCAGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTACTCGGGGCAGTTTCAACCATACTCATTGCAATCttt                   <  2:1032878/151‑1 (MQ=255)
                                                                                                                                aGCCGGACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTACTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTc                 <  1:394155/151‑1 (MQ=255)
                                                                                                                                     gACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTACTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCAt            <  1:566105/151‑1 (MQ=255)
                                                                                                                                     gACGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTACTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCAt            <  1:1160243/151‑1 (MQ=255)
                                                                                                                                       cGGCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTACTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCg          <  1:12404/151‑1 (MQ=255)
                                                                                                                                        ggCCCGATAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTACTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGt         <  1:1143080/151‑1 (MQ=255)
                                                                                                                                              aTAAACTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTACTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTg    >  2:158555/1‑150 (MQ=255)
                                                                                                                                                aaaCTCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTACTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGAt  <  2:1050472/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCAGCGTCGATTTACCGCAGCCGGACGGCCCGATAAACGCCGTTACCTGGTTTTTAGCGATATCCAGGTTGATGTTTTTCAGGGCATGGAATTTGCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTACTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGAT  >  NC_000913/3908100‑3908393

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