Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,902,100 G→A 7.2% A69T (GCC→ACC)  ygcE → putative sugar kinase YgcE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,902,1000GA 7.2% 164.1 / 2.9 65A69T (GCC→ACC) ygcEputative sugar kinase YgcE
Reads supporting (aligned to +/- strand):  ref base G (33/26);  new base A (3/2);  total (36/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.66e-01

CTGGAAGGTAACGTGGTTTGCGAAGGCAAAGGCTTATTACAGCCGATGCACACGCCAGATGCCGATACTGCAGAACATCCTGACGACGATTTATGGGCATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACA  >  NC_000913/2901962‑2902236
                                                                                                                                          |                                                                                                                                        
cTGGAAGGTAACGTGGTTTGCGAAGGCTAAGGCTTTTTACAGCCGGTGCACACGCCAGATGCCGATACTGCAGAACATCCTGACGACGATTTATGGGCATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCgg                                                                                                                                      <  2:122708/143‑1 (MQ=255)
cTGGAAGGTAACGTGGTTTGCGAAGGCAAAGGCTTATTACAGCCGATGCACACGCCAGATGCCGATACTGCAGAACATCCTGACGACGATTTATGGGCATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCgg                                                                                                                                      <  2:143379/143‑1 (MQ=255)
   gAAGGTAACGTGGTTTGCGAAGGCAAAGGCTTATTACAGCCGATGCACACGCCAGATGCCGATACTGCAGAACATCCTGACGACGATTTATGGGCATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGaa                                                                                                                                   >  1:261409/1‑143 (MQ=255)
   gAAGGTAACGTGGTTTGCGAAGGCAAAGGCTTATTACAGCCGATGCACACGCCAGATGCCGATACTGCAGAACATCCTGACGACGATTTATGGGCATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGaa                                                                                                                                   <  1:742259/143‑1 (MQ=255)
         aaCGTGGTTTGCGAAGGCAAAGGCTTATTACAGCCGATGCACACGCCAGATGCCGATACTGCAGAACATCCTGACGACGATTTATGGGCATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGaagctaga                                                                                                                             >  2:1005353/1‑137 (MQ=255)
               gTTTGCGAAGGCAAAGGCTTATTACAGCCGATGCACACGCCAGATGCCGATACTGCAGAACATCCTGACGACGATTTATGGGCATCATTATGTTTTGCCGGCCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGctat                                                                                                                       >  2:858146/1‑143 (MQ=255)
                     gAAGGCAAAGGCTTATTACAGCCGATGCACACGCCAGATGCCGATACTGCAGAACATCCTGACGACGATTTCTGGGCATCATTATGTTTTGATGCTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCgg                                                                                                                 >  1:567040/1‑143 (MQ=255)
                     gAAGGCAAAGGCTTATTACAGCCGATGCACACGCCAGATGCCGATACTGCAGAACATCCTGACGACGATTTATGGGCATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCgg                                                                                                                 <  1:920018/143‑1 (MQ=255)
                      aaGGCAAAGGCTTATTACAGCCGATGCACACGCCAGATGCCGATACTGCAGAACATCCTGACGACGATTTATGGGCATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGt                                                                                                                >  2:904730/1‑143 (MQ=255)
                          cAAAGGCTTATTACAGCCGATGCACACGCCAGATGCCGATACTGCAGAACATCGGGACGACGATTTATGGGCATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTg                                                                                                            <  1:17595/143‑1 (MQ=255)
                            aaGGCTTATTACAGCCGATGCACACGCCAGATGCCGATACTGCAGAACATCCTGACGACGCTTTATGGGCCTCATTCTGTTTTGCCGGTCACGATTTGCTGAGTCAGTTTGCCGGGCATAAAGACGACATTGTCGGTATTCGt                                                                                                          >  2:1089475/1‑143 (MQ=255)
                               gCTTAGTGCAGCCGATGCACACGCCAGTTGCCGATACTGCAGACCATCCTGACGACGATTTATGGGCATCATTATGTTTTGCGGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTg                                                                                                       <  2:114903/143‑1 (MQ=255)
                                 ttattaCAGCCGATGCACACGCCAGATGCCGATACTGCAGAACATCCTGACGACGATTTATGGGCATCATTATGTTTTGCCGGTCACGATTTGATGCGTCAGTTTGCCGGGAATAAAGAAGATATTGGCGGTATTGGTCgggg                                                                                                     >  2:930683/1‑143 (MQ=255)
                                  tatgaCAGCGGATGCACACGCGAGATGCCGATACTGCAGAACATCCTGACGACGATTTATGGGCATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGa                                                                                                    <  1:427485/143‑1 (MQ=255)
                                             aTGCACACGCCAGATGCCGATACTGGAGAATCTCCCGACGACGATTTATGGGCATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTtg                                                                                         <  1:464070/143‑1 (MQ=255)
                                                 acacGCCAGATGCCGATACTGCAGAACATCCTGACGACGATTTATGGGCATCATTCTGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGGCGGTATTGGTCTGGGCTCCATCCGTtgctgc                                                                                     >  2:350689/1‑143 (MQ=255)
                                                 acacGCCAGATGCCGATACTGCAGAACATCCTGACGACGATTTATGGGCATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTtgctgc                                                                                     >  1:347203/1‑143 (MQ=255)
                                                  cacGCCAGATGCCGATACTGCAGAACATCCTGACGACGATTTATGGGCATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGCATAAAGAAGACATTGTCGGTATTGGTCGGGGATCCATCCGTTGCTGcc                                                                                    >  2:256791/1‑143 (MQ=255)
                                                   acGCCAGATGCCGATACTGCAGAACATCCTGACGACGATTTATGGGCATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGACGATATGGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCg                                                                                   >  1:449994/1‑143 (MQ=255)
                                                    cGCTAGATGCCGATACTGCAGAACATCCTTACGACGATTTATGGGCATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGt                                                                                  <  1:866548/143‑1 (MQ=255)
                                                            gCCGATGGTTCAGTACATCCTGATGACGATTTATGGTGATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTAtt                                                                          <  2:741428/143‑1 (MQ=255)
                                                             ccGATACTGCAGAACATCCTGACGACGATTTATGGGCATCATTATGTTTTGCCGGTCACGCTTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCCGGGCTCCATCCGTTGCTGCCGTGCGTTAtcg                                                                         >  2:1039634/1‑141 (MQ=255)
                                                              cGATACTGCAGAACATCCTGACGACGATTCATGGGCATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGCCGGTACCAGTCTGGGAACCACCCCGTGCAGCCGTCCGTGATGGa                                                                        >  2:989846/1‑143 (MQ=255)
                                                                aTACTGCAGAACATCCTGACGACGATTTATGGGCATCATTATGTTTTGCCGGTCACGATTTGATTATTCAGTTTTCTGTGAATAAAGAAGAAATTGTCGGTATTGGTCTGGGATGCATCCGTTGCTGCCGTGCGTTATTGaaa                                                                      <  1:308966/143‑1 (MQ=255)
                                                                aTACTGCAGAACATCCTGAAGACGATTTATGGGCATCATTATGTTTTGCGGGTCACGATTTGATGATTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGaaa                                                                      <  2:232829/143‑1 (MQ=255)
                                                                     gCAGAACATCCTGACGACGATTTATGGGCATCATTATGTTTTGACGGTCGCGATTTCATGGTTCAGGTTGGCGGGAATAAAGAAGATATTGTCGGTGTTGGTCTGGGGTCGATCCGTTGATGGCGTGCGCAATTGaagccgga                                                                 >  2:153611/1‑137 (MQ=255)
                                                                      cAGATCGTCCGGACTACGATGTATGGGCATGATTATGTTTTGCCGGTCACGATTTGCTGAGTCAGTTTGCCGGGAATAAAGAAGATATGGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGAt                                                                <  2:298763/143‑1 (MQ=255)
                                                                      cAGAACATCCTGACGACGATTTATGGGCATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGAt                                                                >  2:43625/1‑143 (MQ=255)
                                                                      cAGAACATCCTGACGACGATTTATGGGCATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGAt                                                                >  2:860820/1‑143 (MQ=255)
                                                                      cAGAACATCCTGACGACGATTTATGGGCATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGAAAATGTAGCCATCGGTCTGGGATACATCCGTTGATGCCGTGCGATATTGAAAGCCGAt                                                                >  2:722293/1‑143 (MQ=255)
                                                                         aaCCTCCTGACGGCGATTTATGGGCCTCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGGc                                                             <  1:609293/143‑1 (MQ=255)
                                                                             tCCTGACGACGATTTATGGGCATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGGCACGc                                                         >  2:375745/1‑143 (MQ=255)
                                                                                  cggacgtTTGATGGCCTTCATTATGTTGTGCCGTTGACGATTTGATTAGTCAGTTTGCCGGGAATAAAGCAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGGCACGCctgct                                                    <  2:195961/141‑1 (MQ=255)
                                                                                   cgacgaTTTATGGGCATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGTCGGGAACAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGGCACGCctgctg                                                   <  1:807560/143‑1 (MQ=255)
                                                                                      cggtttGTGGGCATGATTAGGTTTTGCGGGTCACGATGTGGTGAGTCAGTTTGCCGGGAATAAGGAAGGTGTTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTAATGAAAGCCGATGGCACGCCTGCTgcgc                                                <  2:276776/140‑1 (MQ=255)
                                                                                       gaTTTATGGGCATCATTATGTTTTGCCGGTCACGCTTTGATGAATCAGTTTACCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGAACCATCAGTTGCTACCGTGCGTTATAAAAAGCCGATGGCCAGCCtgctgcaca                                               >  2:608001/1‑140 (MQ=255)
                                                                                       gaTTTATGGGCATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGGCACGCCTGCTGCGcc                                               >  1:912838/1‑143 (MQ=255)
                                                                                         tttATGGGCATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCAGATTGCACGCCTGCTGCTCCGt                                             >  1:421188/1‑143 (MQ=255)
                                                                                          ttATGGGGATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGGCACGCCTGCTGCGCCGtt                                            <  2:986656/143‑1 (MQ=255)
                                                                                          ttATGGGCATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTACCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGatccatc                                                                                              >  2:973640/1‑93 (MQ=255)
                                                                                          ttATGGGCATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTACCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGatccatc                                                                                              <  1:973640/93‑1 (MQ=255)
                                                                                          ttATGGGCATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTACCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGatccatc                                                                                              <  1:966313/93‑1 (MQ=255)
                                                                                          ttATAGCCATGATTACGTAGTGCCGGGCACAACTTGATGAGTCAGTTTGCCGGAAATAAAGAAGATATTCTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGGCACGCCTGCTGCGCCGtt                                            <  2:284767/143‑1 (MQ=255)
                                                                                          ttACGGGCATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTACCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGatccatc                                                                                              >  2:966313/1‑93 (MQ=255)
                                                                                                gCATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGGTTAg                                      >  2:541075/1‑143 (MQ=255)
                                                                                                gCATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAg                                      >  1:426503/1‑143 (MQ=255)
                                                                                                    catTATGTGTTGCCGGTCACGATTTGCTGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATCGGTCTGTGAGCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTgg                                  <  1:102940/143‑1 (MQ=255)
                                                                                                      ttATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTCGCTGGCa                                >  1:961486/1‑143 (MQ=255)
                                                                                                           ttttGCCGGTCACGATTTGATGAGGCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTGTTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATg                           >  2:424345/1‑143 (MQ=255)
                                                                                                            tttGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGGTGGCCGGATGc                          >  1:681022/1‑143 (MQ=255)
                                                                                                            tttGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGc                          >  2:636854/1‑143 (MQ=255)
                                                                                                               gCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAAGAAGCAcg                       >  2:931083/1‑143 (MQ=255)
                                                                                                                cctgTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCAcgc                      <  2:27221/140‑1 (MQ=255)
                                                                                                                ccGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGGCACGCCTGCTGCGCCGGTGATTAGCTCGCAGGATGCCcgc                      >  1:71142/1‑143 (MQ=255)
                                                                                                                ccGGTAACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCAcgc                      <  1:626491/143‑1 (MQ=255)
                                                                                                                caggTCCCGATTTGATGGGTCAGTGGGCCGGGAAGAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCAcgc                      <  1:994035/141‑1 (MQ=255)
                                                                                                                 cGGTCACGATTTGATGAGTCAGTTTGCTGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGAGCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGGCCCGCGTGCTGCGCCGTTGATTAGCTGCCGGCATGCAcgcg                     >  2:189438/1‑143 (MQ=255)
                                                                                                                   gTCACGATTTGATGAGTCAGTTTGCCGGGCATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTTCGTTATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGtt                   >  2:334163/1‑143 (MQ=255)
                                                                                                                          tttGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGcc            >  1:178657/1‑143 (MQ=255)
                                                                                                                            tgatgaGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTTCGTTATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCt          <  2:714278/143‑1 (MQ=255)
                                                                                                                              atgaGTCAGGCTGCCGGCAACAAAGAAGATATTGTCGGTATTGTTCTTGGATCCATCCGTTGCTGCCGTCCGTTATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTa        <  2:567040/143‑1 (MQ=255)
                                                                                                                                 attcggTTTGCCGGGAATACAGATGATATTGGCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGa     <  1:7924/138‑1 (MQ=255)
                                                                                                                                 aGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGa     >  1:174597/1‑143 (MQ=255)
                                                                                                                                   tCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATACGTTGCTGCCGTGCGTTATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAac   >  2:180270/1‑143 (MQ=255)
                                                                                                                                    cAGTGTGCCGGGAATAGCGAGGAGATTGTCGGTATTGGTCTGGGCGCCATCCGTTGCTTCCGTGCGTTATTGAAAGCTGATGGGACGCATGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAaca  <  2:271569/143‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                        
CTGGAAGGTAACGTGGTTTGCGAAGGCAAAGGCTTATTACAGCCGATGCACACGCCAGATGCCGATACTGCAGAACATCCTGACGACGATTTATGGGCATCATTATGTTTTGCCGGTCACGATTTGATGAGTCAGTTTGCCGGGAATAAAGAAGATATTGTCGGTATTGGTCTGGGATCCATCCGTTGCTGCCGTGCGTTATTGAAAGCCGATGGCACGCCTGCTGCGCCGTTGATTAGCTGGCAGGATGCACGCGTTACACGCCCTTACGAACA  >  NC_000913/2901962‑2902236

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: