Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,273,095 A→G 7.6% intergenic (+17/+274) PP_2003 → / ← PP_2004 hypothetical protein/AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,273,0950AG 7.6% 217.4 / 7.1 79intergenic (+17/+274)PP_2003/PP_2004hypothetical protein/AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (31/42);  new base G (3/3);  total (34/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.19e-01

AGCGAGGCGCCAGCGGCGCCGGTGGTGGATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGG  >  NC_002947/2272947‑2273243
                                                                                                                                                    |                                                                                                                                                    
aGCGAGGCGCCAGCGGCGCCGGTGGTGGATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGAc                                                                                                                                                     <  2:1780760/150‑1 (MQ=255)
aGCGAGGCGCCAGCGGCGCCGGTGGTGGATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGAc                                                                                                                                                     <  1:370144/150‑1 (MQ=255)
 gCGAGGCGCCAGCGGCGCCGGTGGTGGATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGAcc                                                                                                                                                    >  1:1778248/1‑150 (MQ=255)
     ggCGCCAGCGGCGCCGGTGGTGGATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCttt                                                                                                                                                <  1:1120155/150‑1 (MQ=255)
      gcgcCAGCGGCGCCGGTGGTGGATGACAGCCAGGGTTCGGAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTc                                                                                                                                               <  1:460021/150‑1 (MQ=255)
      gcgcCAGCGGCGCCGGTGGTGGATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTc                                                                                                                                               <  1:1904848/150‑1 (MQ=255)
         ccAGCGGCGCCGGTGGTGGATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCTGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGcc                                                                                                                                            <  2:1719748/150‑1 (MQ=255)
           aGCGGCGCCGGTGGTGGATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTTAACCGACCCTTTCGCCgg                                                                                                                                          <  2:1142620/150‑1 (MQ=255)
                cgcCGGTGGTGGATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGcc                                                                                                                                     >  2:951207/1‑150 (MQ=255)
                    ggtggtGGATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCcgcg                                                                                                                                 >  1:376086/1‑150 (MQ=255)
                        gtggATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAg                                                                                                                             <  2:1443173/150‑1 (MQ=255)
                         tggATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCGGGCTGCCGGCGAAAgg                                                                                                                            <  1:278079/150‑1 (MQ=255)
                         tggATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAgg                                                                                                                            <  2:1623164/150‑1 (MQ=255)
                          ggATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAggg                                                                                                                           <  1:14059/150‑1 (MQ=255)
                           gATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTTATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGc                                                                                                                          <  1:1194310/150‑1 (MQ=255)
                             tGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCg                                                                                                                        <  1:951207/150‑1 (MQ=255)
                             tGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCg                                                                                                                        <  1:903671/150‑1 (MQ=255)
                             tGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCg                                                                                                                        <  2:564364/150‑1 (MQ=255)
                               aCAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGt                                                                                                                      >  2:808355/1‑150 (MQ=255)
                                 aGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTc                                                                                                                    <  1:721728/150‑1 (MQ=255)
                                        gTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTc                                                                                                                                               <  2:1187853/116‑1 (MQ=255)
                                        gTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTc                                                                                                                                               >  1:1187853/1‑116 (MQ=255)
                                             cAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGGCCCTTTCGCGGGCTGCCGGCGAAAgggtc                                                                                                                         >  2:1435590/1‑131 (MQ=38)
                                             cAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGGCCCTTTCGCGGGCTGCCGGCGAAAgggtc                                                                                                                         <  1:1435590/133‑3 (MQ=38)
                                                 cTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATccc                                                                                                    <  1:934025/150‑1 (MQ=255)
                                                    tgGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTg                                                                                                 >  2:1934883/1‑150 (MQ=255)
                                                      gCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACAt                                                                                                       <  2:1289915/142‑1 (MQ=255)
                                                      gCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACAt                                                                                                       >  1:1289915/1‑142 (MQ=255)
                                                       ccccGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGGCCCTTTCGCGGGCTGCCGGCGaaa                                                                                                                              <  2:1868278/118‑1 (MQ=38)
                                                       ccccGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGGCCCTTTCGCGGGCTGCCGGCGaaa                                                                                                                              >  1:1868278/1‑118 (MQ=38)
                                                       ccccGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTa                                                                                              <  1:134895/150‑1 (MQ=255)
                                                        cccGCCGGCCGAAGCTGCACCCAGTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTAc                                                                                             >  2:141650/1‑150 (MQ=255)
                                                         ccgccgGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACa                                                                                            <  1:1680952/150‑1 (MQ=255)
                                                              ggccgAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGa                                                                                                                                                      >  2:383865/1‑87 (MQ=255)
                                                              ggccgAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGa                                                                                                                                                      <  1:383865/87‑1 (MQ=255)
                                                               gccgAAGCTGCACCCATTCAGGTACAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCcgcg                                                                                                                                 <  1:37964/107‑1 (MQ=255)
                                                               gccgAAGCTGCACCCATTCAGGTACAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCcgcg                                                                                                                                 >  2:37964/1‑107 (MQ=255)
                                                                  gAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCcg                                                                                   >  1:1471574/1‑150 (MQ=255)
                                                                     gCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGa                                                                                <  2:316539/150‑1 (MQ=255)
                                                                          aCCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATccc                                                                                                    <  1:385653/125‑1 (MQ=255)
                                                                          aCCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATccc                                                                                                    >  2:385653/1‑125 (MQ=255)
                                                                          aCCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACt                                                                           >  1:201207/1‑150 (MQ=255)
                                                                              aTTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTCCt                                                                       <  1:1934883/150‑1 (MQ=255)
                                                                               ttCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTg                                                                      <  2:1778248/150‑1 (MQ=255)
                                                                                   ggTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAaca                                                                  >  1:1369216/1‑150 (MQ=255)
                                                                                        agcagGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGc                                                             >  2:377624/1‑150 (MQ=255)
                                                                                                  ccAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCCACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGTCCCGCGATAACTTTCTGAACACCTGCTACGTTTAcc                                                   >  1:415719/1‑150 (MQ=255)
                                                                                                     gccggccgAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAgg                                                                                                                 <  1:603436/85‑1 (MQ=255)
                                                                                                     gccggccgAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAgg                                                                                                                 >  2:603436/1‑85 (MQ=255)
                                                                                                        ggccgAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCAt                                             <  2:201207/150‑1 (MQ=255)
                                                                                                         gccgAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATa                                            >  1:808333/1‑150 (MQ=255)
                                                                                                          ccgAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGATCTTACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATaa                                           <  1:1373327/150‑1 (MQ=255)
                                                                                                          ccgAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTAcc                                                   <  2:203911/142‑1 (MQ=255)
                                                                                                          ccgAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTAcc                                                   >  1:203911/1‑142 (MQ=255)
                                                                                                           cgAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAg                                          >  1:1658845/1‑150 (MQ=255)
                                                                                                           cgAGGCCGCCGCCCCCTCCCTTTGATGGGCCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAg                                          >  2:523266/1‑150 (MQ=255)
                                                                                                             aGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCt                                        <  1:377624/150‑1 (MQ=255)
                                                                                                             aGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCt                                        >  1:759923/1‑150 (MQ=255)
                                                                                                              ggCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCtt                                       <  2:1471574/150‑1 (MQ=255)
                                                                                                              ggCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCtt                                       >  1:117583/1‑150 (MQ=255)
                                                                                                               gccgccgccCCCTCCCTTTGATGGGTCGCCTGAACCGGCCCTTTCGCGGGCTGCCGGCGAAAGGGTCGGTTCAg                                                                                                                  <  1:863406/74‑1 (MQ=14)
                                                                                                               gccgccgccCCCTCCCTTTGATGGGTCGCCTGAACCGGCCCTTTCGCGGGCTGCCGGCGAAAGGGTCGGTTCAg                                                                                                                  >  2:863406/1‑74 (MQ=14)
                                                                                                                 cgccgccCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGtc                                    <  2:639808/150‑1 (MQ=255)
                                                                                                                   ccgccCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGtcat                                  >  1:1131329/1‑150 (MQ=255)
                                                                                                                        ccctccctTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCt                                                            >  2:1502106/1‑119 (MQ=255)
                                                                                                                        ccctccctTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCt                                                            <  1:1502106/119‑1 (MQ=255)
                                                                                                                          ctccctTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCt                           <  1:46315/150‑1 (MQ=255)
                                                                                                                           tccctTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCtt                          <  2:117583/150‑1 (MQ=255)
                                                                                                                              ctTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCAt                       >  1:142390/1‑150 (MQ=255)
                                                                                                                              ctTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTCAACACCCGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCAt                       <  2:376086/150‑1 (MQ=255)
                                                                                                                               tttGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATc                      >  1:1141260/1‑150 (MQ=255)
                                                                                                                                  gATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATcccc                   >  2:1343657/1‑149 (MQ=255)
                                                                                                                                        tCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCt             >  1:1275037/1‑150 (MQ=255)
                                                                                                                                        tCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCt             >  1:740589/1‑150 (MQ=255)
                                                                                                                                          gCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTa           >  2:1918218/1‑150 (MQ=255)
                                                                                                                                           ccTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAg          <  1:808355/150‑1 (MQ=255)
                                                                                                                                            cTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAgg         <  1:523266/150‑1 (MQ=255)
                                                                                                                                            cTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAgg         <  2:740589/150‑1 (MQ=255)
                                                                                                                                            cTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCAGCGCAACGGTGCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAgg         >  2:1611231/1‑150 (MQ=255)
                                                                                                                                                 ccGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGaa    <  2:1658845/150‑1 (MQ=255)
                                                                                                                                                   gACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAgg  <  2:1141260/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGCGAGGCGCCAGCGGCGCCGGTGGTGGATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGG  >  NC_002947/2272947‑2273243

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: