Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,876,944 G→A A122T (GCC→ACC)  iap → aminopeptidase in alkaline phosphatase isozyme conversion

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,876,9440GA98.9% 300.1 / NA 94A122T (GCC→ACC) iapaminopeptidase in alkaline phosphatase isozyme conversion
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (47/46);  total (47/47)
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

ATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTACCGAGTATGGT  >  NC_000913/2876808‑2877078
                                                                                                                                        |                                                                                                                                      
aTTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACAcc                                                                                                                                      <  1:1246355/139‑1 (MQ=255)
            gTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGAt                                                                                                                          >  2:480387/1‑139 (MQ=255)
                     taCCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGcc                                                                                                                 <  2:705115/139‑1 (MQ=255)
                       ccgcccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGa                                                                                                               >  2:1373089/1‑139 (MQ=255)
                         gcccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCCTTATGGCGCATCTGGATACTTCCACCCCGCTGAGCGATGCTGACGCCGAtt                                                                                                             >  1:359235/1‑138 (MQ=255)
                         gcccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATg                                                                                                             >  1:1433188/1‑139 (MQ=255)
                         gcccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATg                                                                                                             >  1:356444/1‑139 (MQ=255)
                         gcccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATg                                                                                                             >  2:1016549/1‑139 (MQ=255)
                         gcccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATg                                                                                                             >  1:1032578/1‑139 (MQ=255)
                          cccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGTTCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGc                                                                                                            >  2:472532/1‑139 (MQ=255)
                          cccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGca                                                                                                                       >  2:1042681/1‑127 (MQ=255)
                          cccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCt                                                                                                                       <  1:1042681/128‑1 (MQ=255)
                          cccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGc                                                                                                                  <  1:282412/133‑1 (MQ=255)
                          cccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGc                                                                                                                  >  2:282412/1‑133 (MQ=255)
                          cccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGc                                                                                                            >  1:479813/1‑139 (MQ=255)
                          cccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGc                                                                                                            >  1:332365/1‑139 (MQ=255)
                          cccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGc                                                                                                            >  2:1060993/1‑139 (MQ=255)
                          cccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGc                                                                                                            >  1:591162/1‑139 (MQ=255)
                          cccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGc                                                                                                            >  2:295593/1‑139 (MQ=255)
                          cccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGc                                                                                                            >  1:1285911/1‑139 (MQ=255)
                          cccgcGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGc                                                                                                            >  2:604945/1‑139 (MQ=255)
                            cgcgATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCg                                                                                                                            <  2:448432/121‑1 (MQ=255)
                            cgcgATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCg                                                                                                                            >  1:448432/1‑121 (MQ=255)
                               gATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCATCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAAtc                                                                                                       >  1:925093/1‑139 (MQ=255)
                                  aaTCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTcg                                                                                                    <  2:764383/139‑1 (MQ=255)
                                         agagCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCt                                                                                             <  2:479813/139‑1 (MQ=255)
                                          gagCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTg                                                                                            <  1:1361077/139‑1 (MQ=255)
                                           agCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGa                                                                                           <  1:1016549/139‑1 (MQ=255)
                                             cTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACg                                                                                         >  2:1290759/1‑139 (MQ=255)
                                             cTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACg                                                                                         <  2:1285911/139‑1 (MQ=255)
                                               ggCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGtt                                                                                       >  1:838518/1‑139 (MQ=255)
                                               ggCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGtt                                                                                       >  2:736128/1‑139 (MQ=255)
                                               ggCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCTCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCt                                                                                             <  2:584244/132‑1 (MQ=255)
                                               ggCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCTCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCt                                                                                             >  1:584244/1‑132 (MQ=255)
                                                  aTAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACa                                                                                    <  1:1373089/139‑1 (MQ=255)
                                                    aaCGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAg                                                                                  <  2:332365/139‑1 (MQ=255)
                                                     aCGTGACGGGAAGTACTGTGATTGCCGCGCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCg                                                                                                                                   <  2:978040/89‑1 (MQ=255)
                                                     aCGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCg                                                                                                                                   >  1:978040/1‑89 (MQ=255)
                                                      cGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACg                                                                                         <  1:1203633/130‑1 (MQ=255)
                                                      cGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACg                                                                                         >  2:1203633/1‑130 (MQ=255)
                                                      cGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGa                                                                                <  1:604945/139‑1 (MQ=255)
                                                       gTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGaa                                                                               <  1:1060993/139‑1 (MQ=255)
                                                          aCGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATgg                                                                            >  2:1608738/1‑139 (MQ=255)
                                                                 gTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCggg                                                                                               >  1:1031297/1‑113 (MQ=255)
                                                                 gTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCggg                                                                                               <  2:1031297/113‑1 (MQ=255)
                                                                 gTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATg                                                                             <  1:1321817/131‑1 (MQ=255)
                                                                 gTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATg                                                                             >  2:1321817/1‑131 (MQ=255)
                                                                 gTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATaa                                                                     <  2:838518/139‑1 (MQ=255)
                                                                     ggTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATTTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAATGAATGGATGATAAcgcc                                                                 <  2:359235/139‑1 (MQ=255)
                                                                     ggTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAAcgcc                                                                 <  2:356444/139‑1 (MQ=255)
                                                                        gATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCa                                                              >  1:1490169/1‑139 (MQ=255)
                                                                        gATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTAACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCa                                                              >  1:1546551/1‑139 (MQ=255)
                                                                         aTTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAg                                                             >  2:1371625/1‑139 (MQ=255)
                                                                                  cATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGgtg                                                    <  1:480387/139‑1 (MQ=255)
                                                                                    tGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTc                                                  <  1:295593/139‑1 (MQ=255)
                                                                                       aGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATg                                               <  2:591162/139‑1 (MQ=255)
                                                                                         gtaaaGCGCCGCAGCAGATCAGCCTTATGGCTCATCTGGATACTAACGCCCCGCTGAGCGATGCTGACGCCGATGCCAGGCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCt                                             <  1:725748/137‑1 (MQ=255)
                                                                                                  cgcAGCAGATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGGATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGc                                    >  2:762071/1‑139 (MQ=255)
                                                                                                    cagcagATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAg                                  <  2:925093/139‑1 (MQ=255)
                                                                                                      gcagATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGaa                                <  1:1290759/139‑1 (MQ=255)
                                                                                                      gcagATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGaa                                <  2:330486/139‑1 (MQ=255)
                                                                                                         gATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCa                                                 <  2:386342/119‑1 (MQ=255)
                                                                                                         gATCATCATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCa                                                 >  1:386342/1‑119 (MQ=255)
                                                                                                          atcatcATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATaa                                                                     >  1:452416/1‑98 (MQ=255)
                                                                                                          atcatcATTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATaa                                                                     <  2:452416/98‑1 (MQ=255)
                                                                                                            catcatTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTg                          <  2:1374575/139‑1 (MQ=255)
                                                                                                               catTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGc                                                                                                                  >  2:1278602/1‑48 (MQ=255)
                                                                                                               catTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGc                                                                                                                  <  1:1278602/48‑1 (MQ=255)
                                                                                                               catTATGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGaaa                       <  1:472532/139‑1 (MQ=255)
                                                                                                                   aTGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATa                   <  1:1608738/139‑1 (MQ=255)
                                                                                                                   aTGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATa                   <  1:762071/139‑1 (MQ=255)
                                                                                                                    tGGCGCATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATAc                  >  2:613261/1‑139 (MQ=255)
                                                                                                                      gcgcATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGc                >  2:487066/1‑139 (MQ=255)
                                                                                                                       cgcATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGa                                                                                >  2:1363498/1‑74 (MQ=255)
                                                                                                                       cgcATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGa                                                                                <  1:1363498/74‑1 (MQ=255)
                                                                                                                       cgcATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATg                                               <  2:1188662/107‑1 (MQ=255)
                                                                                                                       cgcATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATg                                               >  1:1188662/1‑107 (MQ=255)
                                                                                                                       cgcATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGc                <  2:850537/138‑1 (MQ=255)
                                                                                                                       cgcATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGc                >  1:850537/1‑138 (MQ=255)
                                                                                                                         cATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGc                                                                                              <  1:794861/58‑1 (MQ=255)
                                                                                                                         cATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGc                                                                                              >  2:794861/1‑58 (MQ=255)
                                                                                                                         cATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAAcgccgc                                                               >  2:284747/1‑89 (MQ=255)
                                                                                                                         cATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAAcgccgc                                                               <  1:284747/89‑1 (MQ=255)
                                                                                                                         cATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTg                          <  2:1323268/126‑1 (MQ=255)
                                                                                                                         cATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTg                          >  1:1323268/1‑126 (MQ=255)
                                                                                                                         cATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGc                <  1:787923/136‑1 (MQ=255)
                                                                                                                         cATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGc                >  2:787923/1‑136 (MQ=255)
                                                                                                                         cATCTGGATACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTa             >  2:312382/1‑139 (MQ=255)
                                                                                                                                aTACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTACCGAGTa      <  1:312382/139‑1 (MQ=255)
                                                                                                                                aTACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTACCGAGTa      <  1:1371625/139‑1 (MQ=255)
                                                                                                                                aTACTTACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTACCGAGTa      <  1:736128/139‑1 (MQ=255)
                                                                                                                                    ttACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCt                                                                                             <  1:483357/48‑1 (MQ=255)
                                                                                                                                    ttACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCt                                                                                             >  2:483357/1‑48 (MQ=255)
                                                                                                                                    ttACACCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTACCGAGTATGGt  >  1:141758/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                      
ATTTAATAGTCGGTATATTTATACCGCCCGCGATAATCGTAAGAGCTGGCATAACGTGACGGGAAGTACGGTGATTGCCGCTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTAGGTGTCATGCTGGAATTGGCAGAACGCCTGAAAAATACGCCTACCGAGTATGGT  >  NC_000913/2876808‑2877078

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