Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,496,582 T→C 30.7% intergenic (+17/‑339) lpxP → / → ypdK palmitoleoyl‑acyl carrier protein (ACP)‑dependent acyltransferase/inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,496,5820TC30.7% 99.2 / 66.1 85intergenic (+17/‑339)lpxP/ypdKpalmitoleoyl‑acyl carrier protein (ACP)‑dependent acyltransferase/inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (25/33);  new base C (15/11);  total (40/45)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATGCGCGCACCGGAGCAGTACCTCTGGATCCACCGTCGCTTTAAAACGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCTTTTTTCTTAAGGTTTTTTCCTTCGCTGAATATCTAATTATTACCGTCAGAAAT  >  NC_000913/2496484‑2496680
                                                                                                  |                                                                                                  
cATGCGCGCACCGGAGCAGTACCTCTGGATCCACCGTCGCTTTAAAACGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCtt                                                                                                   >  2:1157954/1‑100 (MQ=255)
  tGCGCGCACCGGAGCAGTACCTCTGGATCCACCGTCGCTTTAAAACGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCTTTc                                                                                                 <  2:2082871/100‑1 (MQ=255)
  tGCGCGCACCGGAGCAGTACCTCTGGATCCACCGTCGCTTTAAAACGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCTTTc                                                                                                 <  2:1511718/100‑1 (MQ=255)
   gcgcgcACCGGAGCAGTACCTCTGGATCCACCGTCGCTTTAAAACGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCTTTcc                                                                                                <  2:696220/100‑1 (MQ=255)
    cgcgcACCGGAGCAGTACCTCTGGATCCACCGTCGCTTTAAAACGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCTTTccc                                                                                               <  1:1258015/100‑1 (MQ=255)
    cgcgcACCGGAGCAGTACCTCTGGATCCACCGTCGCTTTAAAACGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCCTTcc                                                                                                >  1:1892049/1‑99 (MQ=255)
     gcgcACCGGAGCAGTACCTCTGGATCCACCGTCGCTTTAAAACGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCTTTccc                                                                                               <  2:952043/99‑1 (MQ=255)
      ngcACCGGAGCAGTACCTCTGGATCCACCGTCGCTTTAAAACGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCTTTCcctc                                                                                             >  2:2114290/2‑100 (MQ=255)
       gcACCGGAGCAGTACCTCTGGATCCACCGTCGCTTTAAAACGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCTTTCcctcc                                                                                            <  1:1856369/100‑1 (MQ=255)
         aCCGGAGCAGTACCTCTGGATCCACTGTCGCTTTAAATCGCGCCCGGTGGGAGTATCGTCGTTGTACATTTAACAAACTGCTCTTACTCTTTccc                                                                                               <  2:314600/95‑1 (MQ=255)
         aCCGGAGCAGTACCTCTGGATCCACCGTCGCTTTAAAACGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCTTTCCctcctc                                                                                          <  1:186166/100‑1 (MQ=255)
         aCCGGAGCAGTACCTCTGGATCCACCGTCGCTTTAAAACGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCCTTCcctcc                                                                                            <  1:209628/98‑1 (MQ=255)
           cGGAGCAGTACCTCTGGATCCACCGTCGCTTTAAAACGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCACTTACTCTTTCCCTCCTCaa                                                                                        >  1:1124266/1‑100 (MQ=255)
                 aGTACCTCTGGATCCACCGTCGCTTTAAAACGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCGGc                                                                                  >  1:1348397/1‑100 (MQ=255)
                  gTACCTCTGGATCCACCGTCGCTTTAAAACGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCGGcc                                                                                 <  2:592739/100‑1 (MQ=255)
                  gTACCTCTGGATCCACCGTCGCTTTAAAACGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCCTTCCCTCCTCAAATCGGcc                                                                                 <  2:381005/100‑1 (MQ=255)
                    aCCTCTGGATCCACCGTCGCTTTAAAACGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCCTTCCCTCCTCAAATCGGcc                                                                                 <  2:1359290/98‑1 (MQ=255)
                     ccTCTGGATCCACCGTCGCTTTAAAACGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCGGCCata                                                                              <  2:1887893/100‑1 (MQ=255)
                      ctctGGATCCACCGTCGCTTTAAAACGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCGGCCataa                                                                             >  1:447453/1‑100 (MQ=255)
                      ctctGGATCCACCGTCGCTTTAAAACGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCGGCCataa                                                                             >  2:1306360/1‑100 (MQ=255)
                      ctctGGATCCACCGTCGCTTTAAAACGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCGGCCataa                                                                             >  1:1372453/1‑100 (MQ=255)
                      ctctGGATCCACCGTCGCTTTAAAACGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCGGCCataa                                                                             >  1:2136589/1‑100 (MQ=255)
                       tctGGATCCACCGTCGCTTTAAAACGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCCAATCGGCCataat                                                                            >  1:50743/1‑100 (MQ=255)
                         tGGATCCACCCTCGCTTTAATACGCGCCCGGTGGGTGAATCGTCGTTGTACATTTAACAAACTGCTCTTTCACTTTCCCTCCTCAATTCGGCCataattg                                                                          <  2:1283253/100‑3 (MQ=255)
                          ggATCCACCGTCGCTTTAAAACGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCGGCCATAATAAc                                                                         >  2:1062014/1‑100 (MQ=255)
                               cACCGTCGCTTTAAAACGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCCTTCCCTCCTCAAATCGGCCATAATAACCTCAt                                                                    >  1:1497187/1‑100 (MQ=255)
                                 ccGTCGCTTTAAAACGCGCCAGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCCTTCCCTCCTCAAATCGGCCATAATAACCTCATgg                                                                  <  1:2236820/100‑1 (MQ=255)
                                 ccGTCGCTCTAAAACGCACCCGGTGAGAGAACCGTCGCTGCACATTTATCAACCTACTCTTACTCTTTCCCTCCTCAATTCGGCCATAAAAAccatga                                                                    <  2:249829/98‑5 (MQ=255)
                                     cGCTTTAAAACGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCGGCCATAATAACCTCATGGtttt                                                              <  1:1472199/100‑1 (MQ=255)
                                     cGCTTTAAAACGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCGGCCATAATAACCTCATGGtttt                                                              >  2:1537530/1‑100 (MQ=255)
                                       cTTTAAAACGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCGGCCATCATAACCTCTTGGTTTTTa                                                            >  2:75262/1‑100 (MQ=255)
                                       cTTTAAAACGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTa                                                            >  1:226488/1‑100 (MQ=255)
                                       cTTTAAAACGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTa                                                            >  1:633525/1‑100 (MQ=255)
                                       cTTTAAAACGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCCTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTa                                                            >  2:1251056/1‑100 (MQ=255)
                                       cTTTAAAACGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCCTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTa                                                            >  2:1214826/1‑100 (MQ=255)
                                        tttAAAACGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTAt                                                           <  1:1157958/100‑1 (MQ=255)
                                        tttAAAACGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTAt                                                           <  2:447454/100‑1 (MQ=255)
                                         ttAAAACGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATg                                                          >  1:789974/1‑100 (MQ=255)
                                           aaaaCGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGc                                                        >  1:148008/1‑100 (MQ=255)
                                            aaaCGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAATTGCTCTTACTCCTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCt                                                       >  1:892513/1‑100 (MQ=255)
                                            aaaCGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCt                                                       >  1:1297935/1‑100 (MQ=255)
                                            aaaCGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCt                                                       >  1:1563032/1‑100 (MQ=255)
                                            aaaCGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCt                                                       >  2:567006/1‑100 (MQ=255)
                                            aaaCGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCCTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATg                                                          >  1:180741/1‑97 (MQ=255)
                                            aaaCGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCCTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCt                                                       >  1:1990014/1‑100 (MQ=255)
                                               cgcgCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCtttt                                                    >  2:2060438/1‑100 (MQ=255)
                                               cgcgCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCCTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCtttt                                                    >  1:1504274/1‑100 (MQ=255)
                                                  gcCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCGGCCATAATATCCTCATGGTTTTTATGGCTTTTTTc                                                 <  2:342013/100‑1 (MQ=255)
                                                  gcCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCTTTTTTc                                                 <  1:2144371/100‑1 (MQ=255)
                                                   cccGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCCTTCCCTCCTCAAATCGACCATAATAACCTCAAGGTTTTTCTCGCtttttt                                                  >  2:264382/1‑98 (MQ=255)
                                                     cGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCTTTTTTCTTa                                              <  1:1395915/100‑1 (MQ=255)
                                                     cGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCCTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCTTTTTTCTTa                                              >  1:1349781/1‑100 (MQ=255)
                                                     cGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCCTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCTTTTTTCTTa                                              <  1:737467/100‑1 (MQ=255)
                                                      ggTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCTTTTTTCTTaa                                             <  2:570020/100‑1 (MQ=255)
                                                      ggTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCTTTTTTCTTaa                                             <  1:1406776/100‑1 (MQ=255)
                                                        tGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCTTTTTTCTTaa                                             <  1:1567996/98‑1 (MQ=255)
                                                        tGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCTTTTTTCTTAAgg                                           <  2:1447337/100‑1 (MQ=255)
                                                            agaATCGTCGTTGTACATTTAACAAACTGCTCTTACTCCTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCTTTTTTCTTAAGGtttt                                       <  1:861741/100‑1 (MQ=255)
                                                             gaATCGTCGTTGTACATTTAACAAACTGCTCTTACTCCTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCTTTTTTCTTAAGGttttt                                      <  2:92528/100‑1 (MQ=255)
                                                              aaTCGTCGTTGTACATTTAACAAACTGCTCTTACTCGTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCTTTTTTCTTAAGGtttttt                                     <  1:456761/100‑1 (MQ=255)
                                                              aaTCGTCGTTGTACATTTAACAAACTGCTCTTACTCCTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCTTTTTTCTTAAGGtttttt                                     <  2:648949/100‑1 (MQ=255)
                                                               aTCGTCGTTGTACATTTAACAAACTGCTCTTACTCCTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCTTTTTTCTTAAGGTTTTTTc                                    >  1:613627/1‑100 (MQ=255)
                                                                tcgtcgTTGTACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCTTTTTTCTTAAGGTTTTTTcc                                   <  1:1151116/100‑1 (MQ=255)
                                                                 cgtcgtTGTACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCTTTTTTCTTAAGGttttttt                                    >  2:415614/1‑97 (MQ=255)
                                                                  gtcgtTGTACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCTTTTTTCTTAAGGTTTTTTCCtt                                 <  1:722817/100‑1 (MQ=255)
                                                                   tcgtTGTACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCTTTTTTCTTAAGGTTTTTTCCTTc                                <  1:988521/100‑1 (MQ=255)
                                                                    cgtTGTACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCTTTTTTCTTAAGGTTTTTTCCTTCg                               >  1:524269/1‑100 (MQ=255)
                                                                     gtTGTACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCTTTTTTCTTAAGGTTTTTTCCTTCGc                              <  1:1591663/100‑1 (MQ=255)
                                                                       tGTACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCTGCCATAATAACCTCATGGTTTTTATGGCTTTTTTCTTAAGGTTTTTTCCTTCGCTg                            <  2:2000730/100‑1 (MQ=255)
                                                                       tGTACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCGGCCATAATAACCTCATGGATTTTATGGCTTTTTTCTTAAGGTTTTTTCCTTCGCTg                            >  2:859992/1‑100 (MQ=255)
                                                                        gTACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCTTTTTTCTTAAGGTTTTTTCCTTCTCTGa                           >  2:209783/1‑100 (MQ=255)
                                                                         tACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCTTTTTTCTTAAGGTTTTTTCCTTCGCTGan                          <  2:2229907/100‑2 (MQ=255)
                                                                           cATTTAACAAACTGCTCTTACTCCTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCTTTTTTCTTAAGGTTTTTTCCTTCGCTGAata                        <  2:398713/100‑1 (MQ=255)
                                                                              ttAACAAACTGCTCTTACTCTTTCCCTTCTCAAATCTGCCATAATAACCTCATGGTTTTTATGGCTTTTTTCTTAAGTTTTTTTCCTTCGCTGAATATCt                     <  2:811846/100‑1 (MQ=255)
                                                                               tAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCTTTTTTCTTAAGGTTTTTTCCTTCGCTGAATATCTa                    <  2:925374/100‑1 (MQ=255)
                                                                                aaCAAACTGCTCTTACTCCTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCTTTTTTATTAAGGTTTTTTCCTTCGCTGAATATCTaa                   >  2:446777/1‑100 (MQ=255)
                                                                                    aaCTGCTCTTACTCTTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCTTTTTTCTTAAGGTTTTTTCCTTCGCTGAATATCTAattat               <  2:1563021/100‑1 (MQ=255)
                                                                                      cTGCTCTTACTCTGTCTCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCTTTTTTCTTATGGTTTTTTCCTTCGCTGAATATCTAAttatta             <  2:391656/100‑1 (MQ=255)
                                                                                       tGCTCTTACTCCTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCTTTTTTCTTAAGGTTTTTTCCTTCGCTGAATATCTAATTATTAc            <  1:1570165/100‑1 (MQ=255)
                                                                                       tGCTCTTACTCCTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCTTTTTTCTTAAGGTTTTTTCCTTCGCTGAATATCTAATTATTAc            <  1:148521/100‑1 (MQ=255)
                                                                                         ctctTACTCTTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCTTTTTTCTTAAGGTTTTTTCCTTCGCTGAATATCTAATTATTACCg          <  1:160884/100‑1 (MQ=255)
                                                                                          tctTACTCCTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCTTTTTTCTTAAGGTTTTTTCCTTCGCTGAATATCTAATTTTTACCGt         >  1:1212753/1‑100 (MQ=255)
                                                                                             tACTCTTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCTTTTTTCTTAAGGTTTTTTCATTCGCTGAATATCTACTTATTACCGTCa       >  2:1905819/1‑99 (MQ=255)
                                                                                             tACTCCTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCTTTTTTCTTAAGGTTTTTTCCTTCGCTGAATATCTAATTATTACCGTCAg      >  1:1081506/1‑100 (MQ=255)
                                                                                             tACTCCTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTGTATGGCTGTTTTCTTAAGGTTTTTTCCTTCGCTGAATATCTAATTATTACCGTCAg      >  1:1219020/1‑100 (MQ=255)
                                                                                              aCTCTTTCCCTCCTCAAATCGGCCATAATAACCTCATGGATTTTATGGCTTTTCTCTTAAGGTTTTTTCCTTCTCTGAATATCTAATTATTACCGTCAGa     >  1:2230598/1‑100 (MQ=255)
                                                                                                 ctTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCTTTTTTCTTAAGGTTTTTTCCTTCGCTGAATATCTAATTATTACCGTCAGAAAt  >  1:416417/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
CATGCGCGCACCGGAGCAGTACCTCTGGATCCACCGTCGCTTTAAAACGCGCCCGGTGGGAGAATCGTCGTTGTACATTTAACAAACTGCTCTTACTCTTTCCCTCCTCAAATCGGCCATAATAACCTCATGGTTTTTATGGCTTTTTTCTTAAGGTTTTTTCCTTCGCTGAATATCTAATTATTACCGTCAGAAAT  >  NC_000913/2496484‑2496680

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