Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,345,478 C→T 5.5% G276G (GGC→GGT acsA‑II → acetyl‑CoA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,345,4780CT 5.5% 213.4 / 2.2 73G276G (GGC→GGTacsA‑IIacetyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base C (37/31);  new base T (2/2);  total (39/34)
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.26e-01

GCGACCTCTGGTATCACGAAGCGACGAAGGATGCCGGCGACGACTGCCCCCCAGAACCGATGGAGGCCGAAGACCCGCTGTTCATCCTCTACACCTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTC  >  NC_002947/5345336‑5345618
                                                                                                                                              |                                                                                                                                            
gcgACCTCTGGTATCACGCAGCGACGAAGGATTCCGGCGACGACTGCCCCCCAGAACCGATGGAGGCCGAAGACCCGCTGTTCATCCTCTACACCTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACctgc                                                                                                                                       >  2:1688622/1‑150 (MQ=255)
     ctctGGTATCACGAAGCGACGAAGGATGCCGGCGACGTCTGCCCCCCAGAACCGATGGAGGCCGAAGACCCGCTGTTCATCCTCTACACCTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCttctgccg                                                                                                                                  >  2:187100/1‑148 (MQ=255)
     ctctGGTATCACGAAGCGACGAAGGATGCCGGCGACGACTGCCCCCCAGAACCGATGGAGGCCGAAGACCCGCTGTTCATCCTCTACACCTCCGGCAGCACCGCCAAACCCAAAGGGCTGCTGCCTACCACCGGCGGCTCCCTGCTGCac                                                                                                                                  >  2:674916/1‑149 (MQ=255)
               aCGAAGCGACGAAGGATGCCGGCGACGACTGCCCCCCAGAACCGATGGAGGCCGAAGACCCGCTGTTCATCCTCTACACCTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGa                                                                                                                        <  2:63588/150‑1 (MQ=255)
                  aaGCGACGAAGGATGCCGGCGACGACTGCCCCCCAGAACCGATGGAGGCCGAAGACCCGCTGTTCATCCTCTACACCTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCt                                                                                                                     >  2:1878527/1‑150 (MQ=255)
                        cgaAGGATGCCGGCGACGACTGCCCCCCAGAACCGATGGAGGCCGAAGACCCGCTGTTCATCCTCTACACCTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAgg                                                                                                               >  2:147916/1‑150 (MQ=255)
                            ggATGCCGGCGACGACTGCCCCCCAGAACCGATGGAGGCCGAAGACCCGCTGTTCATCCTCTACACCTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAggtggt                                                                                                           <  2:1524400/150‑1 (MQ=255)
                             gATGCCGGCGACGACTGCCCCCCAGAACCGATGGAGGCCGAAGACCCGCTGTTCATCCTCTACACCTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGgtggtg                                                                                                          >  1:177514/1‑150 (MQ=255)
                                gCCGGCGACGACTGCCCCCCAGAACCGATGGAGGCCGAAGACCCGCTGTTCATCCTCTACACCTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTc                                                                                                       <  2:805719/150‑1 (MQ=255)
                                 ccGGCGACGACTGCCCCCCAGAACCGATGGAGGCCGAAGACCCGCTGTTCATCCTCTACACCTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCg                                                                                                      >  2:1710904/1‑150 (MQ=255)
                                  cGGCGACGACTGCCCCCCAGAACCGATGGAGGCCGAAGACCCGCTGTTCATCCTCTACACCTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGa                                                                                                     >  1:248182/1‑150 (MQ=255)
                                      gacgacTGCCCCCCAGAACCGATGGAGGCCGAAGACCCGCTGTTCATCCTCTACACCTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTAt                                                                                                 >  1:1306650/1‑150 (MQ=255)
                                       acgacTGCCCCCCAGAACCGATGGAGGCCGAAGACCCGCTGTTCATCCTCTACACCTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATc                                                                                                <  2:1306650/150‑1 (MQ=255)
                                         gacTGCCCCCCAGAACCGATGGAGGCCGAAGACCCGCTGTTCATCCTCTACACCTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGt                                                                                              <  2:1097985/150‑1 (MQ=255)
                                          acTGCCCCCCAGAACCGATGGAGGCCGAAGACCCGCTGTTCATCCTCTACACCTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGt                                                                                              <  2:1878619/149‑1 (MQ=255)
                                          acTGCCCCCCAGAACCGATGGAGGCCGAAGACCCGCTGTTCATCCTCTACACCTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGt                                                                                              >  1:1878619/1‑149 (MQ=255)
                                                 cccAGAACCGATGGAGGCCGAAGACCCGCTGTTCATCCTCTACACCTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGTTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGa                                                                                      >  2:786581/1‑150 (MQ=255)
                                                 cccAGAACCGATGGAGGCCGAAGACCCGCTGTTCATCCTCTACACCTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGTTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGa                                                                                      >  2:786367/1‑150 (MQ=255)
                                                 cccAGAACCGATGGAGGCCGAAGACCCGCTGTTCATCCTCTACACCTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGa                                                                                      >  1:1579436/1‑150 (MQ=255)
                                                 cccAGAACCGATGGAGGCCGAAGACCCGCTGTTCATCCTATACACCTCCGGCAGCACCGGCAAACCCACAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGGAAAGGTGCACTTCAAGGCTGGGTTAGATTATCGTGGCGgaga                                                                                      >  1:396897/1‑150 (MQ=255)
                                                    aGAACCGCTGGCGGCCGAGGACCCGCTGTTCATCCTCTACACCTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTATCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGgt                                                                                   <  2:1003915/150‑1 (MQ=255)
                                                          gATGGAGGCCGAAGACCCGCTGTTCATCCTCTACACCTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTg                                                                             <  1:1839663/150‑1 (MQ=255)
                                                           aTGGAGGCCGAAGACCCGCTGTTCATCCTCTACACCTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGATACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTgg                                                                            <  1:869585/150‑1 (MQ=255)
                                                                 gCCGAAGACCCGCTGTTCATCCTCTACACCTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCACGGTGGTGTTCGATTATCGTGACGGCTAGTTGTTCTGGTGCAcc                                                                      >  2:1450968/1‑150 (MQ=255)
                                                                          ccGCTGTTCATCCTCTACACCTCCGGCAGCACCGGCAACCCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTc                                                             <  1:1450968/150‑1 (MQ=255)
                                                                          ccGCTGTTCATCCTCTACACCTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTc                                                             >  2:1830189/1‑150 (MQ=255)
                                                                          ccGCTGTTCATCCTCTACACCTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTc                                                             <  1:1100111/150‑1 (MQ=255)
                                                                           cGCTGTTCATCCTCTACACCTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCg                                                            <  1:1688622/150‑1 (MQ=255)
                                                                            gCTGTTCCTCCTCTACACCTCCGGCAGCACCGGCACTCCCAAGGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCCAGGTGGTGTTCGATTATCGTGACGGCGAGGTGCTCTGGTGCACCGCCGCTGTCgg                                                           <  2:396897/150‑1 (MQ=255)
                                                                               gTTCATCCTCTACACCTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGTTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTg                                                        <  1:786367/150‑1 (MQ=255)
                                                                               gTTCATCCTCTACACCTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGTTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTg                                                        <  1:786581/150‑1 (MQ=255)
                                                                                ttCATCCTCTACACCTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTgg                                                       <  1:147916/150‑1 (MQ=255)
                                                                                  cATCCTCTACACCTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGt                                                     <  2:1501796/150‑1 (MQ=255)
                                                                                    tCCTCTACACCTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGa                                                   >  2:1479293/1‑150 (MQ=255)
                                                                                    tCCNCTACACCTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGa                                                   >  2:1196318/1‑150 (MQ=255)
                                                                                      ctctACACCTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGAcc                                                 >  2:806596/1‑150 (MQ=255)
                                                                                      ctctACACCTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGAcc                                                 <  2:855956/150‑1 (MQ=255)
                                                                                        ctACACCTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAggtg                                                                                                             <  2:206366/88‑1 (MQ=255)
                                                                                        ctACACCTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAggtg                                                                                                             >  1:206366/1‑88 (MQ=255)
                                                                                            acCTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTgg                                                       >  2:1532112/1‑138 (MQ=255)
                                                                                            acCTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTgg                                                       <  1:1532112/138‑1 (MQ=255)
                                                                                            acCTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCcac                                           <  1:1830189/150‑1 (MQ=255)
                                                                                             ccTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCcaca                                          >  2:1175592/1‑150 (MQ=255)
                                                                                              cTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACTGCCGATGTCGGCTGGGTGACCGGCCACAg                                         >  2:442405/1‑150 (MQ=255)
                                                                                              cTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCgg                                                           <  2:218346/132‑1 (MQ=255)
                                                                                              cTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCgg                                                           >  1:218346/1‑132 (MQ=255)
                                                                                               tCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGa                                                                                      <  1:843997/104‑1 (MQ=255)
                                                                                               tCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGa                                                                                      >  2:843997/1‑104 (MQ=255)
                                                                                               tCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGc                                        <  1:1615176/150‑1 (MQ=255)
                                                                                                ccGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCt                                       <  1:1175592/150‑1 (MQ=255)
                                                                                                      gcaCCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCg                                 >  1:1333067/1‑150 (MQ=255)
                                                                                                        aCCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTc                               <  2:248182/150‑1 (MQ=255)
                                                                                                          cGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACctg                                                                                                                                        <  1:1115987/43‑1 (MQ=255)
                                                                                                          cGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACctg                                                                                                                                        >  2:1115987/1‑43 (MQ=255)
                                                                                                            gCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATg                           <  1:489793/150‑1 (MQ=255)
                                                                                                              aaaCCCAAAGGGGTGCTGCATACCGCCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGc                         >  2:651169/1‑150 (MQ=255)
                                                                                                                  ccAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGc                                                          >  2:1433581/1‑113 (MQ=255)
                                                                                                                  ccAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGc                                                          <  1:1433581/113‑1 (MQ=255)
                                                                                                                  ccAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGc                     >  1:568523/1‑150 (MQ=255)
                                                                                                                    aaaGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTg                   >  2:1831777/1‑150 (MQ=255)
                                                                                                                           tgctgcATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATg            >  1:176946/1‑150 (MQ=255)
                                                                                                                           tgctgcATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATg            <  1:1831777/150‑1 (MQ=255)
                                                                                                                              tgcATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATg            >  1:789019/1‑147 (MQ=255)
                                                                                                                              tgcATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATg            <  2:789019/147‑1 (MQ=255)
                                                                                                                              tgcATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGgcg         >  1:1853735/1‑150 (MQ=255)
                                                                                                                               gcATACCGCCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCAGCGGGCCatagacga        <  1:651169/150‑8 (MQ=255)
                                                                                                                                  tACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAAt     <  1:1878527/150‑1 (MQ=255)
                                                                                                                                   accaccGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAAtc    >  1:178261/1‑150 (MQ=255)
                                                                                                                                   accaccGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAAtc    >  1:1327923/1‑150 (MQ=255)
                                                                                                                                   accaccGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAAtc    >  1:1327974/1‑150 (MQ=255)
                                                                                                                                     caccGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAg                                                                                     >  2:98782/1‑67 (MQ=255)
                                                                                                                                     caccGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAg                                                                                     <  1:98782/67‑1 (MQ=255)
                                                                                                                                     caccGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAAtctc  >  2:1779617/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                            
GCGACCTCTGGTATCACGAAGCGACGAAGGATGCCGGCGACGACTGCCCCCCAGAACCGATGGAGGCCGAAGACCCGCTGTTCATCCTCTACACCTCCGGCAGCACCGGCAAACCCAAAGGGGTGCTGCATACCACCGGCGGCTACCTGCTGCAGGCAACGATGACCTTCAAGGTGGTGTTCGATTATCGTGACGGCGAGGTGTTCTGGTGCACCGCCGATGTCGGCTGGGTGACCGGCCACAGCTACATCGTCTATGGCCCGCTGGCCAATGGCGCAATCTC  >  NC_002947/5345336‑5345618

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