Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,724,095 C→A 7.3% intergenic (+13/‑33) pssA → / → yfiM phosphatidylserine synthase/protein YfiM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,724,0950CA 7.3% 170.3 / 3.8 67intergenic (+13/‑33)pssA/yfiMphosphatidylserine synthase/protein YfiM
Reads supporting (aligned to +/- strand):  ref base C (28/33);  new base A (2/3);  total (30/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01

CACCATCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGAT  >  NC_000913/2723962‑2724237
                                                                                                                                     |                                                                                                                                                
caccaTCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACggg                                                                                                                                         >  1:111228/1‑143 (MQ=255)
  ccaTCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTACAGGACGGGgt                                                                                                                                       <  1:854264/143‑1 (MQ=255)
    agcGTTAGGCACTATCGCGATCTGCGTAGTGTTGCCGTTTATGCGGTGGTGGTTCGTTGTCTCGTGGGCCGTTTGCGCCGTATCCGCATCGACCGATTAGTTGGGCGCATCCTGTAATCACAACCCCGTACAGGACGGGgttg                                                                                                                                     <  2:523484/141‑2 (MQ=255)
              gctctCGCGATCTGCAAAGTATTGCCGGTTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACATCCCCGTCCTGTACGGGGTTTGTTTTTTGGa                                                                                                                         <  2:885958/139‑1 (MQ=255)
              aCTATAGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAAGTCATCCGCCGTTTGCGCCGTATCCGCATCGATCGATTAATTAGCCGCATAATGCCCCTACATCCGCGTACTGTACGGGGTTTGTTTTTTGGa                                                                                                                         <  2:231173/143‑1 (MQ=255)
               cTATCGCGCTCGTCAAAGGATTGCCGATTATCCGGTGAAGGTGCGTAAACTCAGCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAg                                                                                                                        <  1:389983/143‑1 (MQ=255)
               cTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAg                                                                                                                        <  1:761135/143‑1 (MQ=255)
                tATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTGGTCTTTTTgggg                                                                                                                        >  1:887872/1‑141 (MQ=255)
                tATCGCGATCGGCAAAGGATGGCCGATGATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGCCCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAgg                                                                                                                       <  1:64356/143‑1 (MQ=255)
                     cgATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACg                                                                                                                  >  2:307253/1‑143 (MQ=255)
                      gATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACCGGGTCTGTTTGTTGGAGGCAACgg                                                                                                                 >  2:147632/1‑142 (MQ=255)
                          tGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTTGGAGGCCACGtgtt                                                                                                              >  1:415299/1‑143 (MQ=255)
                            cAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGCAGGCCACGTTTTGCg                                                                                                           >  1:1032564/1‑143 (MQ=255)
                            cAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCCGTACGGGGGCCGTCGGTGGGGGGCCACGCTTTGCg                                                                                                           >  2:693205/1‑143 (MQ=255)
                             aaaGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCGACGTTTTGCGt                                                                                                          >  1:406486/1‑143 (MQ=255)
                               aggaTTGCCGGGTATCCGGTGAAGGTTCGTAAACGCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTAt                                                                                                        <  1:137333/140‑1 (MQ=255)
                                    tGCCGATTATCCGGTGAAGGTTCGTAGACTCATCCGGCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGACACGTTTTGCGTATTCttt                                                                                                   <  2:181302/143‑1 (MQ=255)
                                     gCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCtttt                                                                                                  >  1:592685/1‑143 (MQ=255)
                                      tcggttATCCGGTGAAGGGTCGTAAACGCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCttttt                                                                                                 <  1:370431/139‑1 (MQ=255)
                                      cggatgatCCGGTGAAGGTTCGGACACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACACCCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCttttt                                                                                                 <  2:379636/141‑1 (MQ=255)
                                       cGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCTTCTACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTg                                                                                                <  1:491099/143‑1 (MQ=255)
                                       cGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTg                                                                                                >  1:567555/1‑143 (MQ=255)
                                       cGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTACGCCGTATCCGCATCGACCGACTAACTAGCCGCCGACTGTAATCACAACCCCGTCCAGTACGGGGATTGGTATTGGAGGGCCACGATCTGCGTAGTCttgttg                                                                                                >  2:855532/1‑143 (MQ=255)
                                           caTCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTg                                                                                            <  1:225624/142‑1 (MQ=255)
                                               cGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCg                                                                                        <  1:139401/143‑1 (MQ=255)
                                               cGGGGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCg                                                                                        <  2:191389/143‑1 (MQ=255)
                                                    aaGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTAtt                                                                                   <  2:592685/143‑1 (MQ=255)
                                                    aaGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTCTGTTTTTTGGAGGCCACGTTTTGCGTAGTCTTTTTGTCTGTTCGCTAtt                                                                                   >  1:81102/1‑143 (MQ=255)
                                                               aCTCATCCGCCGTTTGCGCCGTATCCGCATCGCCCGATTAATTAGCCGCATACTGTAATCACAACCCCGTCCTGTACGGGGTTTGTCTTTTGGAGGCCACGTTTTGCGTATGCTTTTTCTATGTTCGCTATTGTTACTTTCTg                                                                        >  2:560713/1‑143 (MQ=255)
                                                               aCTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACACCATCAAGGACAGTGTGTGTTTTTTGGAGGCCACGATTTGCGTATTCTTTTTGTCTATTCTCTATTGCTACTTTCTg                                                                        >  2:869079/1‑143 (MQ=255)
                                                                 tCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGCTTTGCGTATTCTTTTTGTCTGCTCGCTATTGTTGCTGTCTGGa                                                                      >  2:848119/1‑143 (MQ=255)
                                                                  cgtctgccgTTTGTTCCGTATCTGCTTTGTCCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGTGGTTTGTTTTTTGGAGGCTACGTTTTTCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGAt                                                                     <  2:82046/138‑1 (MQ=255)
                                                                  cATCCGCCGGTGGCGCCGTATCCGCATCGACCGATAAATCAGCCGCTTCCTGTAGTCACAACCCCGTCCTGTACGGGGTTTGTTTTGTGGAGGCCACGTGTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGAt                                                                     <  1:36762/143‑1 (MQ=255)
                                                                    tCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACACCCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGc                                                                   <  2:52082/143‑1 (MQ=255)
                                                                    tCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGc                                                                   >  1:1013323/1‑143 (MQ=255)
                                                                     cctccgTGTGCGCCGGCTCCGCGGCGACCGATTTCTTTGCCGCGTCCTGTGTGCACATTCCCGTCCTGTACGTGGTTTGTTTTTTGTAGGCTACGTTTTTCGTAGCCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCa                                                                  <  2:763371/140‑1 (MQ=255)
                                                                      cgccgTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAg                                                                 <  1:958702/143‑1 (MQ=255)
                                                                       gccgTTTGCGCGGTATCCGCATCGACCGATTATTTAGCCGCATCCTGTAGTCACAACGCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACTTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGc                                                                <  2:788794/143‑1 (MQ=255)
                                                                          gTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCat                                                             <  1:827913/143‑1 (MQ=255)
                                                                            tcgcgcCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCatat                                                           <  1:696813/141‑1 (MQ=255)
                                                                                  cGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTaa                                                     <  1:445954/143‑1 (MQ=255)
                                                                                        cGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAg                                               >  1:366178/1‑143 (MQ=255)
                                                                                           aTCGACCGATTAGTTCGCCGCATCCTGTAATTACAACCGCGTTCTGTACGGGGTTTGTTGTTCGGAGGCCACTTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTg                                            <  2:742229/143‑1 (MQ=255)
                                                                                           aTCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTg                                            <  2:567555/143‑1 (MQ=255)
                                                                                               aCCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGc                                        <  2:366178/143‑1 (MQ=255)
                                                                                                  gATTATTGAGCCGCATCCTGTAATCACACCCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCggg                                     <  2:503583/143‑1 (MQ=255)
                                                                                                   attaattaGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTTGTGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGc                                    >  1:702911/1‑143 (MQ=255)
                                                                                                    ttatgtggCCGCAGCCTGTAATCACCACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGGTTTGCGTATTCTTTTTGTCTGTCCGCTATGGTTACCTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCa                                   <  2:188429/136‑1 (MQ=255)
                                                                                                       gttgggcGCATCCTGTGAGCCCATCCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGa                                <  1:108811/137‑1 (MQ=255)
                                                                                                        ttaGCCGCATCCTGTAATCACAACCCCGTCCTGTCCGGGGTTTGTTTTTTGGAGGCCACGTCTTGCGTATTCTGTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAAGGATAGCTGGAGCGGGTAGGAt                               >  2:317405/1‑143 (MQ=255)
                                                                                                           gCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATaaa                            >  2:398162/1‑143 (MQ=255)
                                                                                                            ccGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGGTGCAGCCATATGGATAACGATAGCTGGAGCGGGCAGGATAAAg                           >  1:499700/1‑143 (MQ=255)
                                                                                                            ccGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAg                           >  2:455016/1‑143 (MQ=255)
                                                                                                              gcctCCCGTAATCAAAACCCCGTCCGGTACGGGGTTTGTTTTTTGGAGGCCACGCTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCt                         <  1:31038/140‑1 (MQ=255)
                                                                                                              gCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCt                         <  1:455098/143‑1 (MQ=255)
                                                                                                               cATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTc                        <  2:486772/143‑1 (MQ=255)
                                                                                                                 tcgtGTAATCACAGCCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCaa                      <  1:693205/140‑1 (MQ=255)
                                                                                                                    tGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAacac                   >  1:496886/1‑143 (MQ=255)
                                                                                                                    tGTAATCACAACCCCGTCCTGTACGGGAATTGTTATTTGAAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCCGCCATATGGCTAACGATAGCAGGAGCGGGATGGATAAAGCTCAacac                   >  1:152862/1‑143 (MQ=255)
                                                                                                                      aaaaacaaaCCCCGTACAGGACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACtt                 >  1:548283/8‑143 (MQ=255)
                                                                                                                        aTCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGCGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCGACACTGTa               >  1:757352/1‑143 (MQ=255)
                                                                                                                           acaTCCCCGGCCTGTACGGGGTTGGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCg            <  2:387874/143‑1 (MQ=255)
                                                                                                                             aaCCCCGTACAGGACGGGGTTTGTTTTTTGGAGGCCAAAAAGTCAGTAACCTGGTTGTCTGTTCGCTATTGTTCCTTTCTGGATGCAGCCATATGGCTACCGACAGCAGCAGCGGGCACGATAAACATCAACACATTATCGca          >  2:228672/1‑142 (MQ=255)
                                                                                                                               ccccGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTc        >  1:795745/1‑143 (MQ=255)
                                                                                                                                cccGTCCTGTACGGGGTTTGATTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAAAACTTTATCGCCTCg       >  2:85317/1‑143 (MQ=255)
                                                                                                                                 ccGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGCTGCAGCCATGTGGTTAACGCTAGCTGGCGCGGGCAGGATCAAGCTCTACACTTTATCGCGTCgt      >  1:52467/1‑142 (MQ=255)
                                                                                                                                    tCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGa   <  1:455016/143‑1 (MQ=255)
                                                                                                                                     ccTGGAGGGGGGTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGAt  <  1:536027/143‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                                
CACCATCGTTAAGCACTATCGCGATCTGCAAAGTATTGCCGATTATCCGGTGAAGGTTCGTAAACTCATCCGCCGTTTGCGCCGTATCCGCATCGACCGATTAATTAGCCGCATCCTGTAATCACAACCCCGTCCTGTACGGGGTTTGTTTTTTGGAGGCCACGTTTTGCGTATTCTTTTTGTCTGTTCGCTATTGTTACTTTCTGGATGCAGCCATATGGCTAACGATAGCTGGAGCGGGCAGGATAAAGCTCAACACTTTATCGCCTCGGCGAT  >  NC_000913/2723962‑2724237

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