Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,508,454 Δ1 bp coding (137/879 nt) php ← phosphotriesterase homology protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,508,4510C.98.5% 267.0 / ‑2.6 65coding (140/879 nt)phpphosphotriesterase homology protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (33/31);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

TTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAA  >  NC_000913/3508318‑3508583
                                                                                                                                     |                                                                                                                                    
ttCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGcccggg                                                                                                                                 <  1:211176/139‑4 (MQ=255)
   cGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGtca                                                                                                                             <  2:1177607/139‑1 (MQ=255)
     ggAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGtcatc                                                                                                                           <  2:1008624/139‑1 (MQ=255)
                tCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCAt                                                                                                                <  2:29925/139‑1 (MQ=255)
                               gTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGaa                                                                                                 >  1:1119324/1‑139 (MQ=255)
                               gTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGaa                                                                                                 >  1:1205416/1‑139 (MQ=255)
                               gTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGaa                                                                                                 >  2:371014/1‑139 (MQ=255)
                                tACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAAc                                                                                                <  1:371014/139‑1 (MQ=255)
                                tACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAAc                                                                                                <  2:621284/139‑1 (MQ=255)
                                     gCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGccc                                                                                                                                    <  2:561840/99‑1 (MQ=255)
                                     gCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGccc                                                                                                                                    >  1:561840/1‑99 (MQ=255)
                                       caccacGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACt                                                                                         >  1:1244904/1‑139 (MQ=255)
                                          cacGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGAt                                                                                      >  1:928454/1‑139 (MQ=255)
                                          cacGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTAGCAAATGAACGCATACTGAt                                                                                      >  1:904227/1‑139 (MQ=255)
                                            cGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCa                                                                                    >  1:513448/1‑139 (MQ=255)
                                               tGATCCCCGTCTCGCGCACTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAgg                                                                                 >  2:300342/1‑139 (MQ=255)
                                                  tCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAggcgg                                                                              <  2:29030/139‑1 (MQ=255)
                                                       gTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGtcat                                                                                                                            >  1:607453/1‑88 (MQ=255)
                                                       gTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGtcat                                                                                                                            <  2:607453/88‑1 (MQ=255)
                                                        tctcGCGCATTACATCAAGCATAAATTGTGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGcccgggt                                                                                                                                <  2:926242/84‑5 (MQ=255)
                                                        tctcGCGCATTACATCAAGCATAAATTGTGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCgggt                                                                                                                                >  1:926242/1‑83 (MQ=255)
                                                        tctcGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTcc                                                                        >  1:218395/1‑139 (MQ=255)
                                                            gcgcATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAg                                                                                                                         <  2:1283928/86‑1 (MQ=255)
                                                            gcgcATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAg                                                                                                                         >  1:1283928/1‑86 (MQ=255)
                                                                  tACATCACGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGtttt                                                              >  2:1005576/1‑139 (MQ=255)
                                                                  tACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCAt                                                                                                                >  2:229105/1‑89 (MQ=255)
                                                                  tACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCAt                                                                                                                <  1:229105/89‑1 (MQ=255)
                                                                    cATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTa                                                            <  2:1244904/139‑1 (MQ=255)
                                                                     aTCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTaa                                                           <  2:1205416/139‑1 (MQ=255)
                                                                         aGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGcc                                                       <  1:1173986/139‑1 (MQ=255)
                                                                          gCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGtt                                                                                                                   <  1:224148/78‑1 (MQ=255)
                                                                          gCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGtt                                                                                                                   >  2:224148/1‑78 (MQ=255)
                                                                          gCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAAc                                                                                                <  2:1343103/97‑1 (MQ=255)
                                                                          gCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAAc                                                                                                >  1:1343103/1‑97 (MQ=255)
                                                                          gCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCg                                                      >  2:1136670/1‑139 (MQ=255)
                                                                            aTAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGga                                                    >  2:470920/1‑139 (MQ=255)
                                                                            aTAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGga                                                    >  1:1231/1‑139 (MQ=255)
                                                                              aaaTTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgaga                                                  >  2:772286/1‑139 (MQ=255)
                                                                              aaaTTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgaga                                                  <  1:1136670/139‑1 (MQ=255)
                                                                              aaaTTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgaga                                                  >  1:350464/1‑139 (MQ=255)
                                                                               aaTTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgagag                                                 >  1:302747/1‑139 (MQ=255)
                                                                                        ttgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATg                                        <  2:350464/139‑1 (MQ=255)
                                                                                         tgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATg                                                                                                   <  2:1023108/79‑1 (MQ=255)
                                                                                         tgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGc                                       >  1:642295/1‑139 (MQ=255)
                                                                                         tgcgcCCCATGTAACGGTTAGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATg                                                                                                   >  1:1023108/1‑79 (MQ=255)
                                                                                           cgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACgt                                                                    <  2:469590/108‑1 (MQ=255)
                                                                                           cgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACgt                                                                    >  1:469590/1‑108 (MQ=255)
                                                                                           cgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAg                                     <  2:302747/139‑1 (MQ=255)
                                                                                           cgcCCCATGTAACGGTTGGTCATCTCAATCACATTACACACGcccgg                                                                                                                                  <  2:966664/47‑3 (MQ=255)
                                                                                           cgcCCCATGTAACGGTTGGTCATCTCAATCACATTACACACGcccgg                                                                                                                                  >  1:966664/1‑45 (MQ=255)
                                                                                            gcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGa                                    <  2:1231/139‑1 (MQ=255)
                                                                                                   gTAACGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTAGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGATATCAATATGCAGATGCTCAt                             <  2:904227/139‑1 (MQ=255)
                                                                                                      aCGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATggg                          <  1:772286/139‑1 (MQ=255)
                                                                                                      aCGGTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATggg                          <  1:1514960/139‑1 (MQ=255)
                                                                                                        ggTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCg                        <  1:1005576/139‑1 (MQ=255)
                                                                                                        ggTTGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCg                        <  1:300342/139‑1 (MQ=255)
                                                                                                          ttGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                     <  2:1365726/60‑1 (MQ=255)
                                                                                                          ttGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                     >  1:1164184/1‑60 (MQ=255)
                                                                                                          ttGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                     >  1:1365726/1‑60 (MQ=255)
                                                                                                          ttGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                     <  2:1164184/60‑1 (MQ=255)
                                                                                                          ttGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgagaga                                                <  1:839569/113‑1 (MQ=255)
                                                                                                          ttGGTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgagaga                                                >  2:839569/1‑113 (MQ=255)
                                                                                                            ggTCATCTCAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAggg                    >  2:36093/1‑139 (MQ=255)
                                                                                                                  ctcAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTcc                                                                                                            <  2:1110683/45‑1 (MQ=38)
                                                                                                                  ctcAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTcc                                                                                                            >  1:1110683/1‑45 (MQ=38)
                                                                                                                   tcAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGAt                                               >  2:389635/1‑105 (MQ=255)
                                                                                                                   tcAATCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGAt                                               <  1:389635/105‑1 (MQ=255)
                                                                                                                       tCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATa                                                                                           <  1:37899/57‑1 (MQ=255)
                                                                                                                       tCACATTACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATa                                                                                           >  2:37899/1‑57 (MQ=255)
                                                                                                                            ttACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATc    >  1:1353968/1‑139 (MQ=255)
                                                                                                                             tACGCACGCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCa   >  1:670003/1‑139 (MQ=255)
                                                                                                                              acgcacgcCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCaa  <  1:1464597/139‑1 (MQ=255)
                                                                                                                              acgcacgcCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCaa  <  2:642295/139‑1 (MQ=255)
                                                                                                                              acgcacgcCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCaa  <  2:670003/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                    
TTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAA  >  NC_000913/3508318‑3508583

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