Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,004,238 T→A 8.8% V259D (GTC→GAC)  PP_2627 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,004,2380TA 8.8% 217.2 / 10.2 70V259D (GTC→GAC) PP_2627membrane protein
Reads supporting (aligned to +/- strand):  ref base T (28/35);  new base A (3/3);  total (31/39)
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.16e-01

TCCGCGCCGCCCGGCATCCAGGCCCGGTTATGGGAGCACCTGTTGGGCTGGCTGCAGCGTAAGCGAAGCCAGCGGGCGCTCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTGGGGGCGCGCCTGCCGCTGTACGTGGTACTGACCAAATTCGATCTGCTGGACGGCTTCGACCAGTTCTACAGCCAGCTTCCGGCCGCCAAGCGTAAA  >  NC_002947/3004096‑3004386
                                                                                                                                              |                                                                                                                                                    
tCCGCGCCGCCCGGCATCCAGGCCCGGTTATGGGAGCACCTGTTGGGCTGGCTGCAGCGTAAGCGAAGCCAGCGGGCGCTCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTg                                                                                                                                               >  1:1648038/1‑150 (MQ=255)
tCCGCGCCGCCCGGCATCCAGGCCCGGTTATGGGAGCACCTGTTGGGCTGGCTGCAGCGTAAGCGAAGCCAGCGGGCGCTCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTg                                                                                                                                               <  1:1828817/150‑1 (MQ=255)
tCCGCGCCGCCCGGCATCCAGGCCCGGTTATGGGAGCACCTGTTGGGCTGGCTGCAGCGTAAGCGAAGCCAGCGGGCGCTCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGACCGAGCAGCGCGTCGCGCTg                                                                                                                                               >  1:211266/1‑150 (MQ=255)
  cgcgCCGCCCGGCATCCAGGCCCGGTTATGGGAGCACCTGTTGGGCTGGCTGCAGCGTAAGCGAAGCCAGCGGGCGCTCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGc                                                                                                                                             <  1:1411656/150‑1 (MQ=255)
    cgccgcCCGGCATCCAGGCCCGGTTATGGGAGCACCTGTTGGGCTGGCTGCAGCGTAAGCGAAGCCAGCGGGCGCTCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGACCGAGCAGCGCGACGCGCTgctcg                                                                                                                                           <  1:75035/150‑5 (MQ=255)
       cgccCGGCATCCAGGCCCGGTTATGGGAGCACCTGTTGGGCTGGCTGCAGCGTAAGCGAAGCCAGCGGGCGCTCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATg                                                                                                                                        >  2:1165989/1‑150 (MQ=255)
        gccCGGCATCCAGGCCCGGTTATGGGAGCACCTGTTGGGCTGGCTGCAGCGTAAGCGAAGCCAGCGGGCGCTCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATgg                                                                                                                                       <  1:1165989/150‑1 (MQ=255)
              cagccAGGCCCGGTTATGGGAGCACCTGTTGGGCTGGCTGCAGCGTAAGCGAAGCCAGCGGGCGCTCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTgcg                                                                                                                                 <  1:238288/147‑1 (MQ=255)
                       ccGGTTATGGGAGCAACTGTTGGGCTGGCTGCAGCGTAAGCGAAGCCAGCGGGCGCTCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCt                                                                                                                        >  2:348473/1‑150 (MQ=255)
                        cGGTTATGGGAGCACCTGTTGGGCTGGCTGCAGCGTAAGCGAAGCCAGCGGGCGCTCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTgcgc                                                                                                                                >  1:1137381/1‑141 (MQ=255)
                        cGGTTATGGGAGCACCTGTTGGGCTGGCTGCAGCGTAAGCGAAGCCAGCGGGCGCTCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTgcgc                                                                                                                                <  2:1137381/141‑1 (MQ=255)
                         ggTTATGGGAGCACCTGTTGGGCTGGCTGCAGCGTAAGCGAAGCCAGCGGGCGCTCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGt                                                                                                                      >  1:1334000/1‑150 (MQ=255)
                              tGGGAGCACCTGTTGGGCTGGCTGCAGCGTAAGCGAAGCCAGCGGGCGCTCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGaa                                                                                                                 <  2:2183204/150‑1 (MQ=255)
                                ggAGCACCTGTTGGGCTGGCTGCAGCGTAAGCGAAGCCAGCGGGCGCTCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGt                                                                                                               <  2:308459/150‑1 (MQ=255)
                                ggAGCACCTGTTGGGCTGGCTGCAGCGTAAGCGAAGCCAGCGGGCGCTCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGt                                                                                                               <  1:590092/150‑1 (MQ=255)
                                 gAGCACCTGTTGGGCTGGCTGCAGCGTAAGCGAAGCCAGCGGGCGCTCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTg                                                                                                              >  1:167535/1‑150 (MQ=255)
                                  aGCACCTGTTGGGCTGGCTGCAGCGTAAGCGAAGCCAGCGGGCGCTCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGa                                                                                                             <  2:1536027/150‑1 (MQ=255)
                                  aGCACCTGTTGGGCTGGCTGCAGCGTAAGCGAAGCCAGCGGGCGCTCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGa                                                                                                             <  2:1334000/150‑1 (MQ=255)
                                  aGCACCTGTTGGGCTGGCTGCAGCGTAAGCGAAGCCAGCGGGCGCTCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGa                                                                                                             <  2:326098/150‑1 (MQ=255)
                                     aCCTGTTGGGCTGGCTGCAGCGTAAGCGAAGCCAGCGGGCGCTCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGca                                                                                                          <  1:1148424/150‑1 (MQ=255)
                                       cTGTTGGGCTGGCTGCAGCGTAAGCGAAGCCAGCGGGCGCTCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGcagt                                                                                                        <  2:719198/150‑1 (MQ=255)
                                         gTTGGGCTGGCTGCAGCGTAAGCGAAGCCAGCGGGCGCTCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGCGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGGGAGcagtca                                                                                                      <  2:637873/150‑1 (MQ=255)
                                         gTTGGGCTGGCTGCAGCGTAAGCGAAGCCAGCGGGCGCTCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGcagtca                                                                                                      <  1:320857/150‑1 (MQ=255)
                                         gTTGGGCTGGCTGCAGCGTAAGCGAAGCCAGCGGGCGCTCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGcagtca                                                                                                      <  1:1199119/150‑1 (MQ=255)
                                          ttGGGCTGGCTGCAGCGTAAGCGAAGCCAGCGGGCGCTCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGTACACACGCCTGTACGAAGTGAGcagtcag                                                                                                     <  1:1642997/150‑1 (MQ=255)
                                                tggctgCAGCGTAAGCGAAGCCAGCGGGCGCTCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTgggg                                                                                               >  2:1854921/1‑150 (MQ=255)
                                                    tgCAGCGTAAGCGAAGCCAGCGGGCGCTCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTGGGGgcgc                                                                                           <  2:1616735/150‑1 (MQ=255)
                                                          gTAAGCGAAGCCAGCGGGCGCTCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTGGGGGCGCGCCTGc                                                                                     <  2:167535/150‑1 (MQ=255)
                                                                  aGCCAGCGGGCGCTCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTGGGGGCGCGCCTGCCGCTGTAc                                                                             >  2:1760873/1‑150 (MQ=255)
                                                                    ccAGCGGGCGCTCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTGGGGGCGCGCCTGCCGCTGTACGt                                                                           >  2:105871/1‑150 (MQ=255)
                                                                     cAGCGGGCGCTCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCAGCATGGCCGGGCCGAGCAGCGCGCCGCGCTGGCCTATGGGTTGCGCACACGCCTGTACGAAGTGAGAAGTCAGTTGGGGGCGCGCCTGCCGCTGTACGTg                                                                          <  2:1478164/150‑1 (MQ=255)
                                                                      aGCGGGCGCTCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTGGGGGCGCGCCTGCCGCTGTACGTgg                                                                         <  1:2115051/150‑1 (MQ=255)
                                                                           gcgcTCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTGGGGGCGCGCCTGCCGCTGTACGTgg                                                                         <  2:875973/145‑1 (MQ=255)
                                                                           gcgcTCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTGGGGGCGCGCCTGCCGCTGTACGTgg                                                                         >  1:875973/1‑145 (MQ=255)
                                                                           gcgcTCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTGGGGGCGCGCCTGCCGCTGTACGTGGTACTg                                                                    >  1:955744/1‑150 (MQ=255)
                                                                           gcgcTCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTGGGGGCGCGCCCGCCGCTGTACGTGGTACTg                                                                    >  2:1764042/1‑150 (MQ=255)
                                                                              cTCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTGGGGGCGCGCCTGCCGCTGTACGTGGTACTGAcc                                                                 >  2:1826940/1‑150 (MQ=255)
                                                                               tCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTGGGGGCGCGCCTGCCGCt                                                                                 <  1:1804763/133‑1 (MQ=255)
                                                                               tCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTGGGGGCGCGCCTGCCGCt                                                                                 >  2:1804763/1‑133 (MQ=255)
                                                                                 aaTGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGACGCGCTgctc                                                                                                                                            <  1:1836183/72‑4 (MQ=255)
                                                                                 aaTGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGACGCGCTgctc                                                                                                                                            >  2:1836183/1‑69 (MQ=255)
                                                                                      aGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTGGGGGCGCGCCTGCCGCTGTACGTGGTACTGACCAAATTCGa                                                         <  1:1780168/150‑1 (MQ=255)
                                                                                       gtgtTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTGGGGGCGCGCCTGCCGCTGTACGTGGTACTGACCAAATTCGAt                                                        <  1:2353727/150‑1 (MQ=255)
                                                                                       gtgtTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACAAGCCTGTACGAAGTGAGCAGTCAGTTGGGGGCGCGCCTGCCGCTGTACGTGGTACTGAACAAATTCGAt                                                        <  1:1174274/150‑1 (MQ=255)
                                                                                          ttGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTGGGGGCGCGCCTGCCGCTGTACGTGGTACTGACCAAATTCGATctg                                                     >  1:1789660/1‑150 (MQ=255)
                                                                                              tggtggTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTGGGGGCGCGCCTGCCGCTGTACGTGGTACTGACCAAATTCGATCTGCTgg                                                 >  2:2007708/1‑150 (MQ=255)
                                                                                                 tggtCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTGGGGGCGCGCCTGCCGCTGTACGTGGTACTGACCAAATTCGATCTGCTGGACg                                              <  2:2351320/150‑1 (MQ=255)
                                                                                                  ggtCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTGGGGGCGCGCCTGCCGCTGTACGTGGTACTGACCAAATTCGATCTGCTGGACgg                                             <  1:348473/150‑1 (MQ=255)
                                                                                                   gtCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTGGGGGCGCGCCTGCCGCTGTACGTGGTACTGACCAAATTCGATCTGCTGGACGGc                                            >  2:2356464/1‑150 (MQ=255)
                                                                                                      gACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCTTGGGTTGCGCACCTGCCTGTACGAAGTGAGCAGTCAGTTGGGGGCGCGCCTGCCGCTGTACGTGGTACTGGTCAAATTCGATCTGATGGACGGCttt                                         >  1:212808/1‑149 (MQ=255)
                                                                                                      gACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTGGGGGCGCGCCTGCCGCTGTACGTGGTACTGACCAAATTCGATCTGCTGGACGGCTTc                                         >  2:2365874/1‑150 (MQ=255)
                                                                                                        ccTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTGGGGGCGCGCCTGCCGCTGTACGTGGTACTGACCAAATTCGATCTGCTGGACGGCTTCGa                                       <  2:1641260/150‑1 (MQ=255)
                                                                                                         cTGCCGGCGCTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTGGGGGCGCGCCTGCCGCTGTACGTGGTACTGACCAAATTCGATCTGCTGGACGGCtt                                          >  1:2158670/1‑146 (MQ=255)
                                                                                                         cTGCCGGCGCTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTGGGGGCGCGCCTGCCGCTGTACGTGGTACTGACCAAATTCGATCTGCTGGACGGCtt                                          <  2:2158670/146‑1 (MQ=255)
                                                                                                          tGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTGGGGGCGCGCCTGCCGCTGTACGTGGTACTGACCAAATTCGATCTGCTGGACGGCTTCGAcc                                     >  1:1414940/1‑150 (MQ=255)
                                                                                                              ggCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTGGGGGCGCGCCTGCCGCTGTACGTGGTACTGACCAAATTCGATCTGCTGGACGGCTTCGACCAGtt                                 >  2:2288674/1‑150 (MQ=255)
                                                                                                                 gTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTGGGGGCGCGCCTGCCGCTGTACGTGGTACTGACCAAATTCGATCTGCTGGACGGCTTCGACCAGTTCTa                              >  1:603964/1‑150 (MQ=255)
                                                                                                                  ttGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTGGGGGCGCGCCTGCCGCTGTACGTGGTACTGACCAAATTCGATCTGCTGGACGGCTTCGACCAGTTCTAc                             <  2:955744/150‑1 (MQ=255)
                                                                                                                      tgcATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTGGGGgcgcgc                                                                                         >  2:1127644/1‑86 (MQ=255)
                                                                                                                      tgcATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTGGGGgcgcgc                                                                                         <  1:1127644/86‑1 (MQ=255)
                                                                                                                       gcATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTGGGGGCGCGCCTGCCGCTGTACGTGGTACTGACCAAATTCGATCTGCTGGACGGCTTCGACCAGTTCTAcagcca                        <  2:211266/150‑1 (MQ=255)
                                                                                                                       gcATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTGGGGGCGCGCCTGCCGCTGTACGTGGTACTGACCAAATTCGATCTGCTGGACGGCTTCGACCAGTTCTAcagcca                        <  2:1789660/150‑1 (MQ=255)
                                                                                                                            gCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTGGGGGCGCGCCTGCCGCTGTACGTGGTACTGACCAAATTCGATCTGCTGGACGGCTTCGACCAGTTCTACAGCCAGCTTc                   <  1:1273144/150‑1 (MQ=255)
                                                                                                                              aGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTGGGGGCGCGCCTGCCGCTGTACGTgg                                                                         <  1:455716/94‑1 (MQ=255)
                                                                                                                              aGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTGGGGGCGCGCCTGCCGCTGTACGTgg                                                                         >  2:455716/1‑94 (MQ=255)
                                                                                                                                gCCCGAGAAGCGCGTCGCGCTGGACCATGGGTTGGGCACACGCCTGTAGGAAGTGAGCAGTCAGTTGGGGGCGCGCCTGCCGCTGTACGTGGTAATGACCAAATTCGATCTGCTGGACGACTTCGAGTAGGTCTACAGCCAGCTTCAggg               >  2:982081/1‑149 (MQ=255)
                                                                                                                                     gccagcGCGACGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTGGGGGCGCGCCTGCAGCGATCTTTATTTCTTACCAAATTCGATCTGCTGGCCGGCTTCGACCAGTTCTACAGCCAGCTTCCGGCCGCCa          >  1:2223438/3‑150 (MQ=255)
                                                                                                                                      ccagcGCGACGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTGGGGGCGCGCCTGCCGCTGTACGTGGTACTGACCAAATTCGATCTGCTGGACGGCTTCGACCAGTTCTAcagcca                        >  2:591661/2‑135 (MQ=255)
                                                                                                                                      ccagcGCGACGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTGGGGGCGCGCCTGCCGCTGTACGTGGTACTGACCAAATTCGATCTGCTGGACGGCTTCGACCAGTTCTAcagcca                        <  1:591661/134‑1 (MQ=255)
                                                                                                                                             gTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTGGGGGCGCGCCTGCCGCTGTACGTGGTACTGACCAAATTCGATCTGCTGGACGGCTTCGACCAGTTCTACAGCCAGCTTCCGGCCGCCAAGCGTaaa  >  1:570153/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
TCCGCGCCGCCCGGCATCCAGGCCCGGTTATGGGAGCACCTGTTGGGCTGGCTGCAGCGTAAGCGAAGCCAGCGGGCGCTCAATGGAGTGTTGCTGGTGGTCGACCTGCCGGCGTTGCTGCATGGCAGGCCCGAGCAGCGCGTCGCGCTGGCCCATGGGTTGCGCACACGCCTGTACGAAGTGAGCAGTCAGTTGGGGGCGCGCCTGCCGCTGTACGTGGTACTGACCAAATTCGATCTGCTGGACGGCTTCGACCAGTTCTACAGCCAGCTTCCGGCCGCCAAGCGTAAA  >  NC_002947/3004096‑3004386

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