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)
   cGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGtca                                                                                                                             <  2:1177607/139‑1 (MQ=255)
     ggAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGtcatc                                                                                                                           <  2:1008624/139‑1 (MQ=255)
                tCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCAt                                                                                                                <  2:29925/139‑1 (MQ=255)
                               gTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGaa                                                                                                 >  1:1119324/1‑139 (MQ=255)
                               gTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGaa                                                                                                 >  1:1205416/1‑139 (MQ=255)
                               gTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGaa                                                                                                 >  2:371014/1‑139 (MQ=255)
                                tACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAAc                                                                                                <  1:371014/139‑1 (MQ=255)
                                tACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAAc                                                                                                <  2:621284/139‑1 (MQ=255)
                                     gCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGccc                                                                                                                                    <  2:561840/99‑1 (MQ=255)
                                     gCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGccc                                                                                                                                    >  1:561840/1‑99 (MQ=255)
                                       caccacGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACt                                                                                         >  1:1244904/1‑139 (MQ=255)
                                          cacGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGAt                                                                                      >  1:928454/1‑139 (MQ=255)
                                          cacGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTAGCAAATGAACGCATACTGAt                                                                                      >  1:904227/1‑139 (MQ=255)
                                            cGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCa                                                                                    >  1:513448/1‑139 (MQ=255)
                                               tGATCCCCGTCTCGCGCACTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAgg                                                                                 >  2:300342/1‑139 (MQ=255)
                                                  tCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAggcgg                                                                              <  2:29030/139‑1 (MQ=255)
                                                       gTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGtcat                                                                                                                            >  1:607453/1‑88 (MQ=255)
                                                       gTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGtcat                                                                                                                            <  2:607453/88‑1 (MQ=255)
                                                        tctcGCGCATTACATCAAGCATAAATTGTGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGcccgggt                                                                                                                                <  2:926242/84‑5 (MQ=255)
                                                        tctcGCGCATTACATCAAGCATAAATTGTGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCgggt                                                                                                                                >  1:926242/1‑83 (MQ=255)
                                                        tctcGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTcc                                                                        >  1:218395/1‑139 (MQ=255)
                                                            gcgcATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAg                                                                                                                         <  2:1283928/86‑1 (MQ=255)
                                                            gcgcATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAg                                                                                                                         >  1:1283928/1‑86 (MQ=255)
                                                                  tACATCACGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGtttt                                                              >  2:1005576/1‑139 (MQ=255)
                                                                  tACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCAt                                                                                                                >  2:229105/1‑89 (MQ=255)
                                                                  tACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCAt                                                                                                                <  1:229105/89‑1 (MQ=255)
                                                                    cATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTa                                                            <  2:1244904/139‑1 (MQ=255)
                                                                     aTCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTaa                                                           <  2:1205416/139‑1 (MQ=255)
                                                                         aGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGcc                                                       <  1:1173986/139‑1 (MQ=255)
                                                                          gCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGtt                                                                                                                   <  1:224148/78‑1 (MQ=255)
                                                                          gCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGtt                                                                                                                   >  2:224148/1‑78 (MQ=255)
                                                                          gCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAAc                                                                                                <  2:1343103/97‑1 (MQ=255)
                                                                          gCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAAc                                                                                                >  1:1343103/1‑97 (MQ=255)
                                                                          gCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCg                                                      >  2:1136670/1‑139 (MQ=255)
                                                                            aTAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGga                                                    >  2:470920/1‑139 (MQ=255)
                                                                            aTAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGga                                                    >  1:1231/1‑139 (MQ=255)
                                                                              aaaTTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgaga                                                  >  2:772286/1‑139 (MQ=255)
                                                                              aaaTTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgaga                                                  <  1:1136670/139‑1 (MQ=255)
                                                                              aaaTTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgaga                                                  >  1:350464/1‑139 (MQ=255)
                                                                               aaTTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgagag                                                 >  1:302747/1‑139 (MQ=255)
                                                                                        ttgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATg                                        <  2:350464/139‑1 (MQ=255)
                                                                                         tgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATg                                                                                                   <  2:1023108/79‑1 (MQ=255)
                                                                                         tgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGc                                       >  1:642295/1‑139 (MQ=255)
                                                                                         tgcgcCCCATGTAACGGTTAGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATg                                                                                                   >  1:1023108/1‑79 (MQ=255)
                                                                                           cgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACgt                                                                    <  2:469590/108‑1 (MQ=255)
                                                                                           cgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACgt                                                                    >  1:469590/1‑108 (MQ=255)
                                                                                           cgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAg                                     <  2:302747/139‑1 (MQ=255)
                                                                                           cgcCCCATGTAACGGTTGGTCATCTCAATCACATTACACACGcccgg                                                                                                                                  <  2:966664/47‑3 (MQ=255)
                                                                                           cgcCCCATGTAACGGTTGGTCATCTCAATCACATTACACACGcccgg                                                                                                                                  >  1:966664/1‑45 (MQ=255)
                                                                                            gcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGa                                    <  2:1231/139‑1 (MQ=255)
                                                                                                   gTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTAGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGATATCAATATGCAGATGCTCAt                             <  2:904227/139‑1 (MQ=255)
                                                                                                      aCGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATggg                          <  1:772286/139‑1 (MQ=255)
                                                                                                      aCGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATggg                          <  1:1514960/139‑1 (MQ=255)
                                                                                                        ggTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCg                        <  1:1005576/139‑1 (MQ=255)
                                                                                                        ggTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCg                        <  1:300342/139‑1 (MQ=255)
                                                                                                          ttGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                     <  2:1365726/60‑1 (MQ=255)
                                                                                                          ttGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                     >  1:1164184/1‑60 (MQ=255)
                                                                                                          ttGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                     >  1:1365726/1‑60 (MQ=255)
                                                                                                          ttGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                     <  2:1164184/60‑1 (MQ=255)
                                                                                                          ttGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgagaga                                                <  1:839569/113‑1 (MQ=255)
                                                                                                          ttGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgagaga                                                >  2:839569/1‑113 (MQ=255)
                                                                                                            ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAggg                    >  2:36093/1‑139 (MQ=255)
                                                                                                                  ctcAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTcc                                                                                                            <  2:1110683/45‑1 (MQ=38)
                                                                                                                  ctcAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTcc                                                                                                            >  1:1110683/1‑45 (MQ=38)
                                                                                                                   tcAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGAt                                               >  2:389635/1‑105 (MQ=255)
                                                                                                                   tcAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGAt                                               <  1:389635/105‑1 (MQ=255)
                                                                                                                       tCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATa                                                                                           <  1:37899/57‑1 (MQ=255)
                                                                                                                       tCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATa                                                                                           >  2:37899/1‑57 (MQ=255)
                                                                                                                            ttACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATc    >  1:1353968/1‑139 (MQ=255)
                                                                                                                             tACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCa   >  1:670003/1‑139 (MQ=255)
                                                                                                                              acgcacg‑cCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCaa  <  1:1464597/139‑1 (MQ=255)
                                                                                                                              acgcacg‑cCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCaa  <  2:642295/139‑1 (MQ=255)
                                                                                                                              acgcacg‑cCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCaa  <  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: ‑